Добавил:
ac3402546@gmail.com Направление обучения: транспортировка нефти, газа и нефтепродуктов группа ВН (Вечерняя форма обучения) Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

теор / Термодинамика и теплопередача, Калинин, Купцов, Лопатин

.pdf
Скачиваний:
135
Добавлен:
01.06.2021
Размер:
7.08 Mб
Скачать

Термодинамика в технологических процессах…

131

Показатель изотермического процесса:

nt = 0,9585 + 6,280 10-5 T – 0,4498р + 0,1863T0,15 р. (1.379) Показатель изоэнтальпийного процесса:

nh = 0,8929 + 2,447 10-4 T1,1 + 1,480 10-2 р1,4 – 1,035 10-4 T р. (1.380)

Коэффициент сжимаемости:

 

z = 1,015 – 7,080 10-5 T – 0,4708р + 0,1936T0,15 р.

(1.381)

Плотность, кг/м3:

 

ρ = p·106/(pv).

(1.382)

ρ = 107,0 – 682,7T-0,5 – 67,22р0,1 + 4,817 104 T-1,5 р,

(1.382a)

Удельная энтальпия, кДж/кг:

 

h = –16,91 + 2,161T – 181,1р + 72,58T0,15 р.

(1.383)

Удельная энтропия, кДж/(кг.К):

 

s = 8,013 + 7,057 10-3 T – 0,7192р + 0,2519T0,15 р.

(1.384)

В соотношениях (1.375) ÷ (1.384) и далее размерность температуры Т [K], t [oC] и абсолютного давления р [МПа].

Сопоставляя точность эмпирических соотношений для расчета характеристик метана (1.375) ÷ (1.384) с точностью построения термодинамических таблиц, можно отметить, что только определение энтропии немного не отвечает требованиям.

Для наиболее распространенных составов природных газов с содержанием метана 95 100%, отличие природного газа от метана можно учесть, введя в расчеты величину мольного содержания метана в газе rмет , выраженную в долях единицы.

Обобщение экспериментальных данных по метану и двум составам природного газа позволило получить эмпирические зависимости для расчета термодинамических характеристик природного газа с содержанием метана rмет = 0,95 1.

Удельная изобарная теплоемкость, кДж/(кг.К):

cp = 3,347 – 7,561 10-4 Т1,4 + 0,4332р1,1 – 5,460 10-3 Т0,76 р1,05 +

 

+ 7,777 10-10 T3,7 r0,92 .

(1.385)

мет

 

Потенциальная функция, Дж/кг:

 

рv = –180916 + 25033T0,5 – 4,175 104 р +

 

+ 1,224 104 T0,205 р – 2391T0,65 r2,05 .

(1.386)

мет

 

132

Часть 1

Произведение удельной изобарной теплоемкости на коэффициент ДжоуляТомсона, кДж/(кг МПа):

cp Dh = 58,90 – 2,245 10-4 Т2 – rмет(0,4153Т – 1,588p – 21,67rмет +

 

+ 4,277 10-3 T p – 8,144 10-4 Т2 + 2,196 10-2 p2).

(1.387)

Показатель адиабаты:

 

 

k=1,658–6,534·10-3 Т0,79+0,3132·p1,06–0,1563Т0,14 p +1,021 10-5 T1,8 r .

(1.388)

 

мет

 

Показатель изотермического процесса:

 

 

nt = –1,464 + 0,4959Т0,3 – 0,5532р + 0,1559Т0,15 p –

 

1,536 10-8 T2,91 r0,73

+ 0,1466p1,071 r0,5 .

(1.389)

мет

мет

 

Показатель изоэнтальпийного процесса:

nh=2,007+9,446 10-4 Т1,2+1,434 10-2 p1,4–7,752 10-5 Т p1,10,1089T0,5 rмет . (1.390)

Коэффициент сжимаемости:

z = –0,2262 + 9,093T-0,33 – 0,4837p + 0,1988T0,15 p – 766312T-2,7 rмет1,67 . (1.391)

