Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Термодинамика

..pdf
Скачиваний:
12
Добавлен:
12.11.2023
Размер:
13.57 Mб
Скачать

Уравнению (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, видим, что, если торможение сопровождается потерями, конечная точ­ ка процесса смещается от точки вправо вдоль линии 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+ )

 

 

 

Строго говоря, формула (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/wda/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

Соседние файлы в папке книги