CC BY
300
УДК 621.486
ОПРЕДЕЛЕНИЕ ОБЛИКА ТЕРМОАКУСТИЧЕСКОГО ГЕНЕРАТОРА ПО ХАРАКТЕРИСТИКАМ ЭЛЕКТРОДИНАМИЧЕСКОЙ ЧАСТИ НА ОСНОВАНИИ ПАРАМЕТРИЧЕСКИХ ИССЛЕДОВАНИЙ
© 2011 Г. В. Воротников, А. Н. Крючков, Е. А. Зиновьев
Самарский государственный аэрокосмический университет имени академика С.П. Королёва
(национальный исследовательский университет)
Рассмотрен термоакустический генератор, состоящий из термоакустического двигателя на основе бегущей волны и линейного электрогенератора 18102МЛ. С целью оптимизации конструкции проведен параметрический анализ методом электроакустических аналогий. Получен ряд параметрических зависимостей эффективности установки для качественного анализа поведения системы.
Термоакустика, термоакустический двигатель, термоакустический электрогенератор, режимнопараметрический анализ, электромеханическая аналогия.
Термоакустические энергопреобра-
зующие устройства в последние годы испытывают возрастающий интерес со стороны энергетических компаний крупных развитых стран. Этот интерес вполне понятен. В режиме теплового двигателя термоакустические системы способны трансформировать тепловую энергию любого происхождения в электрическую. При этом они имеют минимум подвижных механических частей (или вообще их не имеют), что повышает надёжность системы, увеличивает срок эксплуатации и снижает затраты на изготовление.
Двигатели на основе бегущей волны теоретически считаются самыми эффективными из термоакустических двигателей, что доказано в работах Свифта [1,2]. Так, для формирования необходимого импеданса в системе на основании пульсационной трубы с жиклёром, получившей широкое распространение из-за простоты, вводится элемент с большой акустической резистивностью (кран, жиклёр), на котором рассеивается значительная акустическая мощность. В то же время в системе на основе бегущей волны нужный импеданс системы определяется в основном подбором реактивных составляющих (податливости, инерционности). При этом неиспользуемая нагрузкой часть мощности возвращается по инерционному каналу обратной связи к регенератору, что позволяет повысить КПД системы.
Вторым неоспоримым преимуществом системы на основании бегущей волны является тот факт, что сдвиг фаз между осциллирующим давлением и объёмной скоростью
осцилляции близок к нулю. Это повышает обратимость термодинамического процесса энергопреобразования, протекающего в регенераторе [1-3].
Большинство публикаций, посвящённых термоакустическим двигателям на основе бегущей волны, рассматривают работу такого двигателя с некой абстрактной акустической нагрузкой или вовсе без неё, оценивая эффективность двигателя как отношение мощности, поступающей в резонатор с нагрузкой, к тепловой мощности, подводимой через горячий теплообменник. Оцененная таким образом эффективность может достигать 40% и более. Такой подход, по нашему мнению, не отражает полной картины потерь акустической мощности, так как прежде всего игнорирует фазовые требования самой акустической нагрузки.
В данной работе приводится качественная оценка потерь в термоакустическом двигателе на основе бегущей волны с присоединённым электрическим генератором в качестве нагрузки. Диапазон геометрических размеров установки выбирался из соображений уменьшения массово-габаритных показателей. Анализ выполнен с использованием электромеханических аналогий.
Установка
Термоакустический генератор (рис.1) состоит из двух основных узлов: термоакустического двигателя и альтернатора. Термоакустический двигатель представляет собой замкнутый канал круглого сечения, в котором последовательно размещены: холодный
теплообменник, регенератор, горячий теплообменник и буферная трубка с вспомогательным холодным теплообменником. Петля замыкается узкой “инерционной трубкой”, которая постепенно переходит в широкую трубку с большой акустической податливостью. В основу данной конструкции положена схема термоакустического двигателя на основе бегущей волны ТЛБИЕ, предложенная Свифтом [1]. Такая схема является гибридной, поскольку акустическая нагрузка располагается не в двигателе, а в отдельно стоящем резонаторе. Однако она получила широкое использование благодаря существенно меньшим габаритам, чем чистые системы на основе бегущей волны. Согласно статье Свифта автору удалось получить на своей установке максимальную эффективность 30% при температурном перепаде на регенераторе в 700 К. Эффективность оценивалась как отношение акустической мощности, рассеянной в резонаторе, к подведённой тепловой мощности.
Рис. 1. Схема термоакустического электрогенератора: 1 - регенератор; 2 - холодный теплообменник;
3 - горячий теплообменник; 4 - термическая буферная трубка; 5 - вспомогательный холодный теплообменник; 6 - струйные диоды, 7 - инерционная трубка; 8 - податливость, 9 - альтернатор
Термодинамический узел двигателя, включающий теплообменники и регенератор, определялся на основании установки Свифта, как описано в [5]. При этом площади поперечных сечений двигателя первоначально были уменьшены пропорционально ожидаемой акустической энергии.
Альтернатор представляет собой линейный электродинамический генератор переменного тока, который снабжён поршнем и размещён в специальном корпусе. В качестве альтернатора был выбран 1S102MA производства CFIC Inc (табл.1). Зазор между поршнем и стенкой корпуса в нём составляет 15 мкм.
Таблица 1. Характеристики 1S102MA
Площадь поршня A 2.03e-3 2 м
Предельная амплитуда перемещения поршня f 5.e-3 м
Подвижная масса M 0.5 кг
Механическое трение Rm 3.5 Нс/м
Жёсткость подвески K 4.2e+4 Н/м
Внутренняя электроре-зистивность Re 7 Ом
Индуктивность катушки L 0.11 Гн
Коэффициент электромеханической связи BL 60 Тл-м
Прибегая к определённым упрощениям и используя принципы электроакустической аналогии [4], термоакустический электрогенератор можно представить в виде схемы сосредоточенных элементов (рис. 2). Здесь символом Я обозначены резистивные элементы, С - податливости, Ь - инерционности, 2 - сложные импедансы. Таким образом, термоакустический двигатель на схеме представлен инерционностью ЬА инерционной трубки 7 (рис. 1), податливостью СА широкого участка 8, податливостью Сы буферной трубки 4, а также резистивностью Яд и источником для теплообменного блока 1-3. Акустическая податливость и инерционность определяются следующими выражениями:
V
С=—, (1)
L =
gP
m
PmL
S ''
(2)
где рт и рт -давление и плотность газа в отсутствие колебаний, у - показатель адиабаты,
V - объём податливости, I и £ - длина и площадь сечения трубки.
В качестве объекта анализа рассматривался двигатель с параметрами, приведёнными в табл. 2, кроме оговорённых особо случаев.
Альтернатор 9 представлен акустическим импедансом 2^> электродинамика, сопротивлением щели Яв между поршнем и
корпусом, а также податливостью Св объёма, расположенного за поршнем.
Таблица 2. Характеристики термоакустического
двигателя (рис. 1
Площадь поперечного сечения трубки 7 3.14е-4 2 м
Приведённая длина трубки 7 0.41 м
Площадь поперечного сечения трубки 8 2.3е-3 2 м
Приведённая длина трубки 8 0.2 м
Площадь поперечного сечения трубки 4 1.52е-3 2 м
Длина трубки 4 0.15 м
Среднее давление 3.0е+6 Па
Резистивность теплообменного узла с диодами 9.е+6 Па- с/м3
Электродинамический преобразователь представляет собой четырёхполюсник, к клеммам которого на входе приложен перепад давления Рвх - Рвых (разность давлений перед поршнем и за поршнем). При этом во входной цепи “течёт” объёмная скорость и, равная произведению скорости поршня электродинамика на его площадь.
Рх - Рых = 2,,и + В-1. (3)
В выходном контуре течёт индукционный электрический ток I, а напряжение на выходных клеммах составляет V:
V = 2, -. (4)
■=Р I 1
\Х
<-Л
Рис. 2. Схематическое представление термоакустического электрогенератора в сосредоточенных электрических аналогах
Отсюда акустический импеданс электродинамического преобразователя с внешней электрической нагрузкой 2н определяется выражением
1 Б!2
7 = 7 +-
° ш А 7е + 7Н'
Механическая составляющая
7 = —
ш А2
(5)
(6)
где Ят - механическое трение, М - подвижная масса, К - жёсткость подвески электродинамика, А - площадь поршня.
Собственный электрический импеданс электродинамического преобразователя складывается из резистивного Яе и индуктивного сопротивления катушки Ь:
2е = Яе + 1®Ь • (7)
Величины А, Ят, М, К, Яе, Ь также являются паспортными данными электродинамика.
Выражение (5) позволяет рассчитать акустический импеданс альтернатора по известной электрической нагрузке в выходной цепи.
Оценка импедансов двигателя и альтернатора
В работе [4] мы показали, что акустический импеданс на выходе из двигателя данной конфигурации определяется только параметрами двигателя и описывается выражением (8):
(
1
1
1
Л
7„ =
7ь, 7ы
7ка 7ы
7Ь 7С
7„ 7С
Ка С а
7са 7ы
где Ъ с соответствующим индексом - импе-дансы элементов акустической цепи, приведённой на рис. 2, &»ТТ/Тх - отношение температур газа на горячем и холодном теплообменниках.
Если устройство не имеет никаких промежуточных систем между двигателем и альтернатором, то импеданс на выходе двигателя 7в должен быть равен импедансу на входе в альтернатор 7ац, который из условия соединения составляющих своей цепи находится, как
7 7
7 = ^яв^п + 7
аЫ 7 + 7 ^^Св’
Кв О
(9)
где 7Кв - импеданс щели между поршнем и корпусом альтернатора, 1Св - импеданс полости за поршнем.
Так как узел регенератора и буферная трубка имеют фиксированные конструктивные параметры [5], изменять импеданс сис-
темы можно лишь двумя локальными переменными 7Т и 7С .
ТА СА
На рис. 3, 4 приведены результаты расчёта составляющих импеданса двигателя для различных значений ТА, а на рис 5, 6 - для различных значений СА. На этих же рисунках для сравнения приведены значения составляющих импеданса альтернатора. При этом в электрическую цепь альтернатора была включена активная нагрузка в 110 Ом. Диапазон изменения податливости СА охватывает диапазон допустимых объёмов от 0.1 до 3л. Диапазон изменения инерционности ТА тоже охватывает весь конструктивно приемлемый интервал изменения податливости. Резонансная частота альтернатора при этом составляет примерно 62 Гц.
Рис. 3. Вещественные составляющие импедансов альтернатора (при Я Н =110 Ом) и двигателя. Импеданс двигателя приведён для различных значений инерционности ЬА
Рис. 4. Мнимые составляющие импедансов альтернатора (при К Н =110 Ом) и двигателя. Импеданс двигателя приведён для различных значений инерционности
Та
Рис. 5. Вещественные составляющие импедансов альтернатора (при К Н =110 Ом) и двигателя. Импеданс двигателя приведён для различных значений податливости СА
Как видно из приведённых графиков, импеданс двигателя существенно отличается от импеданса альтернатора во всех рассматриваемых диапазонах изменения параметров
Рис. 6. Мнимые составляющие импедансов альтернатора (при Я Н =110 Ом) и двигателя. Импеданс двигателя приведён для различных значений податливости СА
двигателя. Чтобы приравнять их, необходимо изменить либо выходной импеданс двигателя, либо входной импеданс нагрузки. Импеданс двигателя проще всего изменяется
введением резистивности на входе в теплообменный узел. Импеданс в ветви нагрузки можно изменить введением параллельной нагрузки 7Н. И тот и другой способы требуют введения существенной резистивной составляющей. Все другие способы являются комбинациями этих двух вариантов. Для своей установки мы избрали второй вариант (рис. 7).
і ~г иг СА 71
рс ї
"ІА (
1-А Ъи
- На
г
(0 -1 і иНА І'-
Т Iі
Рис. 7. Схема термоакустического электрогенератора с параллельной системой настройки
Для выбранной схемы эффективность использования нагрузки в контуре В резонатора с нагрузкой определяется выражением
Ев
Ев
Яе 1/ У7В _
Яе 1/ ./7в _
(10)
Как следует из рис. 3 и 5, резистивность системы настройки сравнима с резистивно-стью полезной нагрузки, а значит и необходимые потери мощности будут соизмеримы, так как они пропорциональны резистивной составляющей импеданса. Другими словами, почти половина акустической мощности, входящей в ветвь В нагрузки, должна быть рассеяна в системе настройки. Легко убедиться, что наилучшее значение соотношения (8) в диапазоне 80-90 Гц составляет примерно 0.65.
Определение частоты рабочей точки
На рис. 8-9 приведены графики зависимости внутреннего КПД системы от частоты и характеристик податливости и инерционности при 0 =2,7. Области, в которых вещественная часть импеданса настройки отрицательна, на графиках отсутствуют. Графики построены на основании выражения
1
\7о7кв
0
(7э + 7Кв )■
' 7в7кв
Яе
Яе
-7Г
(11)
7с (7, + 7К ) + 7, 7К
СА \ ЬА КЛ / ЬА КЛ
полученного в [4]. Из этого же источника взята формула для оценки акустической мощности, поглощённой альтернатором
%с (7, + 7„ ) + 7, 7„
СА \ ^А ^А / ^А К
ев = {I Рс\2
7Сл (7ЬЛ + гял )
(12)
7 7
7Св (7В + 7Кв ) + 7В7К
Яе
из которой найденное Ео имеет гораздо большую абсолютную величину, что делает её предпочтительнее для использования. На основании этого выбираем рабочую частоту 80 Гц.
Оценка максимального смещения
Существенным фактором, ограничивающим использование установки вблизи резонанса, может стать ограниченность хода поршня, который в данном случае составляет X = 5мм. То есть должно выполняться неравенство, которое следует из связи амплитуды смещения и скорости для гармонического колебания
р0\<ХсоА . (13)
Отсюда получается выражение для максимальной энергии, которая в принципе может быть поглощена альтернатором с нагрузкой без превышения ограничений по амплитуде:
Еп„, = 2 2В ]. (14)
Из рис. 12, полученного на основании выражения (14) и (12), можно заключить, что при / = 80 Гц амплитуда давления Рс» 1.5е+5 Па, а из рис. 11, полученного из уравнения (12), что Ео »100 Вт.
Эффективность альтернатора
Для данного альтернатора с фиксированной нагрузкой эта величина определяется только рабочей частотой системы. Эта формула объясняет такие расхождения мощности в дорезонансной и зарезонансной областях на рис. 11.
X
2
X
2
1 1 1 1_а=1.20е+3 1_д=3.42е-3
ч І 1 1_д=8.5 1_д=2.5 *0-3 6е-4
і
//
г • у
о •л 0 7 0 9 0 1] ю і: г-1 О >0 1 ГО 19
Рис. 8. Влияние податливости СА на внутренний КПД системы
I Гц
Рис. 9. Влияние инерционности на внутренний КПД системы
I Гц
Рис. 10. Мощность, поглощённая альтернатором, как функция частоты и амплитуды давления Рс на входе в теплообменный блок
I
Э
а
о
с
го
5
I-
5
С
С
(С
Рис
80
180
130
1, Гц
11. Амплитуда поршня альтернатора как функция частоты и амплитуды давления Рс на входе в теплообменный блок
■ 120 Вт 100 Вт
■ 80 Вт 50 Вт
80
-|-
180
Рис. 12. Внутренний КПД электрического генератора как функция электрической нагрузки при максимальном ходе поршня и различной акустической мощности, рассеиваемой на генераторе
До сих пор мы рассматривали акустическую мощность, поглощаемую альтернатором, как акустическую нагрузку. Однако на практике важна лишь та её часть, которая преобразуется в полезную электрическую мощность. Внутренний КПД электрического генератора
130
Рн. Ом
Рис. 13. Электрическая мощность генератора как функция электрической нагрузки при максимальном ходе поршня и различной акустической мощности, рассеиваемой на генераторе
(15)
Е± = И2 Ке [У ^ ]
Ев |^|2Яе [] ‘
Выражая V из (4) и подставляя в (15), получаем
Яе [ 7я ]
ЕГ
\2я + г,! Яе[]
Максимальная акустическая нагрузка, рассеиваемая на генераторе определяется выражением (14).
На рис. 12 приведены графики внутреннего КПД электрического генератора как функции электрической нагрузки Ян для различной акустической мощности. Необходимо обратить внимание, что эти характеристики получены независимо от возможностей термоакустического двигателя (который может не обеспечить необходимой мощности), а основаны лишь на возможностях конструктивных параметров альтернатора.
100 Вт
80 Вт 50 Вт
Ф * 9 Л *
Л * Ф •* фф
30 80 130 180
Рн. Ом
Рис. 14. Частота электрогенератора как функция электрической нагрузки при максимальном ходе поршня и различной акустической мощности, рассеиваемой на генераторе
Из графиков рис. 14 следует, что при поглощаемой акустической мощность в 100 Вт частота 80 Гц поддерживается при электрической нагрузке в 90 Ом.
Из графиков рис. 13 следует, что максимальная электрическая мощность для акустической нагрузки в 100 Вт при нагрузке в 90 Ом составляет 82 Вт.
Выводы
На основании рассмотренных потерь, вследствие согласования импедансов двигателя и нагрузки, а также внутренних потерь в альтернаторе, можно заключить, что для данной конструкции установки максимальное преобразование акустической энергии, выходящей от идеализированного тороидального термоакустического двигателя, в реальной установке при мощности, близкой к номиналу, генератора составит не более
0,65x0,82 » 0,53%. Поэтому, если в идеализированном двигателе отношение мощности, поступающей в резонатор с нагрузкой, к мощности, выходящей из регенератора составляет 40% (что соответствует нашему пе-
репаду температур), то реальный КПД установки составит примерно 21 %.
Большие расхождения в импедансах двигателя и альтернатора с нагрузкой приводят к тому, что в систему вводятся дополнительные узлы, на которых будет рассеваться значительная акустическая мощность, что снижает внутренний КПД системы почти в 2 раза. Такая схема с такой акустической нагрузкой делает акустически совершенный двигатель, подобный использованному в TASHE [2], бессмысленным, так как всё равно приходится увеличивать необходимое резистивное сопротивление. В этой ситуации более оправдана технологически упрощённая модель, подобная установке Northrop Grumman [6].
В целом же напрашивается принципиальный вывод о том, что для повышения эффективности системы двигатель и альтернатор должны составлять единый контур, не содержащий разветвления энергетических потоков (обратной связи и нагрузки), в котором необходимый импеданс настраивается за счёт инерционной обратной связи, а часть энергии, которая не поглощена альтернатором, будет возвращаться в термодинамический узел. Такая конфигурация установки будет рассмотрена нами в дальнейших работах.
Ожидаемый коэффициент преобразования энергии в 20% довольно высок и соответствует КПД лучших машин Стирлинга. Однако ресурс термоакустических систем существенно выше и при использовании современных материалов для подвески якоря и горячих частей теплового узла превышает 100 000 часов против 30 000 часов у двигателей Стирлинга. Приминение таких машин в современной энергетике сулит большие перспективы и потому требует дальнейших исследований.
Библиографический список
1. Backhaus, S. A thermoacoustic-Stirling heat engine: detailed study [Text] / S. Backhaus, G.W. Swift // J. Acoust. Soc. Am., 2000. - P. 3148-3166.
2. Swift, G.W. Acoustic recovery of lost power in pulse tube refrigerators [Text] / G.W. Swift, D.L. Gardner, S. Backhaus // J. Acoust. Soc. Am., 1999. - P. 711-724.
3. Swift, G. W. Thermoacoustics: A unifying perspective for some engines and refrigerators [Text] / G.W. Swift // Acoust. Soc. Am. 2002.
4. Крючков, А.Н. Анализ термоакустического электрогенератора методом электрической аналогии [Текст] / А.Н. Крючков, Г.В. Воротников, Е.А. Зиновьев // Вестн. СГАУ, 2009. - С. 48-51.
5. Довгялло, А.И. Режимно-параметрический анализ работы термоакустического двигателя с внешним подводом тепла [Текст]
/ А.И. Довгялло, Е.А. Зиновьев, Г.В. Воротников // Вестн. СГАУ. 2011. - С. 146-154.
6. Telesz, M.P. Design and testing of a thermoacoustic power converter [Text] / М.Р. Telesz // A Thesis Presented to The Academic Faculty. Georgia Institute of Technology. 2006.
THE DETERMINATION OF THE THERMOACOUSTIC GENERATOR ASPECT FROM ELECTRODYNAMIC UNIT FEATURES BY THE PARAMETRIC-STUDIES
© 2011 G. V. Vorotnikov, A. N. Kruchkov, E. A. Zynovyev
Samara State Aerospace University named after academician S.P. Korolyov (National Research University)
Thermoacoustic generator, consisting of a thermoacoustic travelling-wave engine and linear alternator 1S102MA, is being considered. Parametric analysis is carried out by method electroacoustic analogy for the design optimization. Authors received different device-performance dependencies for analysis of apparatus behavior.
Thermoacoustics, thermoacoustic heat engine, thermoacoustic alternator, parametric operation analysis, acoustic power, electromechanical analogy.
Информация об авторах
Крючков Александр Николаевич, доктор технических наук, доцент кафедры автоматизированных систем энергетических установок, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). Тел.: (846) 335-18-12. E-mail: [email protected]. Область научных интересов: виброакустика машин.
Воротников Геннадий Викторович, аспирант кафедры автоматизированных систем энергетических установок, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). Тел.: (846) 228-50-98. E-mail: [email protected]. Область научных интересов: рабочие процессы термоакустических двигателей и холодильных машин.
Зиновьев Евгений Александрович, аспирант кафедры теплотехники и тепловых двигателей, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). Тел.: (846) 228-50-98. E-mail: [email protected]. Область научных интересов: рабочие процессы термоакустических двигателей и холодильных машин.
Kruchkov Aleksandr Nikolaevich, doctor of technical science, associate professor of department automation technology of power plants, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). Phone: (846) 335-18-12. E-mail: [email protected]. Area of research: the machine vibroacoustics.
Vorotnikov Gennadiy Viktorovich, post-graduate student of department automation technology of power plants, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). Phone: (846) 228-50-98. E-mail: [email protected]. Area of research: work processes of thermoacoustic heat engines and refrigerators, power efficient heat engines and prime movers.
Zynovyev Evgeniy Aleksandrovich, post-graduate student of department heat engineering and heat engines, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). Phone: (846) 228-50-98. E-mail: [email protected]. Area of research: work processes of thermoacoustic heat engines and refrigerators.