Плотность, кг/м3:

ρ = p·106/(pv),

 

 

(1.392)

ρ = –78,72 + 35597T-1,2 + 30,54p0,89 + 2,234 10-2 T1,2 p +

 

+ 3,652 10-4 T2 r

+ 0,3647p2 r

– 0,1328 Tp r0,498 .

(1.392a)

мет

мет

мет

 

Удельная энтальпия, кДж/кг:

 

 

h = –13395 + 2,667T – 34,77p + 7,817 10-2 T p – 29,90 10-4 Т2 0,5882p2 +

+ rмет(26936 – 13636rмет +3,065 10-3 T p + 21,76 10-4 Т2 + 0,5927p2).

(1.393)

Следует отметить, что полученные соотношения (1.385) ÷ (1.393) в целом точнее ныне используемой системы расчетных соотношений для определения термодинамических характеристик. Средние погрешности уменьшены в 1,3 2 раза, величин максимальных ошибок уменьшены в 1,5 3 раза.

Так как точность определения потенциальной функции pv по соотношению (1.386), лучше точности определения коэффициента сжимаемости z и плотности ρ, то при известном составе природного газа целесообразно z и ρ определять через значения потенциальной функции.

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

Термодинамика в технологических процессах…

133

Проведенный анализ показал, что разность энтальпии газа h2 в центробежном нагнетателе целесообразней определять непосредственно по соотношению (1.393).

До настоящего момента для этой цели использовалось

h1 соотношение:

h = cp1 (t2 t1) – (cp Dh)2. (р2 p1).

(1.394)

Оценка точности расчета разности энтальпии метана, по сравнению с эталоном hэ, приведена в таблице 1.7.

В качестве эталона принят расчет по термодинамическому соотношению

(1.394) с определением величин cp1

и (cp Dh)m2 по данным Госстандарта.

 

 

 

 

 

 

 

 

 

 

Таблица 1.7

Расчет разности энтальпии h, кДж/кг при сжатии метана в ЦБН

 

 

 

 

 

 

 

 

 

 

 

 

 

Параметры

 

 

Варианты сравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

3

 

4

 

5

 

6

 

7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Исходные данные:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Температура на входе, К

280

270

275

 

290

 

300

 

280

 

280

 

 

 

 

 

 

 

 

 

 

 

 

Температура на выходе, К

300

290

290

 

310

 

320

 

310

 

310

 

 

 

 

 

 

 

 

 

 

 

 

Давление на входе, МПа

4,0

4,0

3,75

 

4,5

 

4.5

 

5,0

 

5,5

 

 

 

 

 

 

 

 

 

 

 

 

Давление на выходе, МПа

5,0

5,0

4,5

 

5,6

 

5,5

 

7,0

 

7,5

 

 

 

 

 

 

 

 

 

 

 

 

Результаты расчета:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Степень сжатия

1,25

1,25

1,2

 

1,244

 

1,222

 

1,4

 

1,364

 

 

 

 

 

 

 

 

 

 

 

 

Политропный КПД

0,769

0,740

0,83

 

0,775

 

0,714

 

0,787

 

0,702

 

 

 

 

 

 

 

 

 

 

 

 

hэ, кДж/кг

40,00

39,60

29,30

 

40,36

 

42,00

 

58,80

 

60,25

 

 

 

 

 

 

 

 

 

 

 

 

h1, кДж/кг

40,10

39,38

29,17

 

40,64

 

42,00

 

59,24

 

60,62

 

 

 

 

 

 

 

 

 

 

 

 

δ, %

0,25

-0,56

-0,45

 

0,69

 

0

 

0,75

 

0,61

 

 

 

 

 

 

 

 

 

 

 

 

h2, кДж/кг

40,03

39,54

29,29

 

40,38

 

41,82

 

58,71

 

59,93

 

 

 

 

 

 

 

 

 

 

 

 

δ, %

0,08

-0,15

-0,03

 

0,05

 

-0,43

 

-0,15

 

