УДК 621.486
РЕЖИМНО-ПАРАМЕТРИЧЕСКИЙ АНАЛИЗ РАБОТЫ ТЕРМОАКУСТИЧЕСКОГО ДВИГАТЕЛЯ С ВНЕШНИМ ПОДВОДОМ ТЕПЛА
© 2010 А. И. Довгялло, Е. А. Зиновьев, Г. В. Воротников
Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет)
Содержит режимно-параметрический анализ работы регенератора в составе термоакустического двигателя внешнего подвода тепла, основанного на энергопереносе бегущей волной. Представлены рекомендации по определению параметров установки для предпроектной разработки.
Термоакустический двигатель, бегущая волна, регенератор, режимно-параметрический анализ, акустическая мощность, рекомендации по проектированию.
Введение
В настоящее время термоакустические двигатели представляют собой новый класс энергетических систем, где тепловая энергия конвертируется в энергию акустической волны, способную совершать работу.
Такие двигатели обладают большим практическим потенциалом, что проявляется в их сравнительно высокой надёжности, низком уровне стоимости, отсутствии подвижных частей, использовании недорогих материалов и невысоких требованиях по технологичности. Парадигма термоакустики открывает возможность замены широкого класса механических устройств волновыми системами.
Термоакустические двигатели с внешним подводом тепла на сегодняшний день являются новым перспективным направлением в области разработки бортовых систем авиационного и космического назначения, а также энергосберегающих технологий.
В статье представлена оценка основных характеристик термоакустического двигателя на основании анализа работы регенератора с учетом характеристик аналогов.
Двигатель на основе бегущей волны
Питер Хатсон Кеперли в своей работе [1] показал, что газ в регенераторе под воздействием бегущей волны осуществляет термодинамический цикл, подобный циклу Стирлинга. Потенциально такие системы обладают высокой эффективностью.
Схема устройства термоакустического двигателя с внешним подводом тепла показана на рис. 1 [2]. Двигатель представляет собой тороидальный акустический волновод, содержащий регенератор, термическую буферную трубу, горячий и холодный теплообменники. Волновод соединяется с полуволновым резонатором Гельмгольца. В данном устройстве используется регенератор, у которого размеры пор в несколько раз меньше глубины термического и вязкого проникновения. Поэтому для описания термодинамического процесса внутри регенератора можно использовать цикл Стирлинга, требующий, чтобы давление и скорость в регенераторе совпадали по фазе. Кроме того, большие скорости движения осциллирующего газа в порах регенератора вызывают большие потери на вязкое трение. Для получения максимальной акустической мощности в системе необходимо, чтобы на выходе из регенератора давле-
ние и скорость принимали максимальные значения. Поэтому в регенераторе нужно снижать скорость и увеличивать амплитуду колебания давления. Тороидальный волновод пристыковывается к резонатору в точке, где располагается узел скорости (пучность давления).
Рассмотрим регенератор, изготовленный из тонких пластин из нержавеющей стали. По обеим сторонам регенератора располагаются теплообменники. Рабочим газом является гелий.
Рис. 1. Общий вид термоакустического двигателя
(а) - схематическое изображение термоакустического двигателя с внешним подводом тепла;
(б) - тороидальный волновод:
1 - регенератор;
2 - основной холодный теплообменник;
3 - горячий теплообменник;
4 - дополнительный холодный теплообменник;
5 - термическая буферная труба;
6 - акустическая податливость;
7 - инерционная труба;
8 - супрессор;
9 - ответвление к резонатору
Акустическая мощность Жх на холодном конце регенератора возрастает до
значения ЖГ на горячем конце. Разность между Жг и №Х представляет собой акустическую мощность Жр, вырабатываемую регенератором. Тепловая мощность йг подводится к горячему концу регенератора, с холодного конца регенератора отводится мощность йх. Применив закон сохранения энергии, получаем
(1) (2) (3)
Ох = Ж - Н р ,
йг = Нб - Нр ,
Жр = Жг - Жх
где
ЖР - акустическая мощность, вырабатываемая регенератором;
НБ - общая мощность в термической буферной трубе (направлена к ее холодному теплообменнику);
НР - общая мощность в регенераторе (направлена к его холодному теплообменнику).
Введём эффективность регенератора и коэффициент использования акустической нагрузки
Па =
йг
Ж
нагр
(4)
нагр
Эффективность регенератора относительно эффективности цикла Карно записывается в виде
Пр =
Па.
Пк
где Пк =
(5)
Эффективность цикла Карно задаёт верхнюю границу эффективности работы регенератора, так что 0 <ЦР < 1. Поскольку Жр будет больше, чем Жнагр из-за диссипации энергии на стенках узлов системы получаем Пнгр < Па < Пк .
Объектом анализа в данной работе является регенератор как наиболее важный компонент, отвечающий за генерацию акустической волны. Представлены данные о влиянии параметров на эффек-
тивность работы регенератора. Рассматривается двигатель (рис.1) Грегори Свифта и Скотта Бэкхауса [3], в которой в качестве рабочего газа используется гелий со средним давлением 3,1 МПа. Наиболее важные геометрические параметры регенератора - число Лотрека NЛ , пористость ф и длина регенератора Ьр. Наиболее важными рабочими параметрами являются импеданс регенератора 2р , передаточное число Д., равное отношению максимальной амплитуды давления к среднему давлению газа, и перепад температур (ТГ - Тх). На рис. 2 показано влияние параметров на эффективность регенератора. При этом Вг = 5 %, (ТГ - Тх ) = 300 К,
ZР = 30рс. Для оптимальной работы двигателя на основе бегущей волны число Лотрека NЛ должно быть мало, для того чтобы между газом и стенками существовал практически полный термический контакт. Отметим, что число NЛ должно быть меньше единицы, оптимальная вели-
о
чина составляет NЛ ~ 0,275 (при 20 С) (рис. 2, а), что соответствует гидравлическому радиусу 42 мкм [3]. При меньших значениях NЛ диссипация энергии за счет вязкого трения становится значительной, что отрицательно сказывается на работоспособности всей системы в целом.
в)
г)
Рис. 2. График зависимости относительной эффективности регенератора
а) от числа Лотрека NЛ (при 20° С), г) от импеданса 2 р,
б) от пористости ф, д) от передаточного числа д,
в) от длины регенератора ьр, е) от перепада температур (ТГ - Тх ) на концах
регенератора
Для заданного размера регенератора существенное влияние оказывает пористость ф . Из рис. 2, б следует, что с увеличением пористости ф (пластины становятся достаточно тонкими) эффективность регенератора растёт до определённого значения. При дальнейшем увеличении пористости ф эффективность регенератора стремится к нулю, так как наблюдается недостаток теплоёмкости регенератора. Влияние длины регенератора на его эффективность представлено на рис. 2, в. Если регенератор достаточно короткий, то перепад температур будет сравнительно небольшим, следовательно, его эффективность также будет незначительной. Если регенератор слишком длинный, будет иметь место большая диссипация энергии, что снижает его эффективность. Из графика видно, что оптимальная величина длины регенератора для указанных выше условий Ьр = 6,8 см, что довольно близко к значению 7,3 см, указанному в работе [3].
Выше отмечалось, что давление и скорость должны совпадать по фазе, а колебания давления должны быть большими для получения максимальной акустической мощности. Следовательно, необходимо добиваться наиболее оптимального
соотношения скорости и давления путем изменения импеданса. Из рис. 2, г следует, что максимальному значению эффективности соответствует импеданс ZР ~ 16,5рс. На рис. 2, д графически представлено, насколько эффективность регенератора подвержена влиянию передаточного числа. Когда при фиксированном перепаде температур передаточное число увеличивается, мощность на выходе также возрастает, но потери на вязкое трение замедляют темпы её роста. В результате более высокая эффективность наблюдается при увеличении передаточного числа. При дальнейшем увеличении профили эффективности стремятся к КПД по циклу Карно.
Рассматривая воздействие прилагаемого перепада температур, отметим, что он должен быть достаточным для того, чтобы термоакустическое устройство могло работать в качестве двигателя, поскольку при перепаде температур ниже некоторого значения устройство будет работать как холодильная машина, поглощая акустическую энергию. Как следует из рис. 2, е, для данного устройства минимально необходимый перепад температур составляет приблизительно 24 К. Оптимальная эффективность достигается
при перепаде 233 К, а затем начинает снижаться.
При наличии массового потока рабочего газа М Ф 0 эффективность регенератора ПР теоретически должна уменьшаться. Это вызвано тем, что осреднённый по времени поток рабочего газа перемешивает газовые слои с разной температурой, что ухудшает характеристики регенератора. Анализ зависимости относительной эффективности регенератора от осреднён-ного массового потока М / Агаз (рис. 3) показывает, что ухудшение характеристик действительно имеет место. Таким образом, максимальная эффективность достигается при М = 0 и быстро снижается с увеличением абсолютного значения массового потока. Во всех рассмотренных
вариантах установки, описанной в работе [3], в качестве материала регенератора использовалась нержавеющая сталь. На рис. 4 и 5 представлены зависимости относительной эффективности регенератора и произведённой акустической мощности от пористости регенератора ф для различных материалов: нержавеющей стали, корундовой керамики, стекла, майлара (полиэтилен терефталат) и полистирола. Из анализа рис. 4 следует, что многие материалы обеспечивают высокую эффективность регенератора, которая достигается при высокой пористости (тонкие пластины) регенератора. Тем не менее, когда пластины становятся очень тонкими (ф ^ 0 ), эффективность регенератора падает до нуля.
массовый поток |
Рис. 3. Зависимость относительной эффективности регенератора от осреднённого
массового потока М/ А
Чр 0,5
1 • 4 і
1/ •
1 т т ^ 1
✓ у
/ \
алюминий
1 / — • — стекло
1гУ —■— майя ар
—£—полистирол Г 1 | 1 1 1 1 | 1 1 1 1 | 1 II
0 0,2 0,4 0,6 0,-8 Ф 1 1,
Рас. 4. Зависимость относительной эффективности регенератора от пористости ф
для различных материалов
О 0,2 0,4 0,6 0,8 1 1,2
Ф
Рис. 5. Зависимость произведенной акустической мощности
от пористости регенератора ф для различных материалов
Для некоторых материалов, таких, как майлар или пластмассовая соломка, высокая эффективность достигается почти при любой пористости в том случае, если ф не стремится к нулю или единице.
Однако, несмотря на то, что эффективность регенератора может быть достаточно высокой, низкие значения ф приводят к тому, что акустическая мощность почти не производится (рис. 5). Отметим, что для всех материалов наиболее высокой акустической мощности соответствует пористость ф , лежащая в интервале от
0,75 до 0,9.
Из вышесказанного следует, что для практического использования наиболее привлекательными материалами для изготовления регенератора являются стекло и нержавеющая сталь.
Рекомендации
Проектирование термоакустического двигателя - сложная многопараметрическая задача, которая трудно поддаётся оптимизации. Поскольку численные оценки таких устройств ещё несовершенны, их проектирование обычно начинается с рассмотрения аналогов и уже известных зависимостей, а известные алгоритмы оптимизации обеспечивают решение лишь в том случае, если достаточно точно заданы
начальные условия. Исходя из этого, с учётом полученных выше качественных оценок можно сформировать следующие практические рекомендации.
1. Поскольку работа термоакустического устройства опирается на явление резонанса, размер системы вдоль оси распространения волны в первом приближении можно считать пропорциональным частоте акустических колебаний. Отсюда, опираясь на аналог, можно определить продольный размер установки при заданной частоте или рабочую частоту при заданном размере.
2. На основе зависимости, представленной на рис. 6, оценивается КПД установки на основании возможного температурного перепада или необходимый температурный перепад по требуемому КПД. Согласно техническим требованиям по мощности Ж на выходе и КПД системы П определяется количество тепла, подводимого к горячему теплообменнику QГ = Ж / 7]. Распределение тепловой нагрузки между холодными теплообменниками делается в той же пропорции, что и для аналогов.
Рис. 6. График зависимости КПД установки от температурного перепада
3. Так как энергетические характеристики регенератора на низких частотах почти не зависят от материала (рис. 4, 5), материал подбирается на основании рабочего интервала температур в регенераторе и его технологичности. На основе глубины вязкого и термического проникновения вычисляется размер пор и пористость ф регенератора. Значения NЛ и ф определяются на основании рис. 2, 4, 5.
4. Начальные размеры установки в поперечном направлении определяются на основании аналога при изменении размеров пропорционально квадратному корню из соотношения величин энергетических потоков или корню из соотношения комплексов (п • QГ).
5. Так как процедура, указанная в пункте 4, не предусматривает соблюдения условия резонанса, то для его выполнения необходимо откорректировать геометрию резонатора согласно формуле
= Ьа, (6)
аС
где о- частота колебаний рабочего газа, С - акустическая податливость и L -инерционность резонатора.
После того, как будут определены начальные параметры установки, можно приступать к процедуре её оптимизации.
Данные рекомендации могут применяться к любым системам на основе бегущей волны, имеющих с рассматриваемым аналогом одинаковую конструктивную схему и работающих на частотах до 500 Гц при передаточных числах до 10%.
Библиографический список
1. Ceperley, P. H. A pistonless Stirling engine - the travelling wave heat engine [Text]. J. Acoust. Soc. Am. - 1979. - p. 1508-1513.
2. Backhaus, S. Swift, G. W. A thermoacoustic Stirling heat engine [Text]. Nature. - 1999. - p. 335-338.
3. Backhaus, S. Swift, G. W. A thermoacoustic Stirling heat engine: Detailed study [Text]. J. Acoust. Soc. Am. - 2000. - p. 3148-3166.
References
1. Ceperley, P. H. A pistonless Stirling engine - the travelling wave heat engine [Text]. J. Acoust. Soc. Am. - 1979. - p. 1508-1513.
2. Backhaus, S. Swift, G. W. A thermoacoustic Stirling heat engine [Text]. Nature. - 1999. - p. 335-338.
3. Backhaus, S. Swift, G. W. A thermoacoustic Stirling heat engine: Detailed study [Text]. J. Acoust. Soc. Am. - 2000. - p. 3148-3166.
PARAMETRIC ANALYSIS OF THERMOACOUSTIC TRAVELLING-WAVE ENGINE OPERATION WITH EXTERNAL HEAT APPLICATION
© 2010 A. I. Dovgyallo, Ye. A. Zynovyev, G. V. Vorotnikov
Samara State Aerospace University named after academician S. P. Korolyov
(National Research University)
The paper presents parametric operation analysis of a regenerator incorporated into a thermoacoustic travelling-wave engine with external heat application. It focuses on the recommendations for defining the engine parameters in preproject design.
Thermoacoustic heat engine, travelling wave, regenerator, parametric operation analysis, acoustic power, recommendations for design.
Информация об авторах
Довгялло Александр Иванович, доктор технических наук, профессор кафедры теплотехники и тепловых двигателей. Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет). Область научных интересов: рабочие процессы тепловых и холодильных машин, бортовая энергетика, энергосбережение. E-mail: [email protected].
Зиновьев Евгений Александрович, аспирант кафедры теплотехники и тепловых двигателей. Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет). Область научных интересов: рабочие процессы термоакустических двигателей и холодильных машин, энергоэффективные двигатели и движители транспортных средств. E-mail:
Воротников Геннадий Викторович, аспирант кафедры автоматизированных систем энергетических установок. Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет). Область научных интересов: рабочие процессы термоакустических двигателей и холодильных машин. E-mail: vorotnikov [email protected].
Dovgyallo Aleksandr Ivanovich, doctor of technical sciences, professor of the department of heat engineering and heat engines. Samara State Aerospace University named after academician S. P. Korolyov (National Research University). Area of research: work processes of heat engines and refrigerators, airborne on-board electrical devices, energy saving. E-mail: [email protected].
Zynovyev Yevgeny Aleksandrovich, post-graduate student of the department of heat engineering and heat engines. Samara State Aerospace University named after academician S. P. Korolyov (National Research University). Area of research: work processes of thermoacoustic heat engines and refrigerators, power efficient heat engines and prime movers. Email: eazinovev38416@ mail.ru
Vorotnikov Gennadiy Viktorovich, post-graduate student of the department of automation technology of power devices. Samara State Aerospace University named after academician S. P. Korolyov (National Research University). Area of research: work processes of thermoacoustic heat engines and refrigerators, power efficient heat engines and prime movers. E-mail: vorotn [email protected]