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

Круглов А.Б.Харитонов В.С.Радовский И.С. Руководство по технической термодинамике с примерами и задачами

.pdf
Скачиваний:
98
Добавлен:
21.07.2024
Размер:
3 Мб
Скачать
☆

при низком давлении пар можно считать идеальным газом и пре-

небречь удельным объемом

жидкости v' << v" = RT/ρ. Тогда

dp/p = rdT/(RT 2). При r = const

после интегрирования получим

lnp = −r/(RT) + const.

 

П.4.3. Поскольку фазовый переход осуществляется при постоянном давлении, то теплота фазового перехода r равна разности энтальпий пара и воды r = h" – h'. На линии насыщения имеем

ddTr = dTh′′ − dTh′ .

Поскольку dh = Tds + vdp и Tds = cdT, то вдоль линии насыщения (dr/dT) = c" – c' + (v" – v')(dp/dT). Используя соотношение Максвелла, для теплоемкости на линии насыщения можно получить:

 

 

ds

 

∂s

 

 

∂s dp

 

 

 

 

 

∂v

 

dp

 

c = T

 

 

 

= T

 

 

+ T

 

 

 

 

 

= cp − T

 

 

 

 

 

 

.

dT

 

 

 

 

 

 

 

dT

Тогда

 

∂T p

 

 

∂p T dT

 

 

 

 

 

∂T p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr

 

 

 

 

 

 

 

 

 

∂(v′′ −v′)

 

 

dp

 

 

 

 

 

 

= c′′ − c′ + v′′ − v′ −T

 

 

 

 

 

 

 

 

=

 

 

 

dT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂T p dT

 

 

 

 

 

 

 

 

= c′′ − c′ + 1

−

 

∂ln(v′′ −v′)

 

 

r

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂ln T

 

p T

 

 

 

 

 

 

 

 

В частности,

при

низком

давлении

v' << v" = RT/ρ,

и тогда

(dr / dT) =c′p′ −c′p . Теплоемкость жидкости обычно выше, чем теп-

лоемкость пара, поэтому теплота испарения убывает с ростом температуры.

П.4.4. Используя уравнение Клайперона−Клаузиуса, получим

(dp/dT)пл = −13,4 МПа/K.

П.4.5. Работа кругового процесса равна площади цикла в координатах p, v или T, s. Для элементарного цикла из условия задачи имеем (v′′ −v′)dp = (s′′ − s′)dT . Отсюда дифференциальное уравнение для линии насыщения:

141

dp

=

s′′ − s′

=

 

r(Tн)

 

,

 

 

 

 

 

 

 

 

v′′ −v′

T (v′′(T ) −v′(T ))

dT

н

 

 

 

 

 

 

 

н

н

н

 

где индексы «′» «′′» обозначают жидкость и пар соответственно на линии насыщения.

П.4.6. а) Уравнение спинодали в параметрической форме (с обобщенным объемом в качестве параметра):

π =

3ω − 2

; τ = −

(3ω −1)

2

ω2

4ω3

.

 

 

 

б) Согласно второму началу термодинамики работа изотермического кругового процесса abodea равна нулю (см. рис. к задаче) и, следовательно,

π(ω2 − ω1) = ω∫2 3ω8τ−1 − ω32 dω

ω1

или

π(ω

 

− ω ) = 8τ ln 3ω2

−1

+ 3

 

1

−

1

.

(П.4.1)

2

 

 

 

 

 

 

1

3 3ω1

−1

 

 

 

 

 

 

 

 

 

 

 

ω2

 

ω1

 

С другой стороны, для точек a и e (жидкость и пар на линии насыщения) имеем:

π =

8τ

 

−

3

=

8τ

 

−

3

.

(П.4.2)

3ω −1

ω2

3ω

 

−1

 

 

 

 

2

 

ω2

 

 

1

 

 

1

 

 

 

 

 

2

 

 

Из уравнений (П.4.1) и (П.4.2), исключая π и τ, получим трансцендентное отношение между объемами жидкости ω1 и насыщенного пара ω2:

 

3ω

 

−1

 

ω

 

− ω

 

 

1

 

 

 

1

ln

 

2

 

=

 

2

1

2

+

 

 

 

+

 

 

.

 

 

 

3ω1

−1

3ω2

 

 

3ω1 −1

 

ω2 + ω1

 

 

 

−1

Далее определим давление и температуру на линии насыщения:

π =

3ω1ω2 − ω1 − ω2

,

τ =

(ω1 + ω2)(3 ω1 −1)(3ω2 −1) .

 

 

(ω ω

)2

 

 

8(ω ω

)2

 

1

2

 

 

 

1

2

 

П.4.7. Для кругового процесса в окрестности тройной точки (рис. П.4.1) имеем:

142

v∫dQ = rпл + rисп − rсуб = 0 .

Поскольку в соответствии со вторым началом термодинамики работа изотермического цикла равна нулю, отсюда для теплоты сублимации имеем:

rсуб = rисп + rпл .

ж

rпл rис

т г

rсуб

Рис. П.4.1

П.4.8. Подставив численные значения соответствующих параметров в уравнение Клайперона–Клаузиуса, можно получить при Тпл = 273 К (dp/dT)пл = −13,4 МПа/К. Таким образом, с увеличением давления температура плавления льда понижается.

П.4.9. Уравнение связи перегрева поверхности T = Tc − Tн от-

носительно температуры насыщения с критическим радиусом парового зародыша а имеет вид

T = 2σTНv′′ . ar

При низких давлениях пар можно считать идеальным газом,

 

 

 

RT

 

 

 

2σRT 2

 

v′′ =

 

н

 

и

T =

 

н

.

 

 

 

 

 

 

 

 

p

 

 

 

arp

Можно заметить, что перегрев

Т возрастает по мере уменьше-

ния давления.

При

рн = 0,1 МПа,

 

Тн = 373 К, v′′ = 1,72 м3/кг

(R = 462 Дж/(кг К)),

r = 2,26 106 Дж/кг, σ = 60 10–3 Дж/м и

Т = 5° получим

а = 6,8 мкм.

Если

рн = 0,007 МПа, Тн = 313 К,

σ = 69 10–3 Дж/м и а = 6,8 мкм, то

Т = 65°.

П.4.10. При длительном кипении жидкости происходит «истощение» крупных центров, поверхность обедняется центрами парообразования. В действие включаются так называемые «потенциально активные центры» с меньшими характерными размерами, вследствие чего перегрев поверхности увеличивается.

П.4.11. Основное уравнение термодинамики в приложении к магнетику TdS =dU − HdB преобразуем к виду для потенциала F

(свободная энергия) в координатах Т, В: dF =−SdT + HdB .

143

Отсюда, используя свойство потенциала H =∂F /∂B , для критического поля получим

Hкр = FnB−крFs ,

поскольку в сверхпроводнике В = 0.

С другой стороны, S =−∂F /∂T и для разности энтропий в n и s

состояниях

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dHкр

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂(F − F )

 

 

 

 

 

 

 

 

 

 

 

 

 

S

n

− S

p

= −

 

n

 

s

 

= − B

 

.

 

 

 

 

 

 

 

 

 

∂T

 

dT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

кр

 

При H

кр

= H

0

1−T 2

/ T 2

 

и

B

 

= μ

0

H

кр

получаем

 

 

 

 

 

 

(

 

 

 

 

 

кр)

 

кр

 

 

 

 

 

 

S

n

− S

p

= μ

 

H

2

 

2T

1−T

2

/ T

2

,

откуда теплота фазового пере-

0

0

 

 

 

 

 

 

 

 

 

 

 

 

T 2 (

 

 

 

кр)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

кр

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

хода

2

rsn =T(Sn − Sp) = μ0H02 2TT2 (1−T 2 / Tкр2 ).

кр

П.4.12. Поскольку объем жидкости остается постоянным при растяжении поверхности, то дифференциал свободной энергии равен dF = −SdT + σdΩ. Отсюда, используя свойства потенциала, получим

∂∂ΩS = − ddTσ ,

и теплота изотермического увеличения поверхности lΩ = T ∂∂ΩS = −T ddTσ > 0 ,

поскольку поверхностное натяжение σ убывает с ростом температуры.

5. Термодинамика потока

П.5.1. а) Умножая уравнение сохранения импульса на скорость, получим уравнение сохранения механической энергии:

144

 

d w2

 

 

w

 

dp

 

 

 

 

 

 

 

= −

ρ

 

dz

− l

− lтр ,

(П.5.1)

 

dt 2

 

 

 

 

 

 

где l = fw – работа против внешних сил;

lтр = fтрw – работа про-

тив сил трения (в единицу времени и на единицу массы).

В уравнении (П.5.2) перейдем от дифференциала к полной производной:

dh

= q′ −

1

dp

,

(П.5.2)

dt

 

ρ

dt

 

 

где q′ = q + qтр – полная теплота, подводимая к 1 кг газа в единицу

времени от внешнего источника и вследствие трения. Исключая из полученных уравнений слагаемые, содержащие давление, получим

d h + w2 = q − l . dt 2

б) Без совершения полезной работы в адиабатическом потоке сохраняется сумма удельной энтальпии и кинетической энергии:

h + w2 / 2 = h0 = const ,

где h0 – энтальпия «торможения», т. е. энтальпия при давлении и температуре заторможенного потока (p0; T0; w = 0).

В случае идеального газа h0 − h = cp(T0 −T) , для адиабатиче-

ского процесса ρ = ρ0Θn (где Θ = T / T0 , n = cV / R) . Массовый расход газа определяется соотношением M = ρwΩ (где Ω – площадь поперечного сечения канала). Скорость потока

w = 2cpT0(1−Θ) .

Отсюда

 

M

 

 

 

 

 

 

 

 

 

 

 

 

= ρ

0

2c T ΩΘn

(1−Θ) = const .

 

 

 

 

 

 

p

0

 

 

 

 

 

Для минимального критического сечения

 

 

Θ

кр

= 2n / (1+2n) = 2 / (1+k); p / p = Θcp /R

= Θk/(k−1) .

 

 

 

 

 

 

 

 

кр

0

кр

кр

в) Местная скорость течения в критическом сечении

 

 

w

 

=

 

2c T (1−Θ

кр

) =

kRT

/ M .

 

 

 

кр

 

 

p 0

 

 

кр

 

 

145

П.5.2. T = 250 K. Термометр покажет температуру, близкую к начальной температуре.

П.5.3. В системе координат, связанной с самолетом, торможение набегающего потока воздуха описывается уравнением энергии

h + w2 /2 = h0 . Считая воздух идеальным газом, для которого

h0 − h = cp(T0 −T) , получим

 

 

 

 

 

T

w2

 

 

k −1

M 2,

0 = 1 +

 

= 1

+

 

 

2c T

2

 

T

 

 

 

 

 

p

 

 

 

 

 

где M = w/a – число Маха;

a =

kRT –

адиабатическая скорость

звука.

Для воздуха k = 1,40. При M = 1 увеличение относительной температуры при торможении равно

T0/T = 1,2.

П.5.4. pн = 4,0 МПа; T = 250 °С.

П.5.5. а) Из основного уравнения термодинамики следует, что (∂T / ∂p)h = −(λp + v) / cp . С помощью соотношений Максвелла можно показать, что λp = −T(∂v / ∂T) p . Отсюда для коэффициента Джоуля−Томсона получается выражение:

j =

1

T

 

∂v

 

− v

 

=

v

(α T −1) .

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

cp

 

∂T p

 

 

 

cp

б) Для идеального газа αp = 1/T и, следовательно, j = 0.

в) Используя для нахождения производной (∂v/∂T)p уравнение Ван-дер-Ваальса, получим необходимое выражение для j.

г) p j = 2a / (bv j) − 3a / v2j .

д) Используя условие экстремума (∂pj/∂vj)0 = 0, получим v0j = 3b = vкр; p0j = a/(3b2) = 9pкр; T0j = 8a/(9bR) = 3Tкр.

П.5.6.Если кинетическая энергия струи воды, истекающей из щели, мала по сравнению с ее энтальпией, то процесс адиабатиче-

146

ского дросселирования (без совершения полезной работы) можно считать изоэнтальпическим, h = const:

u1 + p1v1 = u2 + p2v2 .

Если v1 v2, то u2 − u1 = (p1 − p2)v1 > 0 , и внутренняя энергия воды при дросселировании увеличивается. Повышение температуры воды составляет величину

 

(p

− p

)v

(20 −0,1)10610

−3

T − T =

1

2

1

=

 

 

= 4, 7 K .

 

c

 

4, 2 103

 

2 1

 

 

 

 

 

 

 

V

 

 

 

 

 

6.Термодинамический анализ процессов преобразования энергии

П.6.1. Указание. Представить произвольный цикл в виде совокупности элементарных циклов Карно; воспользоваться теоремой Карно, представив доказываемое соотношение в виде

Q = qVV = KV (TЯ − Т1) .

П.6.2. Максимальная работа совершается в обратимом процессе. Единственно возможный обратимый процесс, исключающий теплообмен газа с окружающей средой при конечной разности температур, состоит из адиабатического участка 1→2′ и изотермы 2′→2

(рис. П.6.1) при температуре окружающей

T

 

 

среды. Работа расширения в адиабатиче-

 

 

 

ском процессе равна убыли внутренней

T1

 

1

энергии L1→2′ = U1 − U2′ , а в изотермиче-

 

 

 

ском – убыли

свободной

энергии

 

 

 

L2′→2 = F2′ − F2 .

Кроме того, над окру-

2

 

2′

жающей средой производится работа из-

T2

 

 

менения объема при давлении p′, равная

s2

s1

s

p′(V2 – V1). В итоге получается выраже-

 

 

 

ние:

U + T′ S − p′

 

Рис. П.6.1

 

Lmax = −

V.

 

 

 

П.6.3. В изобарическом процессе 3–4–1 (см. рис. 6.25 к задаче) к 1 кг рабочего тела подводится теплота q1 = h1 – h3 = cp(T1 – T2) + r1

147

(r1 – теплота испарения при температуре T1). Изменение энтропии при этом равно

s1 − s3 = cp ln(T1 / T2) + r1 / T1 .

(П.6.1)

Работа в адиабатическом процессе расширения пара в турбине

1–2 равна l = h1 – h2 = cp(T1 – T2) + r1 – x2r2 (r2 – удельная теплота испарения при температуре Т2; x2 – массовое паросодержание на

выходе из турбины в точке 2. Энтропия в точке 2 равна

 

 

 

 

 

 

 

 

 

 

s2 = s3 + r2x2 / T2 .

 

 

 

(П.6.2)

 

Поскольку s1 = s2, то из уравнений (П.6.1) и (П.6.2) можно найти

x2 и l:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

= T2

c

 

ln T1 +

r1

;

l = c

T

− T

− T

ln T1

 

+ r

1 −

T2

.

2

p

 

 

 

 

r2

 

T2

 

 

 

p 1

2

2

 

1

 

 

 

 

 

 

 

T1

 

 

 

 

T2

 

 

 

T1

 

Термодинамический КПД ηt цикла Ренкина получим, подставив

в ηt

=

l

 

значения l и q1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П.6.4. а) Полный дифференциал энергии Гиббса открытой системы, совершающей внешнюю электрическую работу εdQe, равен:

dФ = − SdT + Vdp − εdQe .

При постоянном давлении и постоянной температуре отсюда следует, что эдс топливного элемента определяется выражением

ε = − dФ , dQe

где dQe – элементарный электрический заряд, прошедший через внешнюю цепь. В водородно-кислородном топливном элементе в

расчете на 1 моль окислителя O2

перетекает заряд, равный заряду

4-х молей электронного газа Qe = 4NAe, т.е.

ε = −

Ф

.

 

 

4N e

 

A

б)

 

∂ε

=

S

 

;

 

∂ε

= −

V

,

 

 

 

 

 

 

 

 

 

 

 

4N

 

e

 

4N e

 

 

∂T

p

A

 

 

∂p

T

 

 

 

 

 

 

 

 

 

 

 

A

 

148

где V и S – изменение объема и энтропии веществ, участвующих в химической реакции. Если вода удаляется в парообразной форме, то

 

V = 2V

 

 

− (2V + V ) ≈ − RT

< 0.

 

H O

 

 

H

 

O

 

p

 

 

 

2

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

в)

ηt =

L

=

Ф

=1 +

T S .

 

 

H

 

 

 

 

Qp

 

 

Qp

 

П.6.5. Для изоэнтропного обратимого процесса в открытой системе (dl)s = -(dh)s. При одинаковом перепаде давления внутренний КПД равен отношению работ политропного и изоэнтропного процессов:

η0i = ((∂h / ∂p))x . ∂h / ∂p s

Для идеального газа dh = сpdT и ∂h/∂p = сp(∂T/∂p). Используя уравнение состояния идеального газа pV = RT, уравнение адиабаты

pV cp/cV = const и политропы идеального газа pV n = const , полу-

чим:

η0i = (1−n−1)(1−cV / cp )−1 .

П.6.6. Работа, совершаемая 1 кг рабочего тела в обратимом адиабатическом процессе, равна убыли энтальпии l1,2 = h1 – h2 .

а) Для аналогичного необратимого процесса l1,2′ = h1 − h2′ ;

l1,2′ = l1,2

 

−T2

s , где s –

увеличение энтропии за счет необрати-

мых процессов.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б) Приращение сухости

 

пара

можно

найти

из

уравнения

r2 x = T2

 

s

 

 

(где r2 – теплота парообразования при температуре

насыщенного пара T2):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

s

 

(1−η0i)l1,2

 

 

 

 

l1,2′

 

T

s

 

 

 

x

 

x =

 

2

 

 

=

 

 

;

η

0i

=

 

 

=1 −

2

 

=1

− r

 

 

.

 

 

 

 

r

l

 

 

(h

− h )

 

 

 

r

 

 

 

 

 

 

 

 

l

 

 

2

 

 

 

 

2

 

 

 

2

 

 

 

 

 

1,2

 

 

1,2

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

149

 

 

 

 

 

 

 

 

 

П.6.7. Согласно первому началу ТД работа кругового процесса (на 1 кг) равна l = v∫dq = q1 − q1 . К 1 кг рабочего тела подводится

теплота q1 = h1 − h3 = c′p T + r1 ( T = T1 −T2 , r1 – теплота испарения при температуре Т1) в парогенераторе и отводится теплота q2 = h2 − h3 = r2x2 в конденсаторе. Сухость пара в конце процесса

расширения для обратимого процесса определим из условия изоэнтропичности:

c′

ln T1

+

r1

=

r2x2

, откуда получим

x

2

= T2

c′ ln T1

+

 

r1

.

 

 

 

 

p

T2

 

T1

 

T2

 

 

 

 

 

 

 

 

r2

 

p

T2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T1

Окончательно для КПД цикла Ренкина получим выражение

 

 

 

 

 

 

T

 

c′pT2

 

T

 

 

 

c′p T

−1

 

 

 

 

 

 

 

 

 

 

η =

 

2

+

 

ln

1

1+

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

r1

 

 

 

 

 

r1

 

 

 

 

 

 

 

 

 

 

 

 

 

T1

 

 

T2

 

 

 

 

 

 

 

 

 

 

Работа

 

 

 

расширения

 

 

 

двухфазной

 

 

 

 

смеси

dl = −vdp = −[v′′x + v′(1− x)]dp ,

а

для

пара

dl′′ = −v′′xdp . По-

скольку v′′ >>v′, то влага в турбине практически работы не производит.

П.6.8. Для изобарического процесса в проточной системе первое начало имеет вид dQ = Gdh , где G – массовый расход теплоносителя; h – удельная энтальпия. В однофазной области течения

T2′

T2

 

 

T ′

T1′

T

Tн

 

T1

 

 

Q dQ

Рис. П.6.2

dh = cpdT, dQ = GcpdT ; в двухфазной области h = h′′x + h′(1− x) ( (h′′, h′ – энтальпии воды на линии

насыщения, x – массовое расходное паросодержание) и dh = (h′′− h′)dx =

= rdx, dQ = Grdx, T = Tн (r – тепло-

та испарения, Тн – температура насыщения). Процесс передачи теплоты в парогенераторе может быть представлен диаграммой в коорди-

натах (Т, Q) (рис. П.6.2).

150