-0,53

 

 

 

 

 

 

 

 

 

 

 

 

При хорошей точности, предлагаемые соотношения (1.385) ÷ (1.393) просты в реализации, их удобны в расчетах.

Рассмотренная система расчета термодинамических характеристик действует в указанном диапазоне давлений и температур. Следует отметить, что аналогичные соотношения могут быть получены и для других интервалов давлений и тем-

134

Часть 1

ператур природного газа. Это является актуальной задачей, так как современные магистральные газопроводы спроектированы и работают под гораздо большим давлением газа.

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

Полученные соотношения (1.385) ÷ (1.393) используются в методиках и программах расчета оценки технического состояния и определения технологических показателей работы газоперекачивающих агрегатов, внедренных в ряде газотранспортных предприятий страны, отраслевой методике нормирования топливно-энергетических затрат на газотурбинных компрессорных станциях.

Теплопередача в технологических процессах

135

 

 

Часть 2.

ТЕПЛОПЕРЕДАЧА В ТЕХНОЛОГИЧЕСКИХ ПРОЦЕССАХ НЕФТЯНОЙ И ГАЗОВОЙ ПРОМЫШЛЕННОСТИ

Теплопередача (теория теплообмена) – наука, изучающая процессы передачи теплоты между телами, распространение теплоты в пространстве и распределение температуры в телах [5, 15, 17].

Знание законов теплообмена позволяет определить количество передаваемой теплоты, интенсифицировать теплообмен в одних случаях и затруднить его в других, правильно конструировать и оптимально эксплуатировать машины, аппараты и другие технические устройства, в которых рабочие процессы сопровождаются теплообменом.

2.1. Основные положения теории теплообмена

При соприкосновении двух тел с различной температурой происходит обмен энергии движения частиц этих тел (молекул, атомов, свободных электронов) в результате, которого интенсивность движения частиц тела с меньшей температурой увеличивается, а интенсивность движения частиц тела с более высокой температурой уменьшается: тело с меньшей температурой нагревается, а другое остывает. Следовательно, для возникновения процесса теплообмена между двумя телами или между различными областями пространства, заполненного вещественной средой, необходимо наличие разности температур.

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

Различают три основные формы передачи теплоты: теплопроводность, конвективный теплообмен и лучистый теплообмен.

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

Конвективным теплообменом называется форма передачи теплоты в пространстве, осуществляемая при перемещении макрочастиц жидкости или газа. При перемещении в пространстве различно нагретых частиц жидкости или газа

136

Часть 2

 

 

происходит непосредственное их соприкосновение, поэтому здесь имеет место теплопроводность. Следовательно, конвективный теплообмен представляет собой совокупное действие двух процессов – конвекции и теплопроводности.

В зависимости от причины, вызывающей движение жидкости, различают

конвективный теплообмен при свободном движении (свободная конвекция) и конвективный теплообмен при вынужденном движении жидкости (вынужденная конвекция).

Свободная конвекция возникает вследствие разности плотностей неравномерно

нагретых слоев

жидкости

или газа в

поле

сил

тяготения. Возникновение

и интенсивность

свободной

конвекции

зависит

от

разности температур, рода

ифизических свойств жидкости, геометрии тела и объема пространства в котором протекает процесс.

Вынужденная конвекция возникает под влиянием внешнего воздействия (например, ветра, насоса, компрессора, вентилятора и т.д.), которое создает перепад давления. Интенсивность процесса зависит от рода и физических свойств среды, ее температуры, скорости движения, формы и размеров пространства, в котором происходит движение.

Тепловым излучением называется процесс переноса энергии в пространстве электромагнитными волнами.

Лучистым теплообменом называется форма передачи теплоты излучением между телами, которая включает последовательное превращение внутренней энергии одного тела в энергию излучения, распространение ее в пространстве

ипревращение энергии излучения во внутреннюю энергию другого тела.

Рассмотренные формы передачи теплоты во многих случаях осуществляются совместно, что необходимо учитывать при расчете теплообмена в технологических процессах нефтяной и газовой промышленности.

