Термодинамические свойства воды и водяного пара

Термодинамические свойства воды и водяного пара

Задание

Часть 1. Определить допустимое количество воды, сбрасываемой летом ГРЭС в пруд-охладитель, имеющий размеры 750 х 750 м.

Скорость ветра W, м/с, температура воды и воздуха tв.

Расчеты выполняются для случаев, когда окружающий воздух имеет относительную влажность: а) 𝜑=𝜑1; б)𝜑=𝜑2.

Часть 2. В процессе эксплуатации парового котла его стальные кипятильные трубы диаметром d1/d2снаружи покрылись слоем сажи толщиной 𝛿с , а внутри слоем накипи толщиной 𝛿н .Температура дымовых газов, омывающих кипятильные трубы поперечно - tг , давление пара в котле Рн . Определить температуры на границах между слоями стенки и тепловой поток на 1 пог. м на 1м2 наружной и внутренней поверхности трубы. Представить данные расчета в виде графической зависимости.

Примечание:

накипь богата силикатами.

марка стали Ст. У12.

Оглавление

Введение

Часть 1.Исходные данные

Часть 2.Исходные данные

.1Постановка задачи теплообмена

2.2 Расчет параметров теплообмена

.3 Выполнение нулевого приближения

.4 Выполнение первого приближения

.5 Выполнение второго приближения

Заключение

Перечень ссылок

Приложения

 тепловой поток теплообмен вода

Введение

В дисциплине изучаются закономерности процессов теплообмена <https://ru.wikipedia.org/wiki/%D0%A2%D0%B5%D0%BF%D0%BB%D0%BE%D0%BE%D0%B1%D0%BC%D0%B5%D0%BD>, сопровождающихся переносом вещества, то есть, массообменом <https://ru.wikipedia.org/wiki/%D0%9C%D0%B0%D1%81%D1%81%D0%BE%D0%BE%D0%B1%D0%BC%D0%B5%D0%BD>.

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

В данной курсовой работе предоставлено решение двух проблем.

Во-первых, это расчет допустимого количества воды, сбрасываемой ГРЭС в пруд-охладитель, а, во-вторых, определение теплового потока на 1 м2 трубы.

Часть 1.Исходные данные

Определить допустимое количество воды, сбрасываемой летом ГРЭС в пруд-охладитель, имеющий размеры 750 х 750 м.

Скорость ветра W=12 м/с, температура воды и воздуха tв=21 ̊С.

Расчеты выполняются для случаев, когда окружающий воздух имеет относительную влажность: а) 𝜑1=18; б)𝜑2=65.

Решение

Эта задача относится к типу задач по процессам масообмена при вынужденной конвекции от плоской пластины. Прежде чем подобрать соответствующие безразмерные соотношения для числа Шервуда Sh (т.н. диффузионное число Нуссельта NuD), следует определить, будет ли течение воздуха на поверхности пруда ламинарным или турбулентным. Значение критерия Рейнольдса Re при достижении воздухом конца пруда равно:

значение коэффициента кинематической вязкости воздуха при 21 ̊С;

νt=0 ̊C = 13.3 *10-6 м2/с;

νt=100 ̊C = 23 *10-6 м2/с;

νt=21 ̊C = *21+13,3*= м2/c.=  =  = .

Следовательно, течение воздуха полностью турбулентное(Re˃2300).

Число Шмидта Sc (диффузионный критерий Прандтля Prd) равно:

 Sct=0 ̊C = 0.707=100 ̊C = 0.688=21 ̊C =  * 21+0,707 = 0,703.

Число Шервуда Sh равно

Sh = 0,035*Re0.8 * Sc0.33 =  = 322379.20.

Коэффициент диффузии равен

Dt=0 ̊C = 6.76 м2/c=100 ̊C = 12.1 м2/c=21 ̊C =  *21+6.76*= 7.881* м2/ч = 2,19* м2/c.

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

𝛽 =  =  = 0.0094 м/с.

