книги / Термодинамика
..pdfУравнению (iâ.2). Йолезно отметить, что уравнения (iâ.2) и (16.22) получены из двух различных исходных зависимо стей (механической и термической форм сохранения энер гии).
Из предыдущего видно, что понятия «энтальпия тормо жения» и «температура торможения» родственны, однако области применения уравнений (13.2) и (16.22) различны. Действительно, уравнение (13.2) не ограничено какими-ли бо предположениями о физических свойствах среды, и, следовательно, утверждение, что энтальпия торможения сохраняет свое значение вдоль всего не совершающего внешней работы адиабатного потока, всегда справедливо.
Здесь, однако, необходимо напомнить, что при термоди намическом анализе течение предполагается одномерным и рассматриваются осредненные по сечениям параметры. Со ответственно этому контрольные сечения следует проводить так, чтобы не оказались затронутыми участки, характери зуемые заметными поперечными скоростями или неодно родностями в распределении продольных скоростей. При нарушении этих ограничений могут возникнуть кажущиеся противоречия между требованиями теории и эксперимен тальными фактами. Так, в явном противоречии с получен ными уравнениями сохранения обнаруживается следую щий эффект. Если в поток не слишком малой скорости по местить тело плохо обтекаемой формы (как-то сферу, па раллелепипед, поперечно обтекаемый цилиндр и т. п.), то непосредственно за его кормой образуются зоны, в кото рых температура торможения может быть как ниже, так и выше температуры торможения набегающего потока. При этом конвективные переносы энтальпии торможения, вы численные в контрольных сечениях, проведенных перед те лом и сразу за его кормой, не равны друг другу. Между тем в указанных опытах течение может рассматриваться как стационарное, адиабатное и без производства внешней работы, так что, казалось бы, во всей области течения дол жно соблюдаться условие 2=const.
В § .13.3 разъяснено, что при термодинамическом ана лизе мы всегда полагаем d/Tp= —dqTP. При этом мы опи раемся на опыт, который во многих случаях (даже при не равномерных скоростных полях) подтверждает справедли вость условия i-\-w2/2—const (для идеального газа 7 + -f-ay2/2cp=const). Этот факт означает, что в названных случаях работа касательных напряжений d/Tp практически равна выделению теплоты трения dqTр, так что отдельные малые объемы среды могут рассматриваться как бы энер гетически изолированными от соседних, и для них полная
Энергия ôkà3biBàèïcri эквивалентной энтаЛьпий торМо&ё- ния. Однако это не всегда так: например при значительных поперечных градиентах и высоких скоростях течения исход ные предположения нарушаются настолько сильно, что ис пользование уравнений (13.2) и (16.22) приводит к суще ственным ошибкам.
Изложенное отнюдь не подвергает сомнению обосно ванность термодинамических соотношений. В рассмотрен ном примере, чтобы обеспечить правомерность применения термодинамического анализа, достаточно отодвинуть кон трольное сечение вниз по потоку за корму тела в зону, где течение выравнивается.
Вернувшись к обсуждению свойств параметров тормо жения, укажем, что с их помощью можно исключить из рассмотрения начальную скорость потока w<>. Для этого построим в i, s-диаграмме (рис. 16.4) точку 0' с координа тами iio+^2o/2 и So, соответствующую параметрам адиабат
ного и без потерь торможения в ну левом сечении канала, и все даль нейшие построения будем вести от
нее. |
Из |
сопоставления |
формул |
(16.10) |
и (16.11) видно, что подоб |
||
ный прием |
существенно |
упрощает |
|
анализ. |
|
|
|
Сделанное выше замечание о по |
|||
терях |
касается давления |
торможе |
|
ния, которое в отличие от энтальпии |
|||
торможения |
существенно |
зависит |
|
от наличия |
трения. Действительно, |
обратившись к рис. 16.4, видим, что, если торможение сопровождается потерями, конечная точ ка процесса смещается от точки 0Гвправо вдоль линии io~\- 4~tiy2o/2=const. В пределе, когда потери составят 100%, ко нечная точка процесса торможения 0" оказывается на изо баре р0. В реальных условиях восстановление давления без сколько-нибудь заметных потерь достигается против лобо вой точки помещенного в поток тела, если скорости потока дозвуковые (при сверхзвуковых скоростях набегающего потока торможение сопровождается потерями в ударной волне; более подробно об этом будет сказано ниже).
Согласно формуле (16.21) критическая скорость, т. е. скорость течения, совпадающая с местной скоростью звука, является характерной величиной для процесса течения и может быть поэтому принята в качестве масштаба. Отно шение местной скорости течения к критической Л=ш/а* будем называть приведенной скоростью.
Из формул (16.19) и (16.20) с учетом (16.16) получаем
kpv |
1 k + \ |
|
(k 1)2 £ |
|_ | pQÜ0— 2 Ai — 1» |
|
откуда следует |
|
|
P»V, |
k + 1 |
(16.23) |
|
Если рассматриваемый процесс (напомним, что здесь, как и выше, речь идет об адиабатном процессе) происхо дит без потерь, то его параметры могут быть выражены через приведенную скорость X с помощью справедливых для обратимой адиабаты формул ро/р=(р/ро)~'/к\ р/Ро= =(v/v0)~h и v/vo—pofp, подстановка которых в уравнение (16.23) дает
|
|
k - \ \knk-p |
(16.24) |
|
Ч |
‘ - k + \ K ) |
|||
|
||||
|
|
k - \ \mk-D |
(16.25) |
|
ft |
\ |
k+ 1Л ) |
||
|
||||
|
|
Строго говоря, формула (16.24), так же как и другие формулы аналогичного вида, справедлива только для иде ального газа при £=const. Однако с учетом комментариев к формулам (16.10) и (16.20) эти соотношения с хорошей степенью приближения могут быть использованы в более общих предположениях при k=£const. Несколько иначе формулируются ограничения для формулы, содержащей температуру:
^ 1j I |
(16.26) |
|
k+ 1А‘ |
||
|
Эта формула также допускает возможность непостоянства показателя k, но требует обязательного соблюдения усло вия pv=RT. Обратим внимание, что в отличие от (16.24) и (16.25) формула (16.26) справедлива для необратимой адиабаты, т. е. для процессов с трением. Действительно, формула (16.26) может быть получена непосредственно из (13.2) при единственном предположении о том, что газ является идеальным. В этом случае в формуле (13.2) про изводим замену i=kRTJ(à—1) и учитываем, что а2*= =2kRTol(k-\-\). Дальнейшие преобразования с целью пе рехода к формуле (16.26) очевидны,
Графики |
зависимостей |
(16.24) — |
(16.26) при |
£=1,4 показаны |
на рис. |
16.5. Ввиду важного значения этих со отношений они приводятся в таблич ном виде для различных значений по казателя адиабаты k в руководствах по газодинамике. Широко используют
ся также газодинамические |
функции |
||
(ГДФ), где аргументом является отно |
|||
шение pipa. |
В э т о м случае |
из |
(16.24) |
имеем |
|
|
|
* = У Щ |
\ Г^ (Р /гУ ‘- ' )1к |
(16.27) |
|
Рис. 16.5 |
|
|
|
Положив в (16.27) р = 0, убеждаемся, что при истечении в абсолютный вакуум приведенная скорость достигает ма ксимума
Я*кс = 1/ ■ ( * + ! ) / ( * - 1),
согласующегося, естественно, со значением, которое можно найти по отношению wMaKCfa., где а>ма1(с определяется из формулы (16.12).
Наряду с приведенной скоростью К в качестве безраз мерной характеристики потока часто используется отноше ние скорости потока к местной скорости звука, называемое числом М. Очевидно, М/Л,=а*/а, причем, как нетрудно убедиться, из этой пропорции следует зависимость
M = i y Ç b / / ' - Е Т Т 1'-
Применение различных масштабов для скорости тече ния имеет своим следствием различия не только внешнего вида получаемых соотношений, но и их смысла. При этом числа Х=>ху/а* представляют собой, по сути дела, те же скорости wy только в безразмерном виде, так что зависи мость w=:f(p/po) на рис. 16.2 простым изменением масшта ба вертикальной оси превращается в зависимость Х= =/(р/Ро). В отношении числа М=w/a так просто посту пить нельзя из-за изменяемости знаменателя. Особенно сильно это различие проявляется в предельных условиях при р/ро— *0: здесь величины К и w ограничены уровнями
Ямакс и Домакс, в то время как число М неограниченно воз растает. Под этим углом зрения преимущества, создавае мые применением постоянного масштаба, очевидны. Вме-
284
cfe с тем число М Дает количественную меру того, на сколько близка скорость потока к местной скорости звука.
Выше было отмечено, что малые возмущения распро страняются в среде независимо от ее собственного движе ния (т. е. от основного движения потока). В связи с этим Для неподвижного наблюдателя скорость убегающей по
потоку |
звуковой |
водны равна w-\-a, |
скорость |
встречной |
||
волны есть w—а. |
Очевидно, |
когда |
w |
достигнет |
а (т. е. |
|
М =1), |
звуковые |
волны не |
смогут |
распространяться про |
тив течения, так что через источник звука пройдет перпен
и>=0
Рис. 16.6
дикулярный течению фронт, на котором энергии последо вательных во времени звуковых волн будут складываться. Этот фронт, мощность которого может заметно превышать мощность отдельного звукового импульса, называется ударной волной. При сверхзвуковых скоростях потока (М>1) ударная волна сносится вниз по течению и, как показано на рис. 16.6, образует клин, угол раскрытия кото рого тем меньше, чем больше число M(sin a=ax/wz=a/w= = 1/М), где т — время.
Вернувшись к формулам (16.23) и (16.26), отметим, что при Х=1 значения приведенных параметров определяются исключительно значением показателя адиабаты. Эти зна
чения называются критическими и для течения без потерь определяются по формулам
( р у ) * ____2__ р * ____ / 2 \VO-D |
р*___ / 2 |
y/tt-D |
ч |
А »о f t + 1 ’ />о V* + 1 / |
Ро V* + |
V |
I |
___2_ |
|
|
I |
Г,—fe+1» |
|
|
i |
|
(16.28) |
причем |
последняя формула справедлива только при p v = |
— RT, |
но в отличие от первых трех она сохраняет силу |
также и для потока с потерями. Для того чтобы и осталь ные формулы приведенных параметров использовать для расчета процессов с потерями, обратимся к рассмотренно
му выше принципу суперпозиции. Так, из фрмулы |
(16.24) |
|||
найдем (напомним, что ц — КПД процесса) |
|
|||
k— 1 X2 |
\*/ (*—1) |
(16.29) |
||
k + \ |
Г, |
) |
||
|
||||
Для отношения удельных объемов из (16.23) следует |
||||
о/о. = (1 - |
/ |
Р(Я. ч). |
(16.29а) |
При течении с потерями критическую скорость попрежнему определяем условием Х=1. Соответственно для критического отношения давлений из (16.29) получаем
М'-^тГ^Ч-ттр'’ <16-30>
т. е. вследствие потерь критическое отношение давлений падает. Для критического отношения плотностей имеем
р* |
* + 1 / , |
k — \\kHk-\) |
( 2 |
\4(k-\) |
Qп |
t = |
T - ( ' - i + i ) |
< ( г р |
) |
С6-30») |
Что же касается первой и последней формул (16.28), то они сохраняют свой вид и для течения с потерями.
16.3. Исследование процесса истечения
До сих пор мы не касались вопроса о связи между из менением состояния потока и профилем направляющего канала. В случае течения несжимаемой жидкости (напри мер, воды, нефти и т. п.) задача решается просто. Из урав нения постоянства расхода
G=p/7o/=const |
(16.31) |
Получаем, Принимая бо внимание условие p=const, зави симость
Fw=zconst,
Из которой ясно, что Для увеличения скорости течения не обходимо площадь проходного сечения канала уменьшить, Для уменьшения скорости — увеличить.
Гораздо сложнее дело обстоит при течении сжимаемой
среды, поскольку в этих условиях в (16.31) |
приходится |
Следить за изменениями трех величин: F, р, |
w. Проводя |
Такое исследование, рассматриваем общий случай стацио нарного течения сжимаемого газа с трением, теплообменом Н производством технической работы. При этом остановим ся только на одной стороне вопроса — основных условиях Перехода через критическую скорость. В нашем исследова нии не будем выходить за рамки идеально-газового при ближения, тем более что распространение анализа на ре альные газы не меняет качественных выводов.
При анализе будем исходить из уравнений сплошности (16.31), сохранения энергии (13.5) и сохранения механи ческой энергии (Бернулли) (13.7), к которым присоединим Уравнение состояния p v = R Î и формулу di=kRdT/(k—1). Все уравнения будем рассматривать в дифференциальной форме. Таким образом, имеем
dF/F-\-dw/w-\- ф /р = 0 ;
dq= [k/(k— \)\ RdT-\-d(w2/2) -j-d/T; dp/p+d(w2/2) +Л,+<М,р= 0; dp/p=dp/p+dT/T.
Сопоставляя первую, третью и последнюю формулы, находим
RdT+d(w2/2) -Н/т-М /тР= а 2 (dF/F+dw/w) /к.
Далее с помощью второй формулы исключим dT. Вводя число М, после преобразований окончательно получаем
dF/F—dlv/a2— (к—1) dq/a2~ k d U /a 2= (М2—1 ) dw/w.
(16.32)
Каждое из слагаемых левой части (16.32) отражает не который специфический способ воздействия на скорость потока, стоящую в правой части уравнения. Если все чле ны в левой части, кроме первого, равны нулю и, следова тельно, единственным способом воздействия на скорость остается изменение поперечного сечения канала, то получим формулу
dF /F = (№ ~l)dw iw . |
(16.33) |
Здесь возможны два случая в зависимости от того, ускоряется вниз по потоку движение (dw>0) или замедля ется (dw< 0).
Допустим сначала, что движение ускоряется (dw>0). Положим, что начальная скорость мала (М <1). Тогда со гласно (16.33) получаем, что dF<0, т. е. вниз по течению канал должен сужаться, что согласуется с представления ми гидравлики. Однако если М>1, то из (16.33) следует, что d F > 0, т. е. оказывается, что при сверхзвуковых ско ростях для ускорения потока вниз по течению канал дол жен расширяться. Этот, на первый взгляд, парадоксальный
результат |
следует |
понимать |
как следствие |
того факта, |
|
что в уравнении |
(16.31) |
относительное падение плотности |
|||
начинает |
превалировать |
над |
относительным |
увеличением |
скорости.
При переходе от dF< 0 к dF> 0 сечение канала, оче видно, имеет минимум. Из уравнения (16.33) видно, что минимальная площадь (dF=0) соответствует М =1. Таким образом, при течении без потерь в наиболее узком месте канала (горле) устанавливаются критические параметры. До критического (минимального) сечения канал сужается
(конфузорный участок), |
за |
ним — расширяется |
(диффу- |
зорный участок). Очевидно, |
оба участка — конфузорный и |
||
диффузорный — в месте |
их |
сопряжения должны |
иметь |
общую касательную, так что в окрестностях горла должно
удовлетворяться условие {dF/dx)x^ Xi, = 0, |
где х — коорди |
|
ната, направленная вдоль оси канала. |
(dw<0). В этом |
|
Рассмотрим замедляющееся движение |
||
случае при малых начальных скоростях потока |
(М <1) ка |
|
нал должен быть расходящимся {dF>%0), |
как, |
например, |
для воды. Если же начальная скорость является сверхзву ковой, то для замедления потока необходимо применить
сходящийся профиль |
(конфузор) и, только после того, как |
||
скорость упадет до |
звуковой, перейти на |
расходящийся |
|
канал (диффузор). |
|
{dw> 0 и dw < 0), можно утверж |
|
Обобщая оба случая |
|||
дать, что при переходе |
через критические |
условия для |
обеспечения монотонности процесса знак изменения площа ди канала необходимо изменить на противоположный. При этом при дозвуковых скоростях потока профилирование трубы качественно отвечает принципам гидравлики, при сверхзвуковых скоростях — противоположно им.
В заключение отметим, что короткую конфузорную тру бу с плавным входом принято называть сходящимся или простым соплом. Комбинацию простого сопла с диффузо ром называют соплом Лаваля в честь шведского ннжене-
288
ра, впервые применившего такую конструкцию Для Полупи ния сверхзвуковых скоростей при подаче пара на турбин ное колесо.
Вернемся снова к уравнению (16.32) и допустим те перь, что поперечное сечение канала сохраняется неизмен ным (dF=0), внешний теплообмен и трение отсутствуют {dq^dUр=0). Таким образом, в левой части остается только одно слагаемое, представляющее внешнюю работу, и уравнение (16.32) приобретает вид
—dlT/a2= (М2—1 )dw/w. |
(16.34) |
Отсюда следует, что при М<1 для возрастания скорости (dw>0) необходимо производство потоком работы против внешних сил (d/T> 0 ). Таким образом, если дозвуковой по ток производит внешнюю работу, например вращая турби ну, встроенную в рассматриваемую цилиндрическую трубу, то скорость потока растет. Одновременно падает его тем пература, что видно из записанного для этих условий урав нения сохранения энергии
к/ (k— 1 )RdT= —wdw—d/T, |
(16.35) |
в котором, как только что было указано, |
d w > 0 и dl?>0. |
Поскольку температура падает, падает также и скорость звука. Уменьшение скорости звука и одновременный рост скорости потока — оба эти фактора способствуют увеличе нию числа М, которое в конце концов достигнет значения М =1. С этого момента знак разности М2—1 изменяется на обратный, и если мы по-прежнему стремимся к ускоре нию потока (d w > 0), то, очевидно, знак dlT необходимо изменить на отрицательный. Отрицательная техническая работа (г//т< 0) способствует, например, ситуация, когда встроенный в трубу пропеллер приводится во вращение от постороннего двигателя, создавая некоторую тягу в пото ке. Из (16.35) получим wdw= —di—d/T, отсюда видно, что при сверхзвуковых скоростях увеличение кинетической энергии не может обеспечиваться только за счет уменьше ния энтальпии, а требует еще дополнительного подвода ме ханической энергии извне.
Таким |
образом, |
в цилиндрической трубе |
монотонное |
||
увеличение |
скорости |
с переходом от дозвукового течения |
|||
к сверхзвуковому может быть осуществлено |
с помощью |
||||
механического |
воздействия, причем знак |
воздействия дол |
|||
жен меняться |
при М =1. В качественном |
отношении опи |
санная картина аналогична тому, что происходит в сопле Лаваля, и если сопло Лаваля называют геометрическим
19-3038 |
289 |
Ёоплом, tô цилиндрическую трубу, в которой с itoMoüiUô знакопеременной технической работы осуществляется пе реход через М =1, можно назвать механическим соплом.
Примем теперь, что dF=dlT=dlvp=0. Это значит, что единственным воздействием на поток газа, движущийся ро цилиндрической трубе постоянного сечения, является внешний теплообмен. В этих условиях из (16.32) получим
— [(А—1 )/а2] dq— (М2—1 )dw/w, |
(16.36) |
откуда видно, что увеличение скорости при дозвуковом ре жиме (М2—1<0) требует подвода теплоты к потоку извне (dq>0). Чтобы судить, как при этом изменится число М, необходимо сопоставить относительные изменения скорости потока (dw/w) и скорости звука (da/a). По определению М=w/a, так что
dM./M.=dw/w—da/a=dw/w—0,5 (dT/ T), (16.3?)
где при переходе ко второму равенству использована фор мула a2=kRT. Выразим dw/w через df/T . Для этого обра тимся к уравнению состояния dT/T=dp/p — dp/p, в кото ром в данном случае из уравнений сплошности и Бернулли имеем соответственно dp/p= —dw/w и
dp |
dw |
j ______ dw________ __ dw |
|
|
p Œ |
w |
p |
kpv W |
W • |
Используя эти значения, получаем окончательную фор
мулу |
(16.38) |
d T /T = ( l —kW)dw/w, |
из которой следует, что при малых числах М ускорение по тока, происходящее под влиянием подвода теплоты, сопро
вождается |
ростом температуры. |
При М > 1 /|/й |
знак dT |
||
меняется, |
и дальнейший |
разгон |
потока |
сопровождается |
|
падением |
температуры. |
Таким |
образом, |
при |
М >1 / У k |
увеличение скорости потока сопровождается уменьшением скорости звука, и, следовательно, число М здесь безуслов но растет. Чтобы получить возможность судить о поведе
нии числа М при его значениях, меньших чем 1 /]Пг, под ставим (16.38) в (16.37). Получим
цм |
1 + *м* |
dr |
(16.39) |
|
~ Ж ~ 2(1 — Ш*) |
т> |
|||
|
||||
отсюда ясно, что при М < 1 /|/гЛ |
и dT> 0 число М также |
|||
растет. |
|
|
|
290