Температурное поле

Процесс теплообмена может иметь место только при условии, что в различных точках тела (или системы) температура неодинакова. В общем случае процесс распространение теплоты в телах и теплообмен между телами сопровождается изменением температуры, как в пространстве, так и во времени.

Совокупность значений температуры t в данный момент времени τ для

всех точек пространства, определяемых координатами

x, y, z, называется тем-

пературным полем

 

t = f (x, y, z,τ) .

(2.1)

Уравнение (2.1) является математическим выражением температурного поля, записанным в неявной форме. Различают стационарное (установившееся) и нестационарное (неустановившееся) температурные поля.

Стационарное температурное поле наблюдается в том случае, когда температура в различных точках пространства не изменяется во времени. Если температура изменяется во времени – температурное поле называется нестационарным.

Теплопередача в технологических процессах…

137

 

 

Температурное поле может быть функцией трех, двух и одной координаты. Соответственно оно называется трех-, двух и одномерным.

Уравнение простейшего одномерного стационарного температурного поля имеет следующий вид:

t = f (x) .

(2.1а)

В соответствии с классификацией температурного поля принципиально различают стационарный и нестационарный процессы передачи теплоты.

Аналитическое исследование теплообмена сводится к изучению простран- ственно-временного изменения температуры, т.е. нахождению уравнения (2.1) в явном виде.

Температурный градиент

Геометрическое место точек, имеющих одинаковую температуру называется изотермической поверхностью. Пересечение тела плоскостью в сечении дает семейство изотерм.

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

Рассмотрим две близко расположенные изотермические поверхности с температурами t и t + t (рисунок 2.1). При перемещении вдоль изотермической поверхности с температурой t изменение температуры не наблюдается. Перемещаясь же по направлению x в сторону изотермы, соответствующей значению температуры t + t , мы будем наблюдать изменение температуры. При этом наибольшее изменение температуры на единицу длины будет наблюдаться в направлении нормали n к изотермической поверхности. Интенсивность изменения температуры в направлении по нормали к изотермической поверхности характеризуется градиентом температур.

Рис. 2.1. К определению температурного градиента, изотермических линий и теплового потока

Температурный градиент есть вектор, направленный по нормали к изотермической поверхности в сторону возрастания температуры и численно рав-

138

Часть 2

 

 

ный пределу отношения изменения температуры t термами по нормали n (К/м)

 

 

t

 

t

grad t = lim

 

 

=

 

.

 

 

n0

 

n

 

n

к расстоянию между изо-

(2.2)

В случае трехмерного температурного поля суммарный температурный градиент определяется по правилу сложения векторов

grad t = i grad tx + j grad ty + k grad tz ,

(2.3)

где i, j, k – единичные векторы в направлении x, y, z .

Тепловой поток

Количество теплоты Qτ , проходящее в единицу времени через изотермическую поверхность F, называется тепловым потоком Q, Вт. Тепловой поток, приходящийся на единицу площади поверхности, называется удельным тепловым потоком, плотностью теплового потока q, Вт/м2.

q =

δQτ

.

(2.4)

dF dτ

 

 

 

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

Q = q dF ,

(2.5)

F

где dF – элемент изотермической поверхности, м2.

Тепловой поток Q, а также плотность теплового потока q являются векторами, за положительное направление которых принимается направление по нормали

кизотермической поверхности в сторону уменьшения температуры (рисунок 2.1).

2.2.Теплопроводность при стационарном температурном поле

Вчистом виде теплопроводность наблюдается в твердых телах, а также в тонких прослойках жидкости и газа при отсутствии в них движения. В металлах перенос теплоты осуществляется путем движения (диффузии) свободных электронов, передача теплоты за счет упругих колебаний кристаллической решетки второстепенна. В жидкостях и твердых телах – диэлектриках теплопроводность осуществляется упругими волнами.

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

Теплопередача в технологических процессах…

139

 

 

Стационарная теплопроводность

Необходимым условием распространения теплоты является неравномерность распределения температуры в рассматриваемой среде.

Согласно закону Фурье, количество теплоты проходящей через элемент изотермической поверхности dF за промежуток времени dτ , пропорционально температурному градиенту

δQ = −λ grad t dF dτ = −λ

t

dF dτ ,

(2.6)

τ

n

 

где λ – коэффициент пропорциональности, физический параметр вещества, называемый коэффициентом теплопроводности, Вт/(м·К); dF – элементарная площадь поверхности теплообмена, м2; dτ – временной промежуток, с.

Тепловой поток, согласно закону Фурье

Q = q dF = −λ

t

dF .

(2.6а)

 

 

n

 

FF

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

Q = −λ

 

t

F .

(2.7)

 

 

 

 

n

 

q = −λ

 

t

.

(2.8)

 

 

 

 

 

n

 

Знак минус в правой части уравнений (2.7) и (2.8) указывает на то, что тепловой поток и температурный градиент, как векторы, имеют противоположные направления.

Полное количество теплоты

 

τ

 

 

n

 

 

∫∫

 

 

Q = −

 

 

λ

t

dF dτ .

(2.9)

 

 

 

τ

 

 

 

 

 

 

 

0

F

 

 

 

 

Выражение плотности тепловых потоков в направлении осей x, y, z может быть записано в виде:

q

 

= −λ

t

;

q

 

= −λ

t

;

q

 

= −λ

t

.

(2.10)

x

 

y

y

z

 

 

x

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

Вектор теплового потока с учетом (2.10) для трехмерной задачи составит

q = i qx + j qy + k qz .

(2.10а)

140

Часть 2

 

 

Коэффициент теплопроводности

Под коэффициентом теплопроводности понимают тепловой поток, передаваемый через единичную площадь поверхности при единичном значении температурного градиента

λ =

δQ

 

F grad t .

(2.11)

Коэффициент теплопроводности для каждого тела имеет свое численное значение и зависит от природы, пористости, влажности, давления, температуры и других параметров. При выводе уравнения (2.7) принято, что λ не зависит от температуры и является постоянной величиной. Однако, как показывают опыты, для многих материалов с достаточной для практики степенью точности, зависимость коэффициента теплопроводности от температуры можно принять линейной, во всем рассматриваемом интервале температур

λ = λ0 (1± b t),

(2.12)

где λ0 – коэффициент теплопроводности при температуре t0 ; b – постоянная, характеризующая приращение (уменьшение) λ материала при повышении его температуры на 1 градус.

Наихудшими проводниками теплоты являются газы. Коэффициент теплопроводности газов возрастает с увеличением температуры и изменяется в пределах 0,005–0,5 Вт/(м·К).

Коэффициент теплопроводности жидкостей лежит в пределах 0,07–0,7 Вт/(м·К) и, как правило (за исключением воды и глицерина), уменьшается с увеличением температуры.

Наилучшими проводниками теплоты являются металлы, у которых λ = 10 ÷ 420 Вт/(м·К). У большей части чистых металлов с возрастанием температуры он уменьшается.

Дифференциальное уравнение теплопроводности

Распределение температуры в теле описывается дифференциальным уравнением теплопроводности, которое, при принятых допущениях, а именно – тело однородно и изотропно; физические параметры тела постоянны во времени и пространстве; температурные деформации рассматриваемого элементарного объема малы по сравнению с самим объемом; внутренние источники теплоты распределены в рассматриваемом объеме равномерно; макрочастицы тела неподвижны относительно друг друга; имеет следующий вид:

где a = λ cp ρ

t

 

q

2t

 

2t

 

2t

 

q

 

 

 

= a 2t +

v

= a

 

 

+

 

 

+

 

 

 

+

v

,

(2.13)

τ

cp ρ

x

2

y

2

z

2

cp ρ

 

 

 

 

 

 

 

 

 

 

 

– коэффициент температуропроводности, характеризующий ско-

рость изменения температуры в любой точке тела, м2/c; cp – изобарная тепло-