Определим концентрации паров воды у поверхности озера и в окружающем воздухе. У поверхности воды воздух насыщенный и его относительная влажность равна 100 %. Из таблиц для водяного пара давление насыщения при 21  ̊C равно Рнас = 2,4877кПа. Концентрация водяного пара при 100% влажности равна:

 =  =  = 0,01832 кг/м3.

а) При относительной влажности окружающего воздуха 18% концентрация водяного пара равна

 =  =  = 0,003297 кг/м3.

Массовая скорость испарения воды равна:

m = 𝛽*S( -  ) = 0.0094*7502*(0.01832-0.003297) = 79.52кг/с.

б) При относительной влажности окружающего воздуха 65% концентрация водяного пара в воздухе равна

 =  =  = 0,01191 кг/м3.

Массовая скорость испарения воды равна

m = 𝛽*S( -  ) = 0.0094*7502*(0.01832-0.01191) = 33,94 кг/с.

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

Часть 2. Исходные данные

В процессе эксплуатации парового котла его стальные кипятильные трубы диаметром d1/d2снаружи покрылись слоем сажи толщиной 𝛿с , а внутри слоем накипи толщиной 𝛿н .Температура дымовых газов, омывающих кипятильные трубы поперечно - tг , давление пара в котле Рн . Определить температуры на границах между слоями стенки и тепловой поток на 1 пог. м на 1м2 наружной и внутренней поверхности трубы. Представить данные расчета в виде графической зависимости.

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

Постановка задачи теплообмена

Рис.1 Теплопроводность через цилиндрическую стенку.

В общем случае дифференциальное уравнение теплопроводности для бесконечного цилиндра при отсутствии внутренних источников тепла (qv = 0) имеет вид:

(1)

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

(2)

 Граничные условия имеют вид (см. рис.1):

Граница дымовые газы - сажа:

r = rc

(3)

где: 𝜆с - коэффициент теплопроводности сажи; qконв = 𝞪г(tп - tн) - плотность теплового потока (конвективная составляющая) qизл - составляющая потока излучения дымовых газов.

Граница насыщенный пар - накипь:

r = rн

(4)

где: 𝜆н - коэффициент теплопроводности накипи; qконв2 = 𝞪конв2(tп - tн) - плотность теплового потока(конвективная составляющая омывания насыщенным паром покрытой накипью стенки трубы).

Граница накипь - металл:

r = r2

(5)

где: 𝜆м - коэффициент теплопроводности металла(стальной трубы).

Граница металл - сажа:

r = r1

(6)

Нужно учесть, что рассматривается нелинейная задача, т.к. 𝜆 = 𝜆(t).

Для определения количества тепла, прошедшего через цилиндрическую поверхность (S) в единицу времени, воспользуемся законом Фурье:

С другой стороны, после интегрирования уравнения (2) получим:

Постоянные интегрирования С1 и С2 определим из граничных условий (для металлического слоя цилиндрической стенки):

при r = r1 t = tr1

при r = r2 t = tr2

Тогда имеем систему уравнений с двумя неизвестными:

Решением этой системы будет:

После подстановки С1 и С2 получим общее выражение для распределения температуры по толщине цилиндрической стенки, подчиняющееся логарифмическому закону:

Тогда с учетом выражения получаем выражение для определения плотности теплового потока цилиндрической стенки(в общем случае):

(7)

Полное количество тепла, проходящее через боковую поверхность цилиндрической стенки площадью S = 2𝜋rl (l - длина цилиндрической стенки), удобно отнести к 1 погонному метру стенки, тогда уравнение (7) необходимо умножить на 2𝜋r

(8)

Для многослойной стенки процесс теплопередачи можно рассчитать следующим способом(см. рис.2):

Рис. 2 - Распределение температур

Граница дымовые газы - сажа:

Теплопередача от греющей среды к саже определяется по закону

Ньютона-Рихмана:

конв1 = 𝛼г(tг - td3)  (9)

где 𝛼г - коэффициент теплопередачи от дымовых газов к саже;г - температура дымовых газов;- температура поверхности слоя сажи (на рис. td3 = t4).

