Файл: Численное моделирование процессов энергоразделения в потоках сжимаемого газа.pdf
Добавлен: 12.12.2023
Просмотров: 174
Скачиваний: 3
ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.
13
Достоверность полученных результатов обеспечивается применением со
временных методов математического моделирования, а также подтверждается использованием процедур верификации и валидации. Результаты находятся в соответствии с результатами, полученными другими авторами.
Апробация работы. Основные результаты работы докладывались на:
международной конференции «Современные проблемы теплофизики и энерге
тики» (Москва, 2017), XXIII и XXIV Международных конференциях «Нели
нейные задачи теории гидродинамической устойчивости и турбулентность
Не-За-Те-Ги-Ус» (Звенигород, 2018, 2020); VII Российской национальной кон
ференции по теплообмену РНКТ (Москва, 2018); XII Всероссийском съезде по фундаментальным проблемам теоретической и прикладной механики (Уфа,
2019); XIV, XIX, XX, XXI и XXII школах-семинарах под руководством акаде
мика Леонтьева А.И. (Санкт-Петербург, 2007; Орехово-Зуево, 2013; Звенигород,
2015; Санкт-Петербург, 2017; Москва, 2019); Всероссийских конференциях моло
дых учёных-механиков (Сочи, 2017, 2018); «Ломоносовские чтения» МГУ имени
М.В. Ломоносова (Москва, 2017–2021); 5th International Workshop on Heat/Mass
Transfer Advances for Energy Conservation and Pollution Control (Новосибирск,
2019); V Всероссийской научной конференции «Теплофизика и физическая гидродинамика» (Крым, Ялта, 2020); Двадцать первой международной школе
семинаре «Модели и методы аэродинамики» (Крым, Евпатория, 2021).
Личный вклад. Диссертационная работа выполнена в лаборатории ги
перзвуковой аэродинамики НИИ механики МГУ им. М.В. Ломоносова. Автор принимал участие в разработке общей концепции и методики проведения численных исследований. Автором разработан способ моделирования, позво
ляющий учитывать взаимодействие между потоком и проницаемой стенкой при наличии вдува/отсоса. Автором предложено использовать испарительное охлаждение для повышения давления торможения в сверхзвуковом канале устройства газодинамического энергоразделения. Автором проведены числен
ные исследования течений в устройстве газодинамического энергоразделения и канале с проницаемыми стенками. Автором получены результаты влияния основных факторов на величину энергоразделения.
Публикации. Основные результаты по теме диссертации изложены в 24
печатных изданиях, 2 из которых изданы в журналах, рекомендованных ВАК,
3 в журналах, индексируемых в Scopus, 19 — в тезисах докладов.
14
Объем и структура работы. Диссертация состоит из введения, трёх глав, заключения и трёх приложений. Полный объём диссертации составляет
170 страниц, включая 97 рисунков и 10 таблиц. Список литературы содержит
130 наименований.
Благодарности. Автор выражает благодарность научному руководите
лю — академику РАН, д. т. н., профессору Александру Ивановичу Леонтьеву за долготерпение, постоянное и активное внимание и участие; сотрудникам лабо
ратории №108 НИИ механики МГУ: к. т. н., с. н. с. А. Г. Здитовцу, к. т. н., в. н. с.
Ю. А. Виноградову, с. н. с. М. М. Стронгину , к. т. н., с. н. с. Н. А. Киселёву за всестороннюю поддержку и внимание к работе; к. ф.-м. н. С. А. Чернышеву за ценные замечания и консультации; к. т. н., вед. инж. И. Н. Байбузенко за всесто
роннюю помощь, содействие, ценные советы и замечания.
15
Глава 1. Методы энергоразделения
История методов энергоразделения насчитывает не один десяток лет. Од
ним из первых «зафиксированных» методов был, так называемый, вихревой метод, обнаруженный Ранком в 1931 г. В работе [
13
] Ранк сообщил о самопроиз
вольном разделении вихревого газового потока на холодную и горячую области.
Позднее (1942 г.) Эккерт и Вaйзе [
14
] экспериментально зафиксировали значи
тельную неравномерность в распределении температуры газа на поверхности теплоизолированного цилиндра, поперечно обтекаемого газовым потоком. Да
лее Шпренгер в 1954 г. [
15
] обнаружил существенное повышение температуры газа вблизи торца глухого отверстия при направлении в него недорасширенной газовой струи. В 1986 г. Голдштейн [
2
] сообщил о перераспределении темпера
туры торможения в затопленной недорасширенной струе.
Все приведённые выше случаи представляют собой перераспределение полной энергии потока без совершения им технической работы и теплообме
на с окружающей средой, в результате которого в потоке возникают горячие и холодные области. Подобные явления принято называть энергоразделением
(температурным разделением). На основе некоторых из них созданы устройства для безмашинного температурного разделения потоков [
16
;
17
]. Безмашинного,
так как в этом случае исключается использование машин т.е. получение горяче
го и холодного потоков осуществляется в устройстве, не имеющем подвижных частей, приводимых в движение газом. Иными словами в подобном устройстве газ не совершает техническую работу и не участвует в теплообмене с окру
жающей средой.
1.1 Физические основы энергоразделения
Для анализа механизмов, приводящих к энергоразделению, по аналогии с работой [
18
], рассмотрим уравнение сохранений полной энергии (энтальпии торможения) газового потока:
16
ρ
????ℎ
˚
????????
“
B????
B????
`
B
B????
????
ˆ
λ
B????
B????
????
˙
`
B
B????
????
p????
????
σ
????????
q ,
(1.1)
Первое слагаемое в левой части выражает субстанциональную производ
ную (
????
{
????????
“
B
{
B????
` ????
????
B
{
B????
????
) энтальпии торможения ℎ
˚
потока, которая является суммой статической энтальпии ℎ и кинетической энергии ????
2
{2
, как следует из формулы (
1.2
).
ℎ
˚
“ ℎ `
????
2 2
(1.2)
Субстанциональная производная описывает изменение энтальпии тормо
жения движущейся частицы жидкости во времени. Путь, проходимый частицей для нестационарного потока называется траекторией. Таким образом,
????ℎ
˚
????????
опи
сывает изменение энтальпии торможения частицы жидкости при её движении вдоль траектории, которая для нестационарного потока отличается от линии тока, касательной векторов скорости в фиксированный момент времени.
Первый член в правой части уравнения (
1.1
) выражает частную производ
ную давления ???? потока по времени ????. Два остальных слагаемых в правой части уравнения (
1.1
) выражают перенос энергии теплопроводностью и работой, вы
полненной вязкими напряжениями σ
????????
Рассмотрим три различных случая:
– Стационарное течение при отсутствии трения и теплообмена. В этом случае уравнение (
1.1
) упрощается до:
????ℎ
˚
????????
“ ????
????
Bℎ
˚
B????
????
(1.3)
Для совершенного газа энтальпия является функцией только темпера
туры и в этом случае, как следует из уравнения (
1.3
), температура торможения ????
˚
остаётся постоянной вдоль линии тока, которая в ста
ционарном случае совпадает с траекторией.
ℎ “ ???? p???? q,
????????
˚
????????
“ 0,
????
˚
“ ????????????????????.
Физически тот факт, что температура торможения остаётся постоянной,
легко понять из рис.
1.1
, на котором схематически изображена трубка тока. Для идеальной жидкости тепловые потоки ???? и вязкие напряже
ния τ отсутствуют. Давление жидкости ????, нормальное к поверхности трубки, не выполняет работу, так как поверхность не перемещается в
17
направлении ????, и жидкость в трубке не может получать или терять энергию, так как она движется вдоль. Таким образом, температура тор
можения постоянна вдоль линий тока.
– Стационарное течение в пограничном слое. Уравнение энергии (
1.1
) в этом случае будет выглядеть следующим образом:
ρ
????ℎ
˚
????????
“
B
B????
ˆ
λ
B????
B????
˙
`
B
B????
ˆ
????µ
B????
B????
˙
(1.4)
Нормальные вязкие напряжения σ
????????
в этом случае пренебрежимо малы.
Из уравнения (
1.4
) следует, что энтальпия торможения ℎ
˚
может изме
няться вдоль линии тока либо за счет механизма теплопроводности λ,
либо за счет переноса энергии в качестве работы, выполняемой каса
тельными напряжениями µ
B????
{
B????
. Жидкость в трубке тока (см. рис.
1.1
)
теперь может получать или терять энергию за счет теплопроводности ????
или работы, выполняемой касательными напряжениями τ, потому что поверхность трубки тока движется в направлении этих напряжений.
– Нестационарное течение при отсутствии трения и теплообмена. Уравне
ние (
1.1
) в этом случае примет следующий вид:
ρ
????ℎ
˚
????????
“
B????
B????
,
(1.5)
из которого следует, что энтальпия торможения ℎ
˚
изменятся вдоль линии тока, если давление в фиксированной точке пространства изме
няется во времени.
????
τ
????
Рисунок 1.1 — Схема элементарной трубки тока
Таким образом, из уравнения (
1.1
) следует, что процесс энергоразделе
ния (изменение энтальпии торможения) может происходить под действием трёх факторов:
– нестационарных пульсаций давления;
– механизма теплопроводности;
– работы вязких сил.
18
Количественными характеристиками процесса энергоразделения являют
ся разности между среднемассовыми температурами торможения газового потока на «горячем» ????
˚
ℎ
и «холодном» ????
˚
????
выходах и на входе устройства ????
˚
0
:
Δ????
˚
ℎ
“ ????
˚
ℎ
´ ????
˚
0
;
Δ????
˚
????
“ ????
˚
????
´ ????
˚
0
(1.6)
Ниже приведён обзор устройств, реализующих энергоразделение под дей
ствием одного или нескольких указанных выше факторов.
1.2 Вихревая труба Ранка-Хилша
На данный момент наиболее известным и распространенным методом безмашинного энергоразделения является вихревой метод, предложенный Ж.
Ранком в 1931 г. При измерении температуры воздуха в циклонном пылеуло
вителе Ранк заметил, что на оси вихря температура ниже, чем на периферии.
Он провёл первые экспериментальные исследования зафиксированного им эф
фекта [
13
] и предложил устройства для его реализации (рис.
1.2
). В 1933 г.
Ж. Ранк на заседании Французского физического общества сделал доклад об открытом им явлении разделения сжатого газа, в котором пытался развить цен
тробежную теорию вихревого эффекта. Достигнутый в эксперименте эффект энергоразделения (32 К) примерно в 4 раза превышал значения, получаемые из теоретических исследований. Сообщение изобретателя было встречено с недо
верием, опытные данные объявлены ошибкой измерения, и об эффекте Ранка забыли.
Первое тщательное исследование вихревого эффекта было проведено Ро
бертом Хилшем более чем через десятилетие после доклада Ранка. Результаты работы не вызывали сомнений [
19
]. Устройство для получения холодного и го
рячего потоков из сжатого газа стали называть трубой Ранка (позже она стала называться трубой Ранка—Хилша или вихревой трубой). Устройство и принцип действия трубы Ранка—Хилша (ТРХ) очень просты. Она состоит из цилиндри
ческой трубы 1, сопла 2, диафрагмы 3 и вентиля 4 (рис.
1.2
).
В сопле газ приобретает начальную тангенциальную составляющую ско
рости. Расширяясь в трубе, приосевые слои газа охлаждаются и покидают
19
Рисунок 1.2 — Принципиальная схема трубки Ранка—Хилша [
11
].
1 — цилиндрическая труба; 2 — сопло; 3 — диафрагма; 4 — вентиль;
устройство через отверстие диафрагмы, периферийные слои разогреваются и отводятся через вентиль с противоположного конца. Изменяя положение вен
тиля, можно менять соотношение расходов горячего и холодного потоков. При этом их температуры также изменяются.
Однако несмотря на исключительную простоту данного устройства и большое количество накопленной информации, до сих пор нет единой об
щепризнанной теории происходящего в ней процесса энергоразделения. В
работах [
1
;
20
] собраны основные физические модели, описывающие процесс в устройстве, и приведен их критический анализ. Ниже перечислены основные из них:
1. температурное разделение в результате сжатия и расширения, возника
ющего в поле центробежных сил;
2. акустические течения;
3. гипотеза Фултона о перестройке характера вихревого течения;
4. вторичные течения, формирующие холодильный цикл для переноса энергии от осевых слоев к периферийным;
5. градиент статической температуры.
Вероятно, что в той или иной степени все явления, лежащие в основе разных моделей, вносят свой вклад в суммарный эффект. Наиболее наглядной является модель Фултона [
21
], согласно которой поток газа на входе в трубу закручивается по закону свободного вихря (ω???? “ ????????????????????) с большими угловыми скоростями в центре вихря. Статическая температура изменяется по радиусу вихря, увеличиваясь к периферии, а температура торможения распределена равномерно.
В процессе последующей перестройки потока от свободного вихря (с за
коном ω???? “ ????????????????????) к вращению с постоянной угловой скоростью ω “ ????????????????????
20
возникают потоки теплоты, вызываемые градиентом статической температуры и направленные от периферии к центру, и потоки механической энергии от оси к периферии, обусловленные наличием вязкостных сил. В результате происходят выравнивание термодинамической температуры газа по поперечному сечению трубы и увеличение температуры торможения с ростом радиуса.
Обеспечив отвод холодного газа из приосевой зоны и нагретого из при
стеночной зоны, можно получить температурное разделение в вихревой трубе.
Согласно модели Фултона для получения эффекта энергоразделения (сравни
мого с экспериментальными данными) свободный вихрь преобразовывается в вынужденный вихрь с более высокой, чем у свободного вихря, периферийной скоростью.
На рис.
1.3
приведены данные работы [
22
] иллюстрирующие энергораз
деление в вихревой трубе. Нагрев и охлаждение потоков (
1.6
) показаны в зависимости от соотношения массовой доли холодного потока.
1.3 Резонансная труба Гартмана-Шпренгера
В 1916 г. при измерении давления трубкой Пито на оси недорасширенной сверхзвуковой струи Гартман обнаружил, что на некоторых участках возника
ют резкие колебания трубки, сопровождающиеся интенсивным акустическим излучением. Впоследствии выяснилось, что это явление связано с взаимодей
ствием струи и внутренней полости трубки. Расположив на одной оси сопло и трубку, заглушенную с одной стороны, Гартман создал устройство для возбуж
дения звуковых и ультразвуковых волн высокой интенсивности. В последующих работах долгое время внимание экспериментаторов было сосредоточено на ис
следовании волновых явлений, при этом значительных термических эффектов в данном устройстве не наблюдалось. Шпренгер был первым, кто в 1954 г.
обнаружил существенный разогрев задней стенки трубки [
15
]. В отличие от предыдущих исследователей Шпренгер использовал трубку большой относи
тельной длины ????{???? “ 34, тогда как у Гартмана размеры трубки составляли
????{???? ď 4
. Эффект был настолько значительным (температура стенки более чем в 2 раза превышала начальную температуру газа), что в дальнейшем устройство получило название трубы Гартмана—Шпренгера (ТГШ). Её базовая конструк