CC BY
176
ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2015р. Серія: Технічні науки Вип. 30. Т. 1
ISSN 2225-6733
УДК 533.6.011.72:681.121.089
© Лухтура Ф.И.*
К ВОПРОСУ ОБ УСТАНОВИВШЕМСЯ РЕЖИМЕ ИСТЕЧЕНИЯ ГАЗА ИЗ ОСЕСИММЕТРИЧНЫХ ОТВЕРСТИЙ И СОПЕЛ
Приведен метод расчета действительных параметров газа при истечении из отверстий и сопел с различной конфигурацией сужающейся их части по известным из опыта коэффициентам расхода. Предложен достаточно простой метод расчета величины «второго» критического перепада давлений, зависящей от коэффициентов расхода на запертых режимах истечения. В рамках одномерной модели невязкой жидкости подтверждено, что в выходном сечении отверстия невозможно получить критические скорости.
Ключевые слова: отверстие, сопло, коэффициент расхода, коэффициент скорости, коэффициент сжатия струи.
Лухтура Ф.І. До питання про сталий режим витікання газу з вісосиметричних отворів та сопел. Наведений метод розрахунку реальних параметрів газу при витіканні з отворів і сопел з різноманітними конфігураціями звужувальної їх частини по відомим з опиту коефіцієнтів витрат. Запропоновано достатньо простий метод розрахунку величини «другого» критичного перепаду тиску, який залежить від коефіцієнтів витрати на запертих режимах витікання. В рамках одновимірної моделі нев’язкої рідини підтверджено, що в вихідному перерізу отвору не можна получити критичні швидкості.
Ключові слова: отвір, сопло, коефіцієнт витрати, коефіцієнт швидкості, коефіцієнт стиску струменя.
F.I. Lukhtura. On the question of set gas expiration through axisymmetrical holes and nozzles. Real parameters of the gas expiring out of narrowing nozzles of different configuration have been calculated knowing the discharge coefficients obtained experimentally. Rather a simple method to calculate the «second» critical pressure differential dependent on discharge coefficients at cut-off expiration modes has been offered. Using an inviscid fluid one-dimensional model the author has confirmed that it is impossible to receive critical speeds in the output hole. The method makes it possible to estimate the value of the pressure in the separation point (if it is available), where the Laval nozzle section is minimum and knowing «second» critical pressure differential (or total pressure at cut-off expiration). The received data can be used as initial data to calculate jet flows passing through various kind of holes. The analytical research results of gas discharge regime through the holes (nozzles) have shown, that a one-dimensional model can be used to calculate the parameters of the expiration through holes of different configuration in the nozzle narrowing part precisely enough to estimate changes of the main (basic) flow parameters in the channels of the holes. The analysis of the research indicates that the flow coefficients at flow-off expiration modes for convergent and convergent-divergent holes are identical and the parameters in critical section are identical as well.
The total pressure recovery coefficient is approximately equal to the flow coefficient, there being losses of total pressure. The received results can be used in power engineering, metallurgy etc. developing various blowing devices and burners.
Keywords: hole, nozzle, flow coefficient, coefficient of speed, coefficient of compressibility of a jet.
Постановка проблемы. На протяжении более 170 лет [1], начиная с момента публикаций работ А. Навье (1829 г.) и А.Сен-Венана и Л. Вантцеля (1839 г.), касающихся проблем и
ст. преподаватель, ГВУЗ «Приазовский государственный технический университет», г. Мариуполь, lukhelena@yandex. ru
213
ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2015р. Серія: Технічні науки Вип. 30. Т. 1
ISSN 2225-6733
вопросов истечения газа из отверстий, многие ученые XIX (Ю. Вейсбах, Г. Цейнер, Г.-А. Гирн, Р. Клаузиус, А. Гюгонио, А. Паренти, О. Рейнольдс и др.) и ХХ столетия (С.А. Чаплыгин, Ф.И. Франкль [2, 3] и др.) проявляли достаточно большой интерес к этому процессу, вызывавшему в тот период многочисленные дискуссии на страницах соответствующих изданий и публикаций. С течением времени интерес к данному явлению не снижался и время от времени ученые предпринимали попытки дальнейшего изучения, теоретического и экспериментального исследования [4-26] с целью уточнения механизма истечения газа из отверстий.
Анализ последних исследований и публикаций. Несмотря на длительный период изучения этого явления в современной периодической печати и в пособиях по прикладной газовой динамике [4-26], кроме немногочисленных работ [12, 17, 23, 24], в основном учитывают лишь потери из-за вязкого трения [4, 25, 26 и др.], незначительно изменяющие параметры истечения от теоретически ожидаемых, в большинстве случаев процесс истечения из отверстий и сопл представляют изоэнтропным [8, 26 и многих др.]. Это вызывает некоторое удивление, т.к. во многих практических случаях, как показал исторический опыт изучения, действительный (адиабатный) процесс истечения значительно отличается от теоретического (изоэнтропного). При расчете же параметров истечения несжимаемой жидкости такого недостатка не существует: для определения действительной скорости истечения и расхода несжимаемой жидкости из отверстий теоретическая формула Торричелли используется только лишь с поправочными коэффициентами истечения: коэффициентами сжатия струи є, скорости ф, расхода ц, учитывающие реальность процесса. Кроме того, что касается прикладных аспектов, то, как показал литературный обзор, информация, представленная по истечению газа из отверстий в известных источниках [4-26], достаточно противоречива. Так, например, в одних источниках [4, 21, 24] представлены данные, указывающие на потери скорости (а значит и полного давления) в сужающихся каналах из-за наличия сопротивления, в то же время в других [26] говорится о практически неизменном полном давлении при течении в конфузорах с различным углом наклона его образующей. Наконец, в третьих источниках [8, 9, 22 и др.] о потерях полного давления нет никакой информации. Обычно в работах не приводится анализ причинно-следственной связи между параметрами истечения. А значит, пока отсутствует достаточно полная, соответствующая существующим экспериментальным данным, теория истечения из отверстий различной формы, позволяющая определять действительные параметры на выходе из отверстия и на начальном участке струйного течения: скорость, полное давление - параметры, интересующие, в первую очередь, исследователей и инженерно-технических работников для выполнения прикладных НИР и решения практических задач. К тому же большинство существующих методик расчета [3, 12-14, 17] коэффициента расхода при истечении из отверстий с различными углами наклона образующих в связи со спецификой решения уравнений газовой динамики [27] применимы лишь в области сверхкритических перепадов давлений. Противоречивым, на наш взгляд, является также комментарии по поводу существующего в литературе факта несовпадения, при углах наклона стенок отверстия в диапазоне 65^70°, аналитических [2, 3] (на основе теории Чаплыгина, Франкля), и экспериментальных [9] данных по режимам «запирания» отверстия. Как показано в [23-24], в этом случае не учитываются потери полного давления при истечении газа из отверстий в толстой стенке (сопел) в рассматриваемом диапазоне изменения углов наклона их стенок.
В настоящее время процесс истечения газа из отверстия можно представить следующим образом. С ростом давления перед соплом (или понижением давления в среде, куда истекает поток) расход, следуя выражению Сен-Венана-Вантцеля, растет и теоретически принимает
k
(запертое течение). В
* рх ( 2 1k-1
максимальное значение при критическом перепаде п =--------= I----
KV Ро* \k + 1
действительности из-за неравномерности параметров по сечению (газодинамических потерь) и потерь из-за вязкого трения действительный расход всегда меньше теоретического, а макси-
^ * *
мальное его значение наблюдается только при пкр ниже п , где п и пкр - так называемые
F кр кр F
«первое» и «второе» критические отношения давлений. Таким образом, сопло в действительно-
сти запирается при пкр. Величины п**р
пкр и диапазон давлений пкр
*
< п< п
кр
зависят от формы
214
ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2015р. Серія: Технічні науки Вип. 30. Т. 1
ISSN 2225-6733
отверстия, теплофизических свойств истекающего газа и др. В этом диапазоне с ростом давления перед соплом или по мере увеличения отношения давлений (полного перепада давлений) П0 = 1/п звуковая линия перемещается против потока к срезу отверстия. При этом возрастает коэффициент расхода ц. Следуя наиболее распространенному в настоящее время мнению [2-3, 8-9, 26], изменение расхода происходит вследствие того, что волны давления влияют на линию перехода и дозвуковую часть, т.е. по существу «проходят» внутрь отверстия. Здесь уместно отметить работы А. Паренти (1987 г.) [1], где, в отличие от упомянутого мнения, автор предложил гипотезу роста расхода из-за увеличения размеров сжатого критического сечения в струе. Этого же мнения, на наш взгляд, придерживались авторы работ [18, 19], в которых теоретически показан рост поперечного размера струи и, соответственно, коэффициента расхода вплоть до заполнения (запирания) сечения отверстия потоком. Таким образом, существуют некоторые разногласия и незначительные противоречия в объяснении роста расхода газа при его истечении из различных отверстий и сопел.
Цель статьи - уточнение, по возможности с единых позиций, характеристик установившегося движения сжимаемой жидкости в отверстиях и соплах и теоретическое исследование действительного процесса истечения из отверстий и сопел в рамках квазиодномерного приближения, т.к. возможности одномерной теории далеко не исчерпаны, а также устранение некоторых разногласий в толковании процесса истечения.
Изложение основного материала. Пусть имеется емкость перед отверстиями различной конфигурации (рис. 1, а, б), в которой находится газ, параметры которого Р0, Т0 и р0. Допускаем, что емкость имеет достаточно большие размеры, и скоростью движения газа в ней можно пренебречь. Поэтому параметры газа в емкости можно считать параметрами торможения. Давление среды в окрестности отверстия Рокр.ср.
Также как в случае истечения несжимаемой жидкости, при истечении газа вблизи выходного сечения отверстия образуется сжатое сечение Fcx. Площадь поперечного сечения отверстия - Fome.
- dmn
Р, Т, р
Рис. 1 - Типичная конфигурация отверстий (сопел): а - в толстой стенке; б - в тонкой стенке
При решении этой задачи допускаем:
1. Газ совершенный и подчиняется уравнению состояния идеальных газов:
P р = RT.
(1)
2. Процесс истечения адиабатный, т.е. выполняется условие (уравнение Лапласа-Пуассона) и соблюдается условие Т0 = const:
P
Pl
к
ро
(2)
где Р, Р0 - соответственно абсолютные статическое и полное давление газа;
р , р0 - соответственно статическая плотность и плотность «заторможенного» газа; R = 8314/и - газовая постоянная;
Т - температура газа по абсолютной шкале; р - молекулярная масса газа; к - показатель адиабаты.
3. Во всех точках сечения отверстия и сжатого сечения струи параметры газа по сечению одинаковы (параметры газа по сечению распределены равномерно).
215
ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2015р. Серія: Технічні науки Вип. 30. Т. 1
ISSN 2225-6733
В процессе истечения газа происходит потеря некоторой доли энергии из-за несовершенства поджатия (непараллельности линий тока и неравномерности распределения параметров в поперечном сечении), внутреннего трения, образования ударных волн и др. Учитывая все возможные потери энергии, в т. ч. волновые при образовании течения Дж. Тейлора - симметричного типа околозвукового течения со сверхзвуковыми включениями, примыкающими к стенкам отверстия или сопла, и (или) течения Майера с замыканием этих сверхзвуковых зон скачком уплотнения, тепловую диаграмму течения (процесс расширения) газа из сужающихся отверстий (0 < р< 180°) можно в общем случае представить в виде, приведенном на рис. 2.
Принято обычно считать, что эти потери пропорциональны кинетической энергии. В связи с этим уравнение энергетического баланса (уравнение Бернулли) для двух сечений: первое сечение в емкости, второе - сжатое сечение сразу за отверстием (или выходное сечение отверстия) (рис. 1), можно записать в следующем виде
Рис. 2 - Процесс расширения газа в соплах на /^-диаграмме
w_ +f dp+с w_ = 0
2 J р 2
(3)
или
1
w =
Р+С \
Po dP /г dP
Р
(4)
где
21Р = І 21
P Г II P
w - скорость, которую приобретает газ в струе;
Z - коэффициент сопротивления отверстия (учитывающий в т.ч. газодинамические потери);
Р и р - давление и плотность истекающего газа в сжатом (или в выходном) сечении;
І = , - коэффициент скорости.
Vі+с
Заменив в уравнении (4) р его значением из уравнения (2) и интегрируя, получим
w = І
21РРР VP J ' = І
2k
(
k -1
Pl
Po
P
P
= І
2k P0
к -1 P0
1 -
P
i-i \ k
V P0 J
(5)
Обычно параметры Р и р не даны по условию, и поэтому их нужно предварительно определять по известному давлению окружающей среды в окрестности отверстия. Допустим, что в окрестности выходного сечения давление Р = Рокр.ср, тогда подставляя это значения в (5), находим:
w = І
( 1 J
2k P0 Рокр.ср P0
k -1 Р0 V р0 pk окр.ср J
= І
2k P
k -1 P0
1 -
( P
окр.ср
k-1
'N k
P0
0 J
(6)
Но так как скорость нужно определить, то нельзя решить вопрос: будет ли скорость в струе меньше или больше скорости звука в ней? Поэтому желательно располагать каким-либо признаком, который позволил бы решить этот вопрос. Для нахождения такого признака используют то обстоятельство, что по мере снижения давления в окружающей среде скорость истечения газа возрастает, а скорость звука при этом уменьшается. Очевидно, что совпадения этих скоростей соответствует критической скорости, которую легко определить через парамет-
p
216
ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2015р. Серія: Технічні науки Вип. 30. Т. 1
ISSN 2225-6733
ры торможения:
w =
кр
2k P0
о
k +1 p0 V k +
2k
1RT0 =у/ЩКР
P
kPp.
Ркр
Приравнивая скорости (6) и (7), найдем отношение давлений
P
окр.ср
~Р0~
(7)
= п (или распола-
P0
гаемый перепад давлений
Р,
= П 0), при котором достигается критическая скорость исте-
окр.ср.
чения из отверстия (сопла) и наблюдается запертый режим, т.е. найдем
пкр =
1 -
1 к —1
2/ і • (8)
V ф k + 1)
Другими словами, поток должен быть «разогнан» с целью компенсировать потерянную скорость до такой скорости А, соответствующей отношению пкр при полном расширении потока, чтобы с учетом потерь (коэффициента скорости ф) получить звуковую скорость в потоке. Действительно, если по критическому перепаду Пкр определить скорость Aj при полном расширении газового потока до давления окружающей среды (точка Аі на рис. 2), то получим
А =
к 1 ^ к +1
1 — П к I--------
кр ) к — 1
(
1 к — 1
ф
2 к +1
к +1
к — 1
ф
(9)
Используя уравнение неразрывности, записанного для двух сечений потока (пренебрегая подмешиванием газа из окружающей среды на рассматриваемом участке течения) сжатого и горла отверстия соответственно через теоретически ожидаемые параметры и действительные, получим
Рож-^Аж) = р^ТКт.дАіеІГ) = ,, РоГР^тДятЗ)
4RTo 4RT00 VRTo
Здесь при сравнении параметров теоретического и действительного течений сделаны допущения:
1. Теоретическое и действительное сопла имеют одинаковые площади горла Fome = Fmin.
2. Полные давления P01 и температура торможения Т0 и другие теплофизические свойства в потоках газов на входе действительного и идеального сопел одинаковы.
3. В рассматриваемых сечениях сравниваемых сопл распределение параметров равномерно и соответствует равновесному состоянию рабочего тела.
А/), урав-
(10)
деист '0отв
нение можно представить в виде
Тогда рТОтР = РдотСт . Разделив все части уравнения (10) на величину Fol
М = °с
удеист І ютв
є чАсж) = дур...-
чАтттввр) ІАтп ’
(11)
где
Р0с
рдеист
р0отв
- коэффициент восстановления полного давления в сжатом сечении
потока;
є =
- коэффициент сжатия струи;
q(A) - газодинамическая функция приведенного расхода.
Или, принимая во внимание, что теоретически ожидаемая при сверхкритических перепадах давления приведенная скорость в горле сопла будет АттОр = 1,
М = °сж є ■ д(Асж ) = д(АдтГ ) . (12)
Т. е. по известному значению коэффициента расхода можно определить достаточно про-
к
~ сж
217
ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2015р. Серія: Технічні науки Вип. 30. Т. 1
ISSN 2225-6733
сто продольную составляющую средней скорости в горле отверстия (минимального сечения сопла, отверстия), справедливой при углах наклона образующей меньше 700, и пренебрегая потерями на трение до этого сечения
Я кТ )= И- (13)
Это соотношение совпадает с результатами [24] при течении газа в соплах Лаваля. При этом коэффициент скорости в отверстии фотв = Аотв. При значении п = пкр отверстие (или сопло) «запирается», что означает достижение равенства площадей сжатого и минимального сечений (F^ = FOT^ и є = 1), и коэффициент расхода приближается к максимальному значению ц** для данной геометрии сопла (рис. 3). С понижением п параметры в выходном сечении, в т.ч. и коэффициент расхода, не изменяются, т.е. происходит стабилизация линии перехода.
Уравнение неразрывности, записанное для действительного критического сечения, и теоретические параметры в этом же сечении, пренебрегая подмешиванием газа из окружающей среды на рассматриваемом участке течения, получим
pZopFAk
АТ ц AT
где к - теоретическое значение приведенной скорости в действительном критическом сечении, определяется из соотношения [24]
или ц = С
КР
• я(к)=я кг),
(14)
( Г
к =
ц2 + k2 -1
ц
(15)
у (k -1)2 k -1
Сравнивая выражения (14) и (12), получим
ц = ссж-є-я(лсж) = скр • 4Э (16)
откуда следует, что при є ^ 1 и я(ксж) ^ 1 коэффициенты ссж ^ скр и ц ^ ц**, в сжатом
сечении скорость потока близка к скорости звука. Если п = пкр, то для умеренных значений коэффициента расхода ц**
є = Ь ссж = ц**= Я(ка), (17)
где ц** - коэффициент расхода при запертом режиме истечения.
0,9 U
0,90
0,85
DM
Следует отметить, что каждому значению /і** соответствует определенное значение 7гкр, а структура потока газа в узкой части отверстия зависит, кроме влияния толщины пограничного слоя (числа Re) и степени поджатая, от формы отверстия и в первую очередь от двух факторов: угла конусности [8 сужающейся части (рис. 3 [6, 26], рис. 4 [22, 24]) и относительной кривизны стенок r/Rmin в области узкого сечения, где г - радиус скруглення стенки в узкой части отверстия (сопла), Rmin - радиус минимального сечения (горла) отверстия или сопла (рис. 5 [26]).
На рис. 6 представлены результаты расчета второго критического отношения давления 7гкр к коэффициенту расхода на этих режимах /л**. Здесь же представлены данные авторов [8, 9, 26], которые согласуются с полученными зависимостями. Следует отметить, что при (8 > 65-700 [23] до горла отверстия могут происходить потери полного давления, в т. ч. в случае внезапного сужения канала. Поэтому при малых значениях коэффициента расхода расчетные значения несколько разнятся от экспериментальных данных. Разнятся они и по причине достаточно приближенного принимаемого значения пкр в связи с асимптотической зависимостью ц** от пкр.
Если отношение Рокр ср/Р0 будет больше найденного значения пкр, то выходное сечение
«незапертое» (F^ < F), скорость истечения на выходе из отверстия будет меньше скорости звука в ней и можно пользоваться формулой (6) для ее определения. При этом если Пкр < п < п*кр,
Xft —і і
< і --£j 3
1 / □ 1 """І ) Г
□ ^ н < И —j—і
г
!j Ї *-/3 = 15° □ -J3=15° 0 —fi=U0Q расчет
е в
о о
о о _н_
О 1,0 и,0 5,0 По
Рис. 3 - Изменение коэффициента расхода с изменением полного перепада давления
218
ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2015р. Серія: Технічні науки Вип. 30. Т. 1
ISSN 2225-6733
п
*
кр
k
2 Л k-1
k + 1
(18)
то можно принять, что критическое течение (Асж = Мсж = 1) реализуется в сжатом сечении. Здесь п*кр - первое (теоретическое) критическое отношение давлений определяется по (18) или (8) при ср= 1.
Рис. 4 - Зависимость коэффициента расхода от угла наклона стенок отверстия [24]
о —І
9
) /V 2,0
^ / 4 у=зо° о ■^45° о ~ 75-90°, . Ґ/Ктіп 1,6
п ч 1,2
Рис. 5 - Влияние ра-
о 0,4 0,8 диуса скругления сте-
нок на коэффициент расхода
Рисунок 6 - Зависимость коэффициента расхода при запертом режиме от критического отношения давления для отверстий с различными углами наклона образующей: Д - водяной пар (k =1,3); ▲ - воздух (k =1,4) [8, 9]; • - воздух [26];-,
----- - теория, соответственно для (k =1,3) и (k=1,4)
Коэффициент сжатия струи в этом диапазоне, следуя уравнению неразрывности (10) для потока,
Ксж _ Р0отв „(і ЧІ^отв) _ М
тп Ч \жотв) = = .
Є =
К,,
Р
(19)
отв -* 0сж сж сж
При п = п*кр коэффициент восстановления полного давления в сжатом сечении и коэффициент сжатия
= Р0сж /Р0а = U Є = М = ), (20)
где м - коэффициент расхода при незапертых режимах.
Если п > п*кр, то, пренебрегая потерями на трение, согласно (10), при ст сж = Р0сж /Р0а = 1
219
ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2015р. Серія: Технічні науки Вип. 30. Т. 1
ISSN 2225-6733
е = V ч(4ж) = и І ч(4ж)=чіктв )/ q (h), (21)
где hj - значение приведенной скорости при полном расширении потока для данного отношения давлений ж.
Коэффициент расхода /и при незапертом режиме истечения можно определить по известному коэффициенту и** из соотношения [8]
ж - жкр + (1 - жК
и = -
и*
q теор
-ч($)=
и** q теор
¥) *V
k+1 Г, fi''
к -1
V
£=-
"кр_
1 - ж
кр .
(22)
кр
Чтеор
* + 1 ) к-1
2
ж
V
к +1
(
к -1
к-1 Л
1 -ж *
При совпадении первого и второго критических отношений давления (ж*кр = жкр), что характерно для сопел Лаваля, функция £ = ж, а при ж=жкр £ = ж*кр.
Теоретический приведенный расход чтеор = 1 при ж< ж*кр.
: > :
Предположив, что зависимость q(n) при п > пкр описывается уравнением эллипса
1
q =
1 -
f ^
ж жкр
v1 -жкр j
1 -ж
1 - 2жкр (1 -ж)-ж .
кр
Тогда
и = -
и*
1
1 - 2жкр (1 -ж)-ж^
Чтеор 1 - жкр
Таким образом, решение задачи можно представить в такой последовательности:
1. По справочным данным определяем и**.
2. По формуле (8) и (18) находим значение жкр (или Ркр) и п*кр.
(23)
(24)
3. Находим ж = -
P
окр.ср
Рп
и при пкр < п < п*кр по формуле (22) или (24) ц. При п < п,
кр
и = и**.
4. Далее определяем по известному коэффициенту расхода и (или и**) коэффициент скорости (в сжатом или выходном сечении отверстия) р.
5. Если п > пкр, то скорость находим по (6).
6. Если п < пкр и истечение происходит через отверстие в тонкой стенке (в т. ч. через суживающееся коническое сопло - отверстие с различным углом наклона стенки) или сопло имеет цилиндрический (либо плавно суживающийся) профиль, то скорость истечения, соответственно, за выходным сечением отверстия в струе и на срезе сопла будет равна критической, и ее определяют по формуле (7) для критической скорости. При этом все остальные параметры в этих сечениях также будут критическими.
Следует отметить, что в выходном сечении отверстия в тонкой стенке продольная составляющая средней скорости никогда не достигает критической скорости. Таким образом, при использовании отверстия в тонкой стенке или суживающегося конического сопла получить скорость равной и больше критической на срезе невозможно. Невозможно получить скорость больше критической и при использовании цилиндрических и плавно суживающихся профилированных сопел (без учета сдува пограничного слоя на выходе).
Для получения скорости истечения больше критической (т.е. сверхзвуковой) необходимо использовать сопла Лаваля (рис. 7), имеющие сужающуюся и расширяющуюся части. В общем случае, минимальное сечение сопла Лаваля не соответствует критическому, в котором достигается критическая скорость. Критическое сечение располагается ниже по потоку в расширяющейся части сопла (рис. 8). Для того чтобы внутри сопла параметры потока изменялись непрерывно и в потоке за пределами сопла отсутствовали скачки уплотнения (ударные волны), сопло должно иметь специальный профиль, рассчитанный на данные условия работы (давление окружающей среды, давление и температура торможения, расход газа, необходимая скорость на выходе). В таких соплах коэффициенты истечения (скорости, расхода) близки к единице (и
220
ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2015р. Серія: Технічні науки Вип. 30. Т. 1
ISSN 2225-6733
(или ц**) = р = 1) и пкр = п*кр, критическое сечение практически совпадает с минимальным. Однако на практике чаще всего сопло имеет профиль, соответствующий конической поверхности (рис. 7). При этом угол наклона образующей (угол конусности) обычно короткой сужающейся части к оси сопла /3 может изменяться в широких пределах (50 < /3 < 900), а угол раскрытия 2у конической поверхности расширяющейся части обычно не превышает 10^14° (реже до 30-400). Коэффициенты истечения (ц и ц**, р и с) зависят от угла наклона образующей суживающейся части сопла Д При этом, также как в случае профилированного сопла, пкр = п*кр (в связи с образованием за минимальным сечением зоны отрыва с давлением ниже атмосферного). Некоторое скругление конической образующей сужающейся части в области минимального сечения сопла (рис. 5) или использование цилиндрической поверхности в качестве образующей горла сопла приводит к существенному ослаблению влияния угла наклона образующей суживающейся части на коэффициенты истечения, приближая их значения к соответствующим величинам для профилированного сопла [15-17, 26].
Рис. 7 - Коническое сопло Лаваля Рис. 8 - Процесс расширения га-
за в соплах Лаваля на iS-диаграмме
Необходимо отметить, следуя [24], что коэффициент скорости р может возрастать при наличии ощутимых газодинамических потерь по сравнению с потерями на трение с ускорением потока вдоль сопла Лаваля, что согласуется с мнением [12, 17]. К примеру, при сверхкритических перепадах давления в минимальном сечении сопла Лаваля (как и в выходном сечении суживающегося сопла), пренебрегая потерями полного давления до этого сечения из-за вязких потерь,
Л
pmin = ~Л _ Лтіп и q(2min) _ р, (25)
Лкр
где р - коэффициент расхода сопла Лаваля; Хтт - приведенная скорость в этом сечении (рис. 8).
В критическом сечении, которое находится в расширяющейся части невязкого потока (или сопла Лаваля), в соответствии с данными [24]
221
ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2015р. Серія: Технічні науки Вип. 30. Т. 1
ISSN 2225-6733
л
Фкр =-
_кр_
л
л
и +
k2 -1
4-1
и
(k -1)
2
k-1
Фкр > Фп
(26)
/
Потери полного давления на участке между минимальным и критическим сечениями или, пренебрегая к тому же потерями до выходного сечения, между минимальным и выходным
P
а =^ = и-Ч(л), (27)
P
Omin
где q(A) - приведенный расход в сечении расширяющейся части сопла с теоретической скоростью Л, где достигается критическая скорость и расположено действительное критическое сечение (рис. 8); Р0кр, P0min - полное давление соответственно в критическом и минимальном сечениях сопла.
Тогда в выходном сечении ( k -1 -
Фа =
и
а
1 — ___ Л2
1 k + 1 Ла теор
1 — k - 1 т2Л2
А т ! та^а теор
k + 1
k-1
qW
1 - k-1 Л
1 k + 1 Ла теор
1 - k - 1 ф2 л2 А у ! та^а теор
k + 1
Л*-Г
Фа > Ф
кр
(28)
При известном значении коэффициента расхода сопла р (обычно определяемого экспериментально) расчет коэффициентов скорости ф, восстановления полного давления аа = Р0а/Р0 не вызывает затруднений.
Для определения действительных коэффициентов истечения из сопел в некоторых работах используется выражение, которое справедливо в сугубо теоретическом случае адиабатного расширения газа в соплах до выходного сечения, где давление в потоке должно быть равно давлению окружающей среды («квазирасчетное» истечение) [4 и др.]:
ґ k -1 2 ^
1 k +11 - Ла k +1 а
Pco) 1 4 k -1 о 2 ЛаИ
I V /
k +1 1
1 - —1 Л2 1 k+11 k+1Ла
k -1 Л2 2 k -1
п
2
ЛаИ
(29)
т.е. другими словами, определяются такие потери скорости, чтобы при той же неизменной площади выходного сечения получить давление на выходе, равное давлению в окружающей среде, что некорректно.
Представленную методику можно использовать вместе с полученной из эксперимента величиной второго критического отношения давления при запертом режиме для сопел с различным углом конусности образующих, для определения давления в зоне отрыва потока (рис. 7) за минимальным сечением в соплах Лаваля.
2
+
Выводы
1. Результаты аналитического исследования режимов истечения газа из отверстий (сопел) показали, что одномерная модель может с успехом использоваться для расчета параметров истечения из отверстий (сопел) разной конфигурации его сужающейся части; позволяет правильно оценить изменения основных параметров потока в каналах отверстий.
2. Анализ проведенных исследований и расчетов указывает на то, что коэффициенты расхода на запертых режимах истечения для сужающихся и сужающихся-расширяющихся отверстий (сопел Лаваля) одинаковы, одинаковы при этом и параметры в критическом сечении, наблюдаемом за горлом этих отверстий. При этом существуют потери полного давления, и коэффициент восстановления полного давления примерно равен коэффициенту расхода. При незапертых режимах в диапазоне от второго до первого критических отношений давления пкр < п < п*кр потери полного давления снижаются.
3. Задачами дальнейших исследований являются: уточнение ряда расчетных зависимостей для определения коэффициента сжатия струи и коэффициента восстановления полного давления при незапертых режимах истечения, учет вязких эффектов при определении коэффициентов истечения.
222
ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2015р. Серія: Технічні науки Вип. 30. Т. 1
ISSN 2225-6733
4. Полученные результаты можно использовать в качестве приложений в энергетике, металлургии и др. при разработке различных дутьевых и горелочных устройств.
Список использованных источников:
1. Меркулова Н.М. История механики газа / Н.М. Меркулова. - М. : Наука, 1978. - 232 с.
2. Франкль Ф.И. О задачах С.А. Чаплыгина для смешанных до- и сверхзвуковых течений / Ф.И. Франкль // Изв. АН СССР. Сер. математическая. - 1945. - Т.9, № 2. - С. 121-142.
3. Франкль Ф.И. Истечение сверхзвуковой струи из сосуда с плоскими стенками / Ф.И. Франкль // Докл. АН СССР. - 1947. - Т.58, № 3. - С. 381-384.
4. Газовые турбины авиационных двигателей / Г.С. Жирицкий, В.И. Локай, М.К. Максутова, В.А. Стрункин. - М. : Оборонгиз, 1963. - 608 с.
5. Benson R.S. Compressible flow through a two-dimensional slit / R.S. Benson, D.E. Pool // Int. J. Mech. Sci. - 1965. - V.7, № 5. - P. 315-353.
6. Wehover S. Transonic flow in conical convergent and convergent-divergent nozzles with nonuniform inlet conditions / S. Wehover, W.C. Moger // AIAA Paper. - 1970. - № 70. - Р. 635.
7. Соколов В.Д. Коэффициент расхода осесимметричных сужающихся сопл с произвольным контуром / В.Д. Соколов, С.В. Ягодин // Уч. зап. ЦАГИ. - 1975. - Т.6, № 1. - C. 117-121.
8. Дейч М.Е. Техническая газодинамика / М.Е. Дейч. - М. : Энергия, 1974. - 592 с.
9. О критических режимах истечения перегретого и влажного пара из сопл, щелей и отверстий / М.Е. Дейч, В.К. Шанин, В.И. Соломко, В.А. Дорошенко / Теплоэнергетика. - 1973. -№ 4. - С. 83-85.
10. Fenain M. Calcul des perfor mances d'une tuyere propulsive convergente. Comparaison avec l'ex-perience / M. Fenain, L. Dutouquet, J.-L. Solignac // Recherche Aerospatiale. - 1974. - № 5. -P. 261-276.
11. Танг С.П. Экспериментальное определение коэффициента расхода осесимметричного сверхзвукового сопла / С.П. Танг, Дж.Б. Фенн // РТК. - 1978. - Т.16. - № 1. - С. 53-60.
12. Тагиров Р.К. Теоретическое исследование течения идеального газа в сужающихся соплах / Р.К. Тагиров // Изв. АН СССР. МЖГ. - 1978. - № 2. - С. 198-202.
13. Истечение газов с различными показателями адиабаты / Г.А. Филиппов, Г.А. Салтанов, В.А. Сивобород, Ю.С. Косолапов // Изв. АН СССР. Энергетика и транспорт. - 1982. -№ 1. - С. 121-126.
14. Косолапов Ю.С. Расчетно-теоретическое исследование истечения газа из плоских щелей и осесимметричных отверстий / Ю.С. Косолапов, В.А. Сивобород. // Изв. АН СССР. МЖГ. -1984. - № 2. - С. 109-115.
15. Абалаков Г.В. Экспериментальное исследование критических расходомерных сопел с прямолинейными образующими проточной части / Г.В. Абалаков, В.М. Чефанов, А.П. Герасимов // Изв. ВУЗов. Авиационная техника. - 1987. - № 2. - С. 3-7.
16. Абалаков Г.В. Влияние диаметра входа и радиуса скругления входной кромки горловины на коэффициент расхода сопел с прямолинейными образующими / Г.В. Абалаков, В.М. Чефанов, А.П. Герасимов // Изв. ВУЗов. Авиационная техника. - 1989. - №3. - С. 78-80.
17. Тагиров Р.К. Влияние пограничного слоя на расход и удельный импульс сужающегося сопла / Р.К. Тагиров // Изв. ВУЗов. Авиационная техника. - 1988. - № 1. - С. 77-81.
18. Седов Л.И. Механика сплошной среды. Т.2 / Л.И. Седов. - М. : Наука. - 1988. - 424 с.
19. Черный Г.Г. Газовая динамика / Г.Г. Черный. - М. : Наука. - 1970. - 568 с.
20. Погорелов В.И. Газодинамические расчеты пневматических приводов / В.И. Погорелов. -Л. : Машиностроение, 1971. - 184 с.
21. Пирумов У.Г. Обратная задача теории сопла / У.Г. Пирумов. - М. : Машиностроение, 1988. - 240 с.
22. Ленцов И.А. Определение коэффициента расхода сужающихся и сверхзвуковых сопел с прямолинейными образующими проточной части / И.А. Ленцов, Н.В. Гудкова // Изв. ВУЗов. Черная металлургия. - № 1. - 1991. - С. 88-89.
23. Лухтура Ф.И. О потерях энергии при течении газа в соплах. Часть I / Ф.И. Лухтура // Вісник Приазовського державного технічного університету. - 2004. - Вип. 14. - С. 287-292.
24. Лухтура Ф.И. О потерях энергии при течении газа в соплах. Часть 2 / Ф.И. Лухтура // Вісник Приазовського державного технічного університету. - 2005. - Вип. 15. - С. 167-172.
223
ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2015р. Серія: Технічні науки Вип. 30. Т. 1
ISSN 2225-6733
25. Лухтура Ф.И. О потерях энергии при течении газа в соплах. Часть 3 / Ф.И. Лухтура // Вісник Приазовського державного технічного університету. - 2007. - Вип. 17. - С. 187-185.
26. Абрамович Г.Н. Прикладная газовая динамика. Ч. 1 / Г.Н. Абрамович. - М. : Наука, 1991. -600 с.
27. Годунов С.К. Разностная схема для двумерных нестационарных задач газовой динамики и расчет обтекания с отошедшей ударной волной / С.К. Годунов, А.В. Забродин, Г.П. Прокопов // Ж. вычисл. матем. и матем. физ. - 1961. - Т.1. - № 6. - С. 1020-1050.
Bibliography:
1. Мегкиїоуа N.M. History of a mechanics gas / N.M. Мегкиїоуа. - М. : Nauka, 1978. - 232 р. (Rus.)
2. Frankl F.I. About S.A. Tshapligin problems for blended up to and supersonic flows / F.I. Frankl // Izv. AN SSSR. Ser. matematitsheskaya. - 1945. - Т.9, № 2. - Р. 121-142. (Rus.)
3. Frankl F.I. The expiration of a supersonic jet from a vessel with flat walls / F.I. Frankl // Dоkl. AN SSSR. - 1947. - Т.58, № 3. - Р. 381-384. (Rus.)
4. Gas turbines of aero-engines / G.S. Zhiritskii, V.I. Lokay, М.К. Маksutovа, V.A. Strunkin. - М. : Oborongiz, 1963. - 608 р. (Rus.)
5. Benson R.S. Compressible flow through a two-dimensional slit / R.S. Benson, D.E. Pool // Int. J. Mech. Sci. - 1965. - V.7, № 5. - P. 315-353.
6. Wehover S. Transonic flow in conical convergent and convergent-divergent nozzles with nonuniform inlet conditions / S. Wehover, W.C. Moger // AIAA Paper. - 1970. - № 70. - Р. 635.
7. Sokolov V.D. Flow coefficient of axisymmetrical convergented nozzles with an arbitrary contour / V.D. Sokolov, S.V. Yagodin // Utsh. zap. TsAGI. - 1975. - Т.6, № 1. - Р. 117-121. (Rus.)
8. Deitsh M.E. Technical gas dynamics / М.Е. Deitsh. - М. : Energia, 1974. - 592 р. (Rus.)
9. On critical modes of the expiration an overheated and wet pair from nozzles, slots and holes / М.Е. Deitsh, V.K. Shаnin, V.I. Sоlоmkо, V.A. Dоrоshеnkо / Теplоenеrgеtikа. - 1973. - № 4. -Р. 83-85. (Rus.)
10. Fenain M. Calcul des perfor mances d'une tuyere propulsive convergente. Comparaison avec l'ex-perience / M. Fenain, L. Dutouquet, J.-L. Solignac // Recherche Aerospatiale. - 1974. - № 5. -P. 261-276. (Fr.)
11. Та^ S.P. Experimental definition of a flow coefficient of an axisymmetrical supersonic nozzle /
S.P. Та^, Jh.B. Fenn // RTK. - 1978. - Т.16. - № 1. - Р. 53-60. (Rus.)
12. Таgirov R.K. Theoretical exploration of flow of ideal gas in convergent nozzles / R.K. Таgirov // Izv. AN SSSR. MZhG. - 1978. - № 2. - Р. 198-202. (Rus.)
13. The expiration of gases with various parameters of an adiabat / G.A. Fillipov, G.A. Saltanov, V.A. Sivoborod, Yu.S. Коsоlаpоv // Izv. AN SSSR. Energetika i trasport. - 1982. - № 1. - Р. 121126. (Rus.)
14. Коsоlаpоv Yu.S. Computational theoretical exploration of gas efflux from flat slots and axisymmetrical holes / Yu.S. Коsоlаpоv, V.A. Sivoborod // Izv. AN SSSR. MZhG. - 1984. - № 2. -Р. 109-115. (Rus.)
15. Abаlаkоv G.V. Experimental exploration of critical Venturi-type nozzles with straight-line generator a flowing part / G.V. Abаlаkоv, V.M. Tshefanov, A.P. Gеrаsimоv // Izv. VUZov. Aviat-sionnaya tekhnika. - 1987. - № 2. - Р. 3-7. (Rus.)
16. Abаlаkоv G.V. Influence of an inlet diameter and radius of rounding of a lip of a port on a flow coefficient of nozzles with straight-line generator / G.V. Abаlаkоv, V.M. Tshefanov, A.P. Gеrаsimоv // Izv. VUZov. Aviatsionnaya tekhnika. - 1989. - №3. - Р. 78-80. (Rus.)
17. Таgirov R.K. Boundary layer effect on the consumption and specific impulse of a convergented nozzle / R.K. Таgirov // Izv. VUZov. Aviatsionnaya tekhnika. - 1988. - № 1. - Р. 77-81. (Rus.)
18. Sedov L.I. Mechanics of continuum. Т.2 / L.I. Sedov. - M. : Nauka. - 1988. - 424 р. (Rus.)
19. Tshernii G.G. Gas dynamics / G.G. Tshernii. - M. : Nauka. - 1988. - 424 р. (Rus.)
20. Pogorelov V.I. Gas dynamics calculations of pneumatic traction mechanisms / V.I. Pogorelov. -
L. : Mаshinostroeniye, 1971. - 184 р. (Rus.)
21. Pirumov U.G. Revertive problem of the theory of a nozzle / U.G. Pirumov. - M. : Mаshinostroe-niye, 1988. - 240 р. (Rus.)
22. Lеntsоv I.A. Definition of a flow coefficient of convergented and supersonic nozzles with straight-
224
ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2015р. Серія: Технічні науки Вип. 30. Т. 1
ISSN 2225-6733
line generator a flowing part / I.A. Lentsov, N.V. Gudkova // Izvestiya VUZov. Tshyomaya Metal-lurgia. - 1991. - № 1. - Р. 88-89. (Rus.)
23. Lukhtura F.I. About energy losses at flow of gas in nozzles. Part 1 / F.I. Lukhtura // Reporter of the Priazovskyi state technical university. - 2004. - Issue 14. - Р. 287-292. (Rus.)
24. Lukhtura F.I. About energy losses at flow of gas in nozzles. Part 2 / F.I. Lukhtura // Reporter of the Priazovskyi state technical university. - 2005. - Issue 15. - Р. 167-172. (Rus.)
25. Lukhtura F.I. About energy losses at flow of gas in nozzles. Part 3 / F.I. Lukhtura // Reporter of the Priazovskyi state technical university. - 2007. - Issue 17. - Р. 187-185. (Rus.)
26. Abramovitsh G.N. Applied gas dynamics. Part 1 / G.N. Abramovitsh. - М. : Nauka, 1991. -600 p. (Rus.)
27. Godunov S.K. The incremental scheme for the two-dimensional non-stationary problems of gas dynamics and calculation of flowing around with a stepping back shock wave / S.K. Godunov, A.V. Zabrodin, G.P. Prokopov // J. Vitsheslitelnaya matematika i matematitsheskaya fizika. -1961. - Т.1. - № 6. - Р. 1020-1050. (Rus.)
Рецензент: В.А. Масгов
д-р техн. наук, пpoф., ГВУЗ «ПГТУ»
Статья вступила 27.05.2015
225