Полное количество тепла, проходящее через d = d3 длиной 1 м определяется следующим образом:

(10)

В данном случае 𝛼г состоит из конвективной и лучистой составляющей, т.е.

𝛼г = 𝛼конв1 + 𝛼изл

Теплопроводность через сажу:

Количество тепла через слой сажи определяется на основании уравнения (8)

(11)

Теплопроводность через металл:

Количество тепла через металл трубы определяется из уравнения (8)

(12)

Теплопроводность через накипь:

Количество тепла, проходящее через накипь, определяется следующим образом:

(13)

Теплопроводности от накипи к насыщенному пару определяется законом Ньютона:

(14)

На основании закона сохранения энергии при стационарном тепловом режиме

(15)

тогда

 (16)

Решая совместно систему уравнений (16) получим:

(17)

Выражение (17) можно записать в следующем виде:

Q = K(tг - tп) (18)

где К - коэффициент теплопередачи

(19)

Выражение (18) можно записать в другом виде:

(20)

где Rп - термическое сопротивление теплопроводности от накипи к пару,

Rн - термическое сопротивление теплопроводности накипи,

RM - термическое сопротивление материала трубы,

с - термическое сопротивление теплопроводности сажи,

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

Коэффициент теплоотдачи 𝛼г состоит из конвективной и лучистой составляющей, т.е.:

(21)

Коэффициент теплоотдачи 𝛼конв1 можно определить из критериального уравнения

(22)

где Nu -критерий Нуссельта; Re - критерий Рейнольдса; Prtг - критерий Прандтля при t = tг, Prtст - критерий Прандтля при t = tст.

Коэффициенты с, m, n выбираются в зависимости от формы обтекаемого тела и от диапазонов чисел Рейнольдса. Коэффициент с зависит также от температурного фактора и учитывается сомножителем в уравнении (22).

Число Рейнольдса (характер набегающего потока при внешнем обтекании цилиндра - ламинарный, переходной или турбулентный) можно определить по формуле

(23)

где dc - определяющий размер, в данном случае dc = d3 , νг - коэффициент кинематической вязкости газа.

Критерий Прандтля (характеризует соотношение между скоростями формирования полей скоростей и температур движущегося потока) можно определить по формуле:

(24)

где 𝛼г - коэффициент температуропроводности дымовых газов.

Помимо уравнения(22), критерий Нуссельта можно определить следующим образом:

(25)

Отсюда

(26)

Лучистая составляющая 𝛼изл определяется из закона Стефана-Больцмана:

(27)

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

(28)

где 𝜀пр - приведенная степень черноты стенки 𝜎0 = 5,67* Вт/м2К4.

Подставив выражение для 𝜀пр , окончательно получим:

(29)

𝜀ст - степень черноты стенки, 𝜀г - степень черноты газа при температуре газа. Для упрощения расчетов при заданных температурах дымового газа tг и стенки можно сделать допущение, приняв

Здесь 𝜀СО2 и 𝜀Н20 - степень черноты углекислоты и водяного пара, содержащихся в объеме дымовых газов одновременно, зависят от средней длины луча РіSi,

𝜉 - поправочный коэффициент, учитывающий неодинаковую степень влияния толщины газового слоя Si и парциального давления Рі на єнергию излучения водяніх паров. Излучательная способность газов, молекулы которых состоят из одного или двух атомов, столь незначительна, что в теплотехнических расчетах они принимаются лучепрозрачными.

Для расчета коэффициента теплоотдачи 𝛼конв2 при пузырьковом кипении воды в окисленных стальных трубах можно воспользоваться формулой (2)

(30)

где qn - плотность теплового потока, Вт/м2, Рн - давление насыщенного пара, бар (1 бар = 100кПа).

Расчет параметров теплообмена.

Расчет диаметров:= 0,18 м;= 0,162 м;= d1 + 2𝛿c = 0,18+2*0,0011 = 0,1822 м;= d2 - 2𝛿н = 0,162 - 2*0,0013 = 0,1594 м.

Определение теплопроводности (из данных [1]):

𝜆с = 0,09 Вт/мК, табл.1; 𝜆н = 0,175 Вт/мК

Определяем температуры насыщения ( из данных [2])

Рнас = 1,82 МПа;нас = 207,5  ̊С (табл.4)

Для температуры дымовых газов tг:

Теплопроводность смеси дымовых газов (путем интерполяции данных [1])

𝜆N2 = 66,45* Вт/мК

𝜆H20 = 106,625* Вт/мК

𝜆C02 = 73,5* Вт/мК

 Вт/мК

Плотность смеси дымовых газов ( путем интерполяции данных [1]):

𝜌N2 = 0.3265 кг/м3

𝜌Н2О = 0.2095 кг/м3

𝜌СО2 = 0,5152 кг/м3

𝜌г = 0,3265*0,64+0,2095*0,165+0,5152*0,175) = 0,333 кг/м3

Теплоемкость смеси дымовых газов (путем интерполяции данных [1]):

СN2 = 1.166*103 Дж/кг* с

СН2О = 2,327*103 Дж/кг* с

ССО2 = 1.2427*103 Дж/кг* с

СPг = (1.166*0.64+0.165*2.327+0.175*1.2427)*103 = 1.347*103

Кинематическая вязкость смеси дымовых газов газов ( путем интерполяции данных [1]):

νN2 = 128*10-6 м2/c

νН2О = 140,75*10-6 м2/c

νСО2 = 81,82*10-6 м2/c

νг = (128*0,64+140,75*0,165+81,82*0,175)*10-6 = 1,19*10-4 м2/c

Коэффициент температуропроводности смеси дымовых газов:

𝛼г =  м2/c

Критерий Прандтля:

Prг =  =  = 0,732.

2.4 Выполнение нулевого приближения

Принимаем температуры

Расчет характеристик газов для tст = 400 ̊С [1]

Теплопроводность смеси газов ( путем интерполяции данных [1]):

= 0,0558 Вт/мК

Плотность смеси газов( путем интерполяции данных [1]):

 = 0,519 кг/м3.

Теплоемкость смеси газов ( путем интерполяции данных [1]):

 = 1,228 * 103 Дж/кгК

Кинематическая вязкость смеси газов ([1]):

 = 55,38 * 10-6 м2/с

Коэффициент температуропроводности смеси дымовых газов:

 =  = 7,97 * 10-5 м2/с

 =  = 0,694

Reг =  = 16678,91(0) = 0.035* * Reг0,8 * ( = 0,035*0,7320,33 *16678,910,8 * = 76,36

=  =  = 35,15 Вт/м2К.

Определение парциального давления трехатомных газов:

РСО2 =  = 0,175 * 105 Па

РН20 =  = 0,165 * 105 Па

Эффективная длина луча:

l = 3,6  = 3,6 = 3,6 = 1.08 м.

Средняя длина луча:

 l РСО2 = 1,08*0,175*105 = 18900 м*Н/м2РН2О = 1,08*0,165*105 = 17820 м*Н/м2

Определение степени черноты излучающих трехатомных компонентов при температуре tг и tст [3]

 tг = 775  ̊C

tст = 400  ̊C

𝜀со2 = 0,1397

 = 0,12

𝜀Н2О = 0,4307

 = 0,48

поправочный коэффициент 𝜉 = 1,1 [3]

Степень черноты смеси газов при температурах tг и tст

Определение коэффициента теплоотдачи излучением

 = 65,08 Вт/м2К.

 Вт/м2К.

 Принимая за нулевое приближение плотность теплового потока q = 10000Вт/м, определяем коэффициент теплоотдачи от накипи к пару:

𝛼П = 3,14q0.7 *  = 3,14*100000,7 *18,20,15 = 3061,54 Вт/м2К.

tм =  =  = 300  ̊C

𝜆м = 41,4 Вт/мК.

Определение термических сопротивлений:

 = 0,0174 м  ̊С/Вт

 =0,00065 м  ̊С/Вт

 = 0.0215 м  ̊С/Вт

 = 0.000405 м  ̊С/Вт

 = 0,0147 м  ̊С/Вт

∑R(0) = 0.0174+0.00063+0.0215+0.000405+0.0147=0.0547 м  ̊С/Вт

Количество тепла, проходящее через 1 метр цилиндрической поверхности.

q =  =  = 10372,393 Вт/м

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

 =  = 19979.38 Вт/м2.

=  =  = 17479,22 Вт/м2.

Расчет температур на границах между слоями стенки:

 = 775-10372,393*0,0174 = 594,11  ̊С

 = 775-10372,393*(0.0174+0.0215) = 371,17  ̊С

 = 775-10372,393*(0.0174+0.0215+0.000405)=

= 366,97  ̊С

 = 775-10372,393 * (0.0174 +0.0215 +0.000405 +0,0147)) = 214,27  ̊С

 Определение погрешности вычислений нулевого приближения:

𝜀 =  * 100% =  *100% = -3,7 ˃2%

2.5 Выполнение первого приближения

Принимаем за начальные температуры нулевое приближение:

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

q = 10372,393 Вт/м.

Определение среднеинтегральной теплопроводности металла трубы: (3)

 =  =  = 369,07   = 37,927 Вт/мК

Расчет теплофизических характеристик газов для tст = 594   [1]

Теплопроводность смеси дымовых газов (путем интерполяции, аналогично нулевого приближения):

 Вт/мК

Плотность смеси дымовых газов (путем интерполяции ):

 0,395*0,64+0,6228*0,175+0,2276*0,165=0,3993 кг/м3

Теплоемкость смеси дымовых газов (путем интерполяции ):

 = (1,1276*0,64+1,1898*0,175+2,269*0,165)*103 = 1,3043*103 Дж/кг/К

Кинематическая вязкость смеси дымовых газов (путем интерполяции ):

 = (93,256*0,64+57,634*0,175+120,668*0,165)*10-6 = 89,68*10-6 м2/с

Расчетный коэффициент температуропроводности смеси газов:

Критерий Прандтля:

 =  = 0,715

Reг =  = 16678,91(1) = 0.035* * Reг0,8 * ( = 0,035*0,7320,33 *16678,910,8 * = 75,8

=  =  = 34,88 Вт/м2К.

Определение степени черноты излучающих трехатомных компонентов при температуре tст. [3]

tст = 594  ̊С

 = 0,1302

 = 0,4545

Степень черноты смеси газов при температуре tст.

Определение коэффициента теплоотдачи излучением при температуре tст.

 = 82,99 Вт/м²К.

 Вт/м²К.

Принимая за нулевое приближение плотность теплового потока q = 10372,39Вт/м, определяем коэффициент теплоотдачи от накипи к пару:

𝛼_П = 3,14q^0.7 *  = 3,14*10372,39^0,7 *18,2^0,15 = 3140,91 Вт/м²К.

Определение термических сопротивлений:

 = 0,0148 м  ̊С/Вт

 =0,00063  м  ̊С/Вт

 = 0.0215 м  ̊С/Вт

 = 0.000442 м  ̊С/Вт

 = 0,0147 м  ̊С/Вт

∑R(1) = 0.0148+0.00063+0.0215+0.000442+0.0147=0.0521 м  ̊С/Вт

Количество тепла, проходящее через 1 метр цилиндрической поверхности.

q =  =  = 10887,859 Вт/м

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

 =  = 21753,27 Вт/м².

=  =  = 19031,13 Вт/м².

Расчет температур на границах между слоями стенки:

 = 775-10887,859*0,0148 = 613,55  ̊С

 = 775-10887,859*(0.0148+0.0215) = 379,53  ̊С

 = 775-10887,859*(0.0148+0.0215+0.000442)=

= 374,72  ̊С

 = 775-10887,859 * (0.0148 +0.0215 +0.000442 +0,0147))  =  214,43  ̊С

Определение погрешности вычислений нулевого приближения:

𝜀 =  * 100% =   *100% = 4,7 ˃2%

2.6 Выполнение второго приближения:

Принимаем за начальные температуры первое приближение:

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

q = 10887,859 Вт/м.

Определение среднеинтегральной теплопроводности металла трубы: (3)

 =  =  = 377,13    =  37,686Вт/мК

Расчет теплофизических характеристик газов для tст = 614   [1]

Теплопроводность смеси дымовых газов (путем интерполяции, аналогично нулевого приближения):

 Вт/мК

Плотность смеси дымовых газов (путем интерполяции ):

 0,3864*0,64+0,6092*0,175+0,2229*0,165=0,3906 кг/м³

Теплоемкость смеси дымовых газов (путем интерполяции ):

 = (1,1321*0,64+1,1963*0,175+2,2831*0,165)*103 = 1,310*103 Дж/кг/К

Кинематическая вязкость смеси дымовых газов (путем интерполяции ):

= (96,918*0,64+60,134*0,175+125,5*0,165)*10-6 = 93,258*10-6 м2/с

Расчетный коэффициент температуропроводности смеси газов:

 = = 12,981* 10-5 м2/с

Критерий Прандтля:

 =  = 0,718

Re_г =  = 16678,91

Nu⁽²⁾ = 0.035* * Re_г^0,8 * ( = 0,035*0,732^0,33 *16678,91^0,8 * = 75,72

=  =  = 34,851 Вт/м²К.

Определение степени черноты излучающих трехатомных компонентов при температуре tст. [3]

t_ст = 614  ̊С

 = 0,1312

 = 0,4519Степень черноты смеси газов при температуре tст.

Определение коэффициента теплоотдачи излучением при температуре tст.

 = 84,56 Вт/м²К.

 Вт/м²К.

Принимая за нулевое приближение плотность теплового потока q = 10887,859Вт/м, определяем коэффициент теплоотдачи от накипи к пару:

𝛼_П = 3,14q^0.7 *  = 3,14*10887,859^0,7 *18,2^0,15 = 3249,377 Вт/м²К.

Определение термических сопротивлений:

 = 0,0146 м  ̊С/Вт

 =0,000615  м  ̊С/Вт

 = 0.0215 м  ̊С/Вт

 = 0.000445 м  ̊С/Вт

 = 0,0147 м  ̊С/Вт

∑R⁽²⁾ = 0.0146+0.000615+0.0215+0.000445+0.0147=0.0519 м  ̊С/Вт

Количество тепла, проходящее через 1 метр цилиндрической поверхности.

q =  =  = 10931,56 Вт/м

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

 =  = 21753,27 Вт/м².

=  =  = 19031,12 Вт/м².

Расчет температур на границах между слоями стенки:

 = 775-10931,56*0,0146 = 496,42  ̊С

 = 775-10931,56*(0.0146+0.0215) = 380,03  ̊С

 = 775-10931,56*(0.0146+0.0215+0.000445)=

= 375,16  ̊С

 = 775-10931,56 * (0.0146 +0.0215 +0.000445 +0,0147))  =  214,22  ̊С

Определение погрешности вычислений нулевого приближения:

𝜀 =  * 100% =   *100% = 0,4 ˃2%

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

Заключение

В результате выполнения курсовой работы были рассчитаны следующие параметры:

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

Тепловой поток через цилиндрическую стенку.

Тепловой поток через внешнюю поверхность труб.

Тепловой поток через внутреннюю поверхность труб.

Температуры на границах слоев.

Перечень ссылок

Казанцев Е.И. Промышленные печи: - М.: Металлургия 1975 - 368, с.

Ривкин С.Л. Термодинамические свойства воды и водяного пара: - М.: Энергоатомиздат 1984 - 80 с.

Кутателадзе С.С. Справочник по теплопередаче: - М.: Госэнергоиздат 1958 - 415 с.

Приложения

Таблица 1 Физические свойства дымовых газов

Таблица 2

Таблица 3

Таблица 4 Физические свойства водяного пара на линии насыщения