Курсовая работа: Расчет охладителя конденсата пара

Министерство образования Российской Федерации

Архангельский государственный технический университет

Факультет промышленной энергетики, III-2

Кафедра промышленной теплоэнергетики

Курсовой проект по курсу:

«Тепломассобменное оборудование предприятий»

Расчет охладителя конденсата пара

Архангельск 2007

1. Краткое описание конструкции аппарата

Установка охладителя конденсата греющего пара какого-либо подогревателя, приводит к уменьшению количества отбираемого из турбины пара на этот подогреватель и соответствующему увеличению расхода пара из отбора с меньшим давлением. Это несколько увеличивает тепловую экономичность установки. С другой стороны, увеличивается стоимость устанавливаемого оборудования. Таким образом, определение поверхности нагрева охладителя (минимальной разности температур теплообменивающихся сред), как, впрочем, и поверхности нагрева собственно подогревателя, является технико-экономической задачей.

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

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

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

Схема теплообменного аппарата приведена на рисунке 1.1. горизонтальный, двухходовой по конденсату пара и воды на ХВО. Движение потоков в охладителе применяется противоточное. Конденсата движется в межтрубном пространстве, вода на ХВО-в трубном.

Горячий агент (конденсат) поступает в обечайку (1) через входной патрубок (9) и заполняет межтрубное пространство. Выводится конденсата через выходной патрубок(10).Вода на ХВО входит в водяную камеру (4) через входной патрубок (7), проходит по теплообменным трубкам (3), совершает поворот и возвращается обратно в водяную камеру и выводится через выходной патрубок (8).Необходимое число ходов в аппарате создаётся за счёт перегородки в водяной камере (6) и перегородки в межтрубном пространстве (5).

2. Расчет недостающих параметров в аппарате

Определяем теплофизические свойства теплоносителей по их средним температурам.

Средняя температура греющего теплоносителя:

;

где ^o C находим по таблице 12 [1] при Р=0,4 МПа;

^o C,

^o C.

Средняя температура нагреваемого теплоносителя:

;

где ^o C;

^o C,

^o C.

По таблице 11 [1] определяем теплофизические свойства теплоносителей и сводим их в таблицу 1.

Таблица 1. Теплофизические свойства теплоносителей

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

_,

где Q – тепловая нагрузка аппарата (тепловая производительность), кВт;

G₁ - расход греющего теплоносителя (конденсат пара), кг/с;

G₂ – расход нагреваемого теплоносителя (вода на ХВО), кг/с;

, - теплоемкости греющего и нагреваемого теплоносителей соответственно, взятые по средним температурам, кДж/(кг K);

, - температуры на входе и выходе из аппарата греющего теплоносителя соответственно, ^o C;

,- температуры на входе и выходе из аппарата нагреваемого теплоносителя соответственно, ^o C;

- коэффициент удержания теплоты изоляцией;

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

,

где ,

,

^o C,

^o C,

принимаем равным 0,98,

.

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

,

где,

,

^o C,

^o C,

Определяем среднелогарифмическую разность температур:

,

где - наибольшая и наименьшая разница температур.

^o C

^o C

из этих двух значений выбираем наибольшее и наименьшее

^o C

^o C

^o C

Определяем температуру стенки

,

^o C

по этой температуре определяем Pr_ст =2,851 по таблице 11 [1].

3. Сравнение поверхностей теплообмена по энергетическим характеристикам

Необходимо выбрать оптимальную поверхность теплообмена из гладких труб (; ) и труб с кольцевыми выступами (d/D=0,94; t/D=0,5 ст. 375 [2]).

Рисунок 2. Труба с кольцевыми выступами и гладкая труба

Расчет теплообменника с гладкими трубами.

Принимаем скорость нагреваемого теплоносителя , равной 2 м/с.

Необходимое сечение канала можно определить из уравнения сплошности:

,

где G ₂ – расход греющего теплоносителя, кг/с;

-принятая скорость нагреваемого теплоносителя, м/с;

- плотность греющего теплоносителя, взятая по средней температуре, ^o C.

Тогда необходимое сечение канала будет:

,

где G ₂ =121,5 кг/с;

м/с;

кг/м³

м² .

Определяем приблизительное число труб в одном ходу:

,

где м² ;

м, внутренний диаметр труб.

шт.

Найдем общее число трубок:

,

где число ходов в аппарате.

шт.

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

Точное число трубок определяем исходя из табл. 23.1 [6] шт.

Окончательное число труб принимаем:

,

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

шт.

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

,

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

мм

коэффициент заполнения площади трубной решетки трубами (зависит от числа ходов по трубному пространству), т. к. то .

мм

Конечно диаметр принимаем по табл. 15.1 [6] D_вн =650 мм.

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

,

где количество труб в одном ходу

,

шт.

м/с.

Определяем площадь межтрубного сечения для греющего теплоносителя:

,

где мм толщина перегородки в межтрубном пространстве, принятая конструктивно.

м² .

Определяем скорость греющего теплоносителя в межтрубном пространстве:

,

м/с.

Определяем смоченный периметр по греющему теплоносителю:

,

мм.

Определяем эквивалентный диаметр по греющему теплоносителя:

,

мм.

В результате перерасчета задаемся другой температурой стенки ^o C по этой температуре определяем Pr_ст =3,848 по таблице 11 [1].

Определим число Рейнольдса для нагреваемого теплоносителя:

,

.

Определим число Рейнольдса для греющего теплоносителя:

,

.

Определим числа Нуссельта для греющего и нагреваемого теплоносителей по формуле Михеева, так как режим течения турбулентный :

,

где - число Прандтля, принимается по таблице 1;

;

.

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

,

Вт/(м² К).

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

,

Вт/(м² К).

Проверяем температуру стенки:

,

^o C.

Полученная температура незначительно отличается от предварительно принятой.

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

,

где Вт/(мК) коэффициент теплопроводности трубки по табл. 7 [1],

м, толщина стенки трубки,

коэффициент загрязнения.

Вт/(м² К).

Определим расчетную поверхность теплообмена аппарата;

,

м² .

Определим активную длину трубок:

,

где средний диаметр,

м.

м.

Определим конструктивность аппарата:

,

условие соблюдается.

Принимаем скорость нагреваемого теплоносителя , равной 1 м/с.

Необходимое сечение канала можно определить из уравнения сплошности:

,

где G ₂ – расход греющего теплоносителя, кг/с;

-принятая скорость нагреваемого теплоносителя, м/с;

- плотность греющего теплоносителя, взятая по средней температуре, ^o C.

Тогда необходимое сечение канала будет:

,

где G ₂ =121,5 кг/с;

м/с;

кг/м³

м² .

Определяем приблизительное число труб в одном ходу:

,

где м² ;

м, внутренний диаметр труб.

шт.

Найдем общее число трубок:

,

где число ходов в аппарате.

шт.

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

Точное число трубок определяем исходя из табл. 23.1 [6] шт.

Окончательное число труб принимаем:

,

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

шт.

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

,

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

мм

коэффициент заполнения площади трубной решетки трубами (зависит от числа ходов по трубному пространству), т.к. то .

мм

Конечно диаметр принимаем по табл. 15.1 [6] D_вн =900 мм.

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

,

где количество труб в одном ходу

,

шт.

м/с.

Определяем площадь межтрубного сечения для греющего теплоносителя:

,

где мм толщина перегородки в межтрубном пространстве, принятая конструктивно.

м² .

Определяем скорость греющего теплоносителя в межтрубном пространстве:

,

м/с.

Определяем смоченный периметр по греющему теплоносителю:

,

мм.

Определяем эквивалентный диаметр по греющему теплоносителя:

,

мм.

Определим число Рейнольдса для нагреваемого теплоносителя:

,

.

Определим число Рейнольдса для греющего теплоносителя:

,

.

Определим числа Нуссельта для греющего и нагреваемого теплоносителей по формуле Михеева, так как режим течения турбулентный :

,

где - число Прандтля, принимается по таблице 1;

;

.

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

,

Вт/(м² К).

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

,

Вт/(м² К).

Проверяем температуру стенки:

,

^o C.

Полученная температура незначительно отличается от предварительно принятой.

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

,

гдеВт/(мК)) коэффициент теплопроводности трубки по табл. 7 [1],

м, толщина стенки трубки.

коэффициент загрязнения

Вт/(м² К).

Определим расчетную поверхность теплообмена аппарата;

,

м² .

Определим активную длину трубок:

,

где средний диаметр,

м.

м.

Определим конструктивность аппарата:

,

условие соблюдается.

Гидравлический расчет для гладких труб.

Для скорости нагреваемого теплоносителя , равной 1 м/с:

Определим полную длину трубок:

,

где м толщина трубной решетки принятая конструктивно.

м высота выступа трубок принятая конструктивно.

м.

При турбулентном режиме движения воды коэффициент трения по трубному пространству находим по формуле Блазиуса:

,

.

Определим потери давления на трение по трубному пространству:

,

где количество ходов по трубному пространству.

Па.

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

,

где сумма коэффициентов местных сопротивлений, где

x_вх – коэффициент местного сопротивления при входе потока в камеру,

принимаем x_вх =1,5;

x_вых – коэффициент местного сопротивления при выходе потока из камеры, принимаем x_вых =1,5;

x_п – коэффициент местного сопротивления при повороте потока на 180°, принимаем x_п =2,5;

x_вх.тр – коэффициент местного сопротивления при входе потока в трубки, принимаем x_вх =0,5;

x_вых.тр – коэффициент местного сопротивления при выходе потока из трубок, принимаем x_вых =1,0;

Па.

Определим потерю давления по трубному пространству:

,

Па.

Определим мощность, потребляемую насосом для перемещения воды по трубному пространству:

,

Вт.

При турбулентном режиме движения воды коэффициент трения по межтрубному пространству находим по формуле Блазиуса:

,

.

Определим потери давления на трение по межтрубному пространству:

,

где количество ходов по межтрубному пространству.

Па.

Определим потери давления на местные сопротивления в аппарате по межтрубному пространству:

,

где сумма коэффициентов местных сопротивлений

Па.

Определим потерю давления по межтрубному пространству:

,

Па.

Определим мощность, потребляемую насосом для перемещения конденсата по межтрубному пространству:

,

Вт.

Для скорости нагреваемого теплоносителя , равной 2 м/с:

Определим полную длину трубок:

,

где м толщина трубной решетки принятая конструктивно.

м высота выступа трубок принятая конструктивно.

м.

При турбулентном режиме движения воды коэффициент трения по трубному пространству находим по формуле Блазиуса:

,

.

Определим потери давления на трение по трубному пространству:

,

где количество ходов по трубному пространству.

Па.

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

,

где сумма коэффициентов местных сопротивлений, где

x_вх – коэффициент местного сопротивления при входе потока в камеру, принимаем x_вх =1,5;

x_вых – коэффициент местного сопротивления при выходе потока из камеры, принимаем x_вых =1,5;

x_п – коэффициент местного сопротивления при повороте потока на 180°, принимаем x_п =2,5;

x_вх.тр – коэффициент местного сопротивления при входе потока в трубки, принимаем x_вх =0,5;

x_вых.тр – коэффициент местного сопротивления при выходе потока из трубок, принимаем x_вых =1,0;

Па.

Определим потерю давления по трубному пространству:

,

Па.

Определим мощность, потребляемую насосом для перемещения воды по трубному пространству:

,

Вт.

При турбулентном режиме движения воды коэффициент трения по межтрубному пространству находим по формуле Блазиуса:

,

.

Определим потери давления на трение по межтрубному пространству:

,

где количество ходов по межтрубному пространству.

Па.

Определим потери давления на местные сопротивления в аппарате по межтрубному пространству:

,

где сумма коэффициентов местных сопротивлений

Па.

Определим потерю давления по межтрубному пространству:

,

Па.

Определим мощность, потребляемую насосом для перемещения конденсата по межтрубному пространству:

,

Вт.

Расчет теплообменника с кольцевыми выступами.

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

Для определения среднего коэффициента теплоотдачи (для капельных жидкостей) при d/D=0,935 и t/D=0,5 используется уравнение подобия для критерия Нуссельта следующего вида (формула 6.33 [2]):

,

где .

В результате перерасчета задаемся другой температурой стенки ^o C по этой температуре определяем Pr_ст =4,865 по таблице 11 [1].

Определим число Нуссельта для нагреваемого теплоносителя:

,

.

Определим число Нуссельта для греющего теплоносителя по формуле Михеева, так как режим течения турбулентный :

,

.

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

,

Вт/(м² К).

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

,

Вт/(м² К).

Проверяем температуру стенки:

,

^o C.

Полученная температура незначительно отличается от предварительно принятой.

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

,

гдеВт/(мК) теплопроводность латуни,

м

коэффициент загрязнения.

Вт/(м² К).

Определим расчетную поверхность теплообмена аппарата;

,

м² .

Определим активную длину трубок:

,

где средний диаметр,

м.

м.

Определим конструктивность аппарата:

,

условие соблюдается.

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

Определим число Нуссельта для нагреваемого теплоносителя:

,

.

Определим число Нуссельта для греющего теплоносителя по формуле Михеева, так как режим течения турбулентный :

,

.

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

,

Вт/(м² К).

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

,

Вт/(м² К).

Проверяем температуру стенки:

,

^o C.

Полученная температура незначительно отличается от предварительно принятой.

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

,

где Вт/(мК) теплопроводность латуни,

м

коэффициент загрязнения.

Вт/(м² К).

Определим расчетную поверхность теплообмена аппарата;

,

м² .

Определим активную длину трубок:

,

где средний диаметр,

м.

м.

Определим конструктивность аппарата:

,

условие соблюдается.

Гидравлический расчет для кольцевых выступов.

Для скорости нагреваемого теплоносителя , равной 1 м/с:

Определим полную длину трубок:

,

где м толщина трубной решетки принятая конструктивно.

м высота выступа трубок принятая конструктивно.

м.

Коэффициент трения вычисляем по формуле 6.35 [2]:

,

где

.

Определим потери давления на трение по трубному пространству:

,

где количество ходов по трубному пространству.

Па.

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

,

где сумма коэффициентов местных сопротивлений, где

x_вх – коэффициент местного сопротивления при входе потока в камеру, принимаем x_вх =1,5;

x_вых – коэффициент местного сопротивления при выходе потока из камеры, принимаем x_вых =1,5;

x_п – коэффициент местного сопротивления при повороте потока на 180°, принимаем x_п =2,5;

x_вх.тр – коэффициент местного сопротивления при входе потока в трубки, принимаем x_вх =0,5;

x_вых.тр – коэффициент местного сопротивления при выходе потока из трубок, принимаем x_вых =1,0;

Па.

Определим потерю давления по трубному пространству:

,

Па.

Определим мощность, потребляемую насосом для перемещения воды по трубному пространству:

,

Вт.

При турбулентном режиме движения воды коэффициент трения по межтрубному пространству находим по формуле Блазиуса:

,

.

Определим потери давления на трение по межтрубному пространству:

,

где количество ходов по межтрубному пространству.

Па.

Определим потери давления на местные сопротивления в аппарате по межтрубному пространству:

,

где сумма коэффициентов местных сопротивлений

Па.

Определим потерю давления по межтрубному пространству:

,

Па.

Определим мощность, потребляемую насосом для перемещения конденсата по межтрубному пространству:

,

Вт.

Для скорости нагреваемого теплоносителя , равной 2 м/с:

Определим полную длину трубок:

,

где м толщина трубной решетки принятая конструктивно.

м высота выступа трубок принятая конструктивно.

м.

Коэффициент трения вычисляем по формуле 6.35 [2]:

,

где

Определим потери давления на трение по трубному пространству:

,

где количество ходов по трубному пространству.

Па.

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

,

где сумма коэффициентов местных сопротивлений, где

x_вх – коэффициент местного сопротивления при входе потока в камеру, принимаем x_вх =1,5;

x_вых – коэффициент местного сопротивления при выходе потока из камеры, принимаем x_вых =1,5;

x_п – коэффициент местного сопротивления при повороте потока на 180°, принимаем x_п =2,5;

x_вх.тр – коэффициент местного сопротивления при входе потока в трубки, принимаем x_вх =0,5;

x_вых.тр – коэффициент местного сопротивления при выходе потока из трубок, принимаем x_вых =1,0;

Па.

Определим потерю давления по трубному пространству:

,

Па.

Определим мощность, потребляемую насосом для перемещения воды по трубному пространству:

,

Вт.

При турбулентном режиме движения воды коэффициент трения по межтрубному пространству находим по формуле Блазиуса:

,

.

Определим потери давления на трение по межтрубному пространству:

,

где количество ходов по межтрубному пространству.

Па.

Определим потери давления на местные сопротивления в аппарате по межтрубному пространству:

,

где сумма коэффициентов местных сопротивлений

Па.

Определим потерю давления по межтрубному пространству:

,

Па.

Определим мощность, потребляемую насосом для перемещения конденсата по межтрубному пространству:

,

Вт.

Сравнение поверхностей теплобмена по энергетической эффективности.

Для гладких труб при скорости м/с:

Определим удельные затраты мощности на прокачку:

,

где суммарная мощность на прокачку обоих теплоносителей

Вт

Вт/м² К.

Определим коэффициент энергетической эффективности:

,

Для гладких труб при скорости м/с:

Определим удельные затраты мощности на прокачку:

,

где суммарная мощность на прокачку обоих теплоносителей

Вт

Вт/м² К.

Определим коэффициент энергетической эффективности:

,

Т.к. графики строим в логарифмических координатах то:

Для труб с кольцевыми выступами при скорости м/с:

Определим удельные затраты мощности на прокачку:

,

где суммарная мощность на прокачку обоих теплоносителей

Вт

Вт/м² К.

Определим коэффициент энергетической эффективности:

,

Для труб с кольцевыми выступами при скорости м/с:

Определим удельные затраты мощности на прокачку:

,

где суммарная мощность на прокачку обоих теплоносителей

Вт

Вт/м² К.

Определим коэффициент энергетической эффективности:

,

Т.к. графики строим в логарифмических координатах то:

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

Из этого графика нельзя точно судить о том, какие трубки эффективнее т.к. линии практически сливаются, поэтому проведем сравнение по другим характеристикам.

Сравнение поверхностей теплообмена по габаритной характеристике.

Для гладких труб при скорости м/с:

Определим коэффициент компактности:

,

где объем

м³ .

м² /м³ .

Для гладких труб при скорости м/с:

Определим коэффициент компактности:

,

где объем

м³ .

м² /м³ .

Для труб с кольцевыми выступами при скорости м/с:

Определим коэффициент компактности:

,

где объем

м³ .

м² /м³ .

Для труб с кольцевыми выступами при скорости м/с:

Определим коэффициент компактности:

,

где объем

м³ .

м² /м³ .

Найдем отношение для гладких трубок при скорости м/с:

;

Теперь найдем логарифм от этого соотношения:

Найдем отношение для гладких трубок при скорости м/с:

;

Теперь найдем логарифм от этого соотношения:

Найдем отношение для трубок с кольцевыми выступами при скорости м/с:

;

Теперь найдем логарифм от этого соотношения:

Найдем отношение для трубок с кольцевыми выступами при скорости м/с:

;

Теперь найдем логарифм от этого соотношения:

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

Определим на сколько гладкие трубки эффективнее трубок с кольцевыми выступами:

проекция на ось абсцисс для гладких трубок.

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

100%

100%=12,9%.

Гладкие трубки выгоднее.

4. Тепловой расчет аппарата

Принимаем скорость нагреваемого теплоносителя , равной 1,5 м/с.

Необходимое сечение канала можно определить из уравнения сплошности:

,

где G ₂ – расход греющего теплоносителя, кг/с;

-принятая скорость нагреваемого теплоносителя, м/с;

- плотность греющего теплоносителя, взятая по средней температуре, ^o C.

Тогда необходимое сечение канала будет:

,

где G ₂ =121,5 кг/с;

м/с;

кг/м³

м² .

Определяем приблизительное число труб в одном ходу:

,

где м² ;

м, внутренний диаметр труб.

шт.

Найдем общее число трубок:

,

где число ходов в аппарате.

шт.

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

Точное число трубок определяем исходя из табл. 23.1 [6] шт.

Окончательное число труб принимаем:

,

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

шт.

данные расчеты подтверждает компоновка приложение 1.

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

,

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

мм

коэффициент заполнения площади трубной решетки трубами (зависит от числа ходов по трубному пространству), т. к. то .

мм

Конечно диаметр принимаем по табл. 15.1 [6] D_вн =800 мм.

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

,

где количество труб в одном ходу

,

шт.

м/с.

Определяем площадь межтрубного сечения для греющего теплоносителя:

,

где мм толщина перегородки в межтрубном пространстве, принятая конструктивно.

м² .

Определяем скорость греющего теплоносителя в межтрубном пространстве:

,

м/с.

Определяем смоченный периметр по греющему теплоносителю:

,

мм.

Определяем эквивалентный диаметр по греющему теплоносителя:

,

мм.

В результате перерасчета задаемся другой температурой стенки ^o C по этой температуре определяем Pr_ст =3,848 по таблице 11 [1].

Определим число Рейнольдса для нагреваемого теплоносителя:

,

.

Определим число Рейнольдса для греющего теплоносителя:

,

.

Определим числа Нуссельта для греющего и нагреваемого теплоносителей по формуле Михеева, так как режим течения турбулентный :

,

где - число Прандтля, принимается по таблице 1;

;

.

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

,

Вт / (м² К).

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

,

Вт/(м² К).

Проверяем температуру стенки:

,

^o C.

Полученная температура незначительно отличается от предварительно принятой.

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

,

где Вт/(мК) коэффициент теплопроводности трубки по табл. 7 [1],

м, толщина стенки трубки,

коэффициент загрязнения.

Вт/(м² К).

Определим расчетную поверхность теплообмена аппарата;

,

м² .

Определим активную длину трубок:

,

где средний диаметр,

м.

м.

Определим конструктивность аппарата:

,

условие соблюдается.

5. Конструктивный расчет аппарата

Гидравлический расчет патрубков.

Определим внутренний диаметр патрубков.

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

,

где площадь поперечного сечения патрубка, по которому движется теплоноситель, м² ;

плотность теплоносителя в патрубке, кг/м³ ;

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

Тогда площадь внутреннего поперечного сечения патрубка находим по формуле:

,

Внутренний диаметр патрубка определим по формуле:

,

Определим диаметр патрубка для греющего теплоносителя (если теплоноситель не меняет свое агрегатное состояние, то диаметры труб на входе и выходе равны):

м.

Расчётный диаметр патрубка округляем до ближайшего стандартного значения. Принимаем трубы стальные прямошевные для прямых участков трубопроводов по табл. 20.9 [6]. Условный диаметр d_у = 200 мм; наружный диаметр d_н = 219 мм; толщина стенки трубы S = 6 мм, материал сталь 10Г2.

Определим диаметр патрубка для нагреваемого теплоносителя (если теплоноситель не меняет свое агрегатное состояние, то диаметры труб на входе и выходе равны):

м.

Расчётный диаметр патрубка округляем до ближайшего стандартного значения. Принимаем трубы стальные прямошевные для прямых участков трубопроводов по табл. 20.9 [6]. Условный диаметр d_у = 300 мм; наружный диаметр d_н = 325 мм; толщина стенки трубы S = 8 мм, материал стальВСт3.

Подбор фланцев.

Стандартные фланцы принимаются по условному диаметру d_y и давлению, и их размеры представлены в таблице 5.1. Патрубок с присоединенным к нему фланцем показан на рис. 5.1.

Размеры фланца выбираются по табл. 20.9 [6] (Тип а).

Таблица 5.1

Рис. 5.1. Фланец плоский приварной

Определение вылетов фланцевых штуцеров.

Каждый патрубок в зависимости от его размера имеет определённую длину вылета, которую назначаем по табл. 25.1 [6].

Для патрубка конденсата пара: l = 130 мм.

Для патрубка воды на ХВО: l = 140 мм.

Крепление труб в трубной решетке.

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

Рис. 5.2. Развальцовка с отбортовкой труб в трубной решётке.

l_в – высота вальцовки;

S-толщина трубной решётки;

h_в –высота выступа;

d_о –диаметр отверстия под трубку в решётке;

d_н –наружный диаметр трубки;

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

Проверка патрубков на прочность.

Условие прочности:

S_прин >S_исп ,

где S_прин – принятая толщина патрубка, мм;

S_исп – исполнительная толщина стенки штуцера, мм, определяемая по формуле:

S_исп =S_р +С, мм

где S_p – расчетная толщина стенки штуцера, определяемая по формуле:

,

С – прибавка к расчетной толщине штуцера, определяемая по формуле:

С = C₁ + С₂ =1+3 =4 мм.

,

где t – время эксплуатации теплообменника, лет; t = 10 лет;

u_к – скорость коррозии, мм/год, u_к = 0,3 мм/год;

мм

Патрубок для конденсата пара.

При температуре °С, 154,76 МПа из табл1.3 [7], тогда

, мм

Выбираем стандартную толщину листа из справочника [4] табл. 2.22:

S_исп =0,3+4=4,3 мм.

Т.к. S_прин = 6 мм, то S_исп < S_прин – условие прочности выполняется.

=S– ()=6–4,3=1,7 мм

Патрубок для конденсата пара.

При температуре °С, 138,88 МПа из табл1.3 [7], тогда

, мм

S_исп =0,7+4=4,7 мм.

Т.к. S_прин = 8 мм, то S_исп < S_прин – условие прочности выполняется.

=S– ()=8–4,7=3,3 мм.

Расчет толщины стенки цилиндрической обечайки.

Расчет допускаемых напряжений.

Для изготовления обечайки теплообменника выбираем сталь 10Г2 ГОСТ 1050–74.

Допускаемое напряжение [s], МПа вычисляем по формуле:

, МПа

где s* – нормативное допускаемое напряжение, МПа; при расчетной температуре

t_H = 143,62 s* = 154,76 МПа, по табл. 1.3. [7];

h – поправочный коэффициент; h= 1;

[s] = 154,76*1 =154,76 МПа.

Расчет толщины стенки обечайки.

Расчет выполняем по формуле:

,

где S_p – расчетная толщина стенки обечайки, мм;

Р_р – расчетное давление, МПа; Р_р = P₁ = 0,4 МПа;

D_BH – внутренний диаметр обечайки, мм; D_ВН = 800 мм;

j – коэффициент прочности сварных соединений; j = I.;

мм.

Исполнительную толщину стенки обечайки определяем по формуле:

S = S_p + С, мм

где С – прибавка к расчетной толщине стенки, м, вычисляемая по формуле:

, мм

где C₁ – прибавка для компенсации минусовых допусков на толщину листа, мм; С₁ = 1 мм; С₂ =3 мм; =1,03+1+3=5,03 мм

Выбираем стандартную толщину листа из справочника [4] табл. 2.22: S=8 мм.

Расчет толщины стенки цилиндрической части водяной камеры

Используем сталь 10Г2 ГОСТ 1050–74

t = 35 °С 138,875 МПа Р_раб =Р₂ =0,6МПа

мм

S==1,73+1+3+0,7=5,73 мм

Выбираем стандартную толщину листа из справочника [4] табл. 2.22:

S=8 мм. Длину выбираем исходя из формулы:

,

где ширина укрепляемой зоны

мм

мм.

Расчет фланцев для обечайки

В качестве фланца проектируемого аппарата выбираем фланец плоский приварной. Материал фланца: сталь 20. Базовые размеры выбранных фланцев сведены в таблицы 5.2 Эскиз фланцев для аппаратов с указанием базовых размеров представлен на рисунке 5.2

Таблица 5.2 Базовые размеры плоского приварного фланца для аппарата

Рис. 5.3. Эскиз плоского приварного фланца для аппарата

Расчет толщин днищ и крышек цилиндра теплообменного аппарата.

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

Выбираем эллиптические днища табл. 16.1 [6], эскиз которых представлен на рис. 5.3.

Рис. 5.4. Эллиптическое днище.

Расчетную толщину эллиптических отбортованных днищ, определяем по формуле:

, мм

где - рабочее давление, МПа;

R– радиус кривизны в вершине днища, мм; R = D_BH

[s] – допускаемое напряжение для материала днища при его рабочей температуре

[s] = 138,88 МПа;

мм.

Выбираем стандартную толщину днища: S=8 мм.

Высоту выпуклой части крышки h_в , м находим по формуле:

H_в = 0,25 D_вн = 0,25×800=200 мм.

Таблица 5.3.

Укрепление отверстий.

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

Чтобы скомпенсировать потерю прочности, существуют 2 принципиальных способа:

1. Повышение толщины стенки всей оболочки аппарата.

2. Повышение толщины оболочки только в некоторой зоне, расположенной вокруг отверстия. Последний способ более целесообразен, так как требует меньших затрат металла. Нормы и методы расчета на прочность, укрепление отверстий установлены ГОСТ 24755–81.

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

Рис. 5.5. Отверстие, укрепленное с помощью кольца.

В этой схеме (на рис. 5.4) толщина укрепляющего кольца примерно равна толщине оболочки. Кольцо имеет сигнальное отверстие резьбой М10, которое предназначено для испытания сварочных швов на плотность. В это отверстие с резьбой при испытании подводится воздух давлением 0,6 МПа.

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

где S- исполнительная толщина укрепляемого элемента, мм;

S_р =S- c-расчетная толщина оболочки, мм;

D_R - расчетный диаметр оболочки, для оболочки D_R =D_в ,

c₂ - прибавка на коррозию;

Для обечайки:

мм

Если , то отверстия можно не уплотнять.

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

,

где S_p -расчетная толщина оболочки;

d_в.шт. -внутренний диаметр штуцера, мм;

– длина штуцера, мм;

– расчетная толщина стенки штуцера, мм;

-исполнительная толщина стенки штуцера, мм;

c₂ - прибавка на коррозию, мм;

c₃ - прибавка для округления расчетной толщины штуцера до размера по сортаменту, мм;

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

Из расчетов следует, что отверстия нужно укреплять.

Диаметр укрепляющего кольца находится по формуле:

,

где d_н.шт. - наружный диаметр штуцера, м;

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

м.

Действительную толщину укрепляющего кольца определяют из неравенства:

,

где S_K -толщина укрепляющего кольца, мм;

D_K -наружный диаметр кольца, мм;

d– диаметр укрепляемого отверстия, мм;

S-исполнительная толщина оболочки (обечайки), мм;

d₀ ^пред - предельный допускаемый диаметр, мм;

Из этого неравенства следует:

Таким образом, для патрубков конденсата пара:

мм

Исполнительную толщину кольца принимаем 8 мм.

Рассчитаем диаметр укрепляющего кольца для патрубков воды на ХВО:

м.

Действительную толщину укрепляющего кольца определяют из неравенства:

,

где S_K - толщина укрепляющего кольца, мм;

D_K - наружный диаметр кольца, мм;

d– диаметр укрепляемого отверстия, мм;

S- исполнительная толщина оболочки (обечайки), мм;

d₀ ^пред - предельный допускаемый диаметр, мм;

Из этого неравенства следует:

Таким образом, для патрубков воды на ХВО:

мм

Исполнительную толщину кольца принимаем 8 мм.

Расчет трубной решетки

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

Толщину трубной решетки определяем:

,

где средний диаметр под прокладку,

мм

– допускаемое напряжение для материала решетки;

C=4 мм – прибавка к расчетной высоте;

коэффициент ослабления трубной решетки отверстиями под трубки;

диаметр отверстия под трубку;

м,

м.

м

Дополнительные условия проверки прочности трубной решетки.

Условие прочности безтрубной зоны:

где d_E – максимальный диаметр окружности, вписанной в безтрубную зону решетки меж ду обечайкой и трубками, м; принимаем по компоновке D_e = 10 мм;

Р_т – давление в трубном пространстве, МПа; Р_т = 0,6 МПа;

Р_м – давление в межтрубном пространстве, МПа; Р_м = 0,4МПа;

- допускаемое напряжение для материала решетки, МПа; =154,76МПа;

м

S³0,00418 м.

Толщина трубной решетки в сечении уплотнительной канавки должна быть не менее:

где (S – С) – расчетная толщина трубной решетки, м; (S – С) = 26 мм;

t_p – шаг разбивки трубной решетки t_p = 31,4 мм;

t_п – расстояние между ближайшими рядами труб, принимаем по компоновке t_п = 62,5 мм

b_п – ширина паза в трубной решетке; b_n = 10 мм;

d_o – диаметр отверстия под трубки в трубной решетке.

φ_р – коэффициент прочности трубной решетки,

мм.

Минимальная толщина трубной решетки, обеспечивающая надежность вальцовки труб:

где l_В – высота вальцовки, м;

[q] – допускаемая нагрузка, приходящаяся на единицу условной поверхности развальцовки и зависящая от способа развальцовки, МПа; для развальцовки с отбортовкой [q] = 30 МПа;

N_T – продольная сила, действующая в месте закрепления трубы в решетке, МН.

Так как рассчитываемый аппарат является аппаратом нежесткого типа, осевое усилие вычисляется по формуле:

Где Р=0,4 МПа – давление в межтрубном пространстве

мм

Минимальная толщина решетки под кольцевую прокладку.

,

где D_с.п. - средний диаметр прокладки, м;

допускаемое напряжение для материала решетки, МПа;

D_в - внутрений диаметр обечайки;

(м)

Условие выполняется.

Дополнительное условие:

мм

Все условия выполняются.

Расчет теплообменных трубок на прочность.

Расчетные напряжения в осевом направлении s_1т , МПа, в трубах составляют:

,

где S_T – толщина стенки трубки, м .

.

Расчетные напряжения в окружном направлении s_2т , МПа, в трубах составляют:

,

.

Условие статической прочности труб:

,

где [s]^Т – допускаемое напряжение для материала трубки, МПа, принимаем [s]^Т = 230МПа.

.

Таким образом, условие статической прочности труб выполняется.

Расчет прокладочной обтюрации

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

Для герметичности соединения обечайки с жидкостными крышками используем прокладки типа 1. Тип обтюрации – I-А. Обтюрация типа I-А изображена на рисунке 13.1. Расчетная сила осевого сжатия для прокладок типа I определяется по формуле:

,

где D _сп – средний диаметр прокладки, м ;

P – расчетное давление среды, воздействующей на фланцевое соединение, МПа,

K – коэффициент, зависящий от материала прокладки и ее конструкции, для паронита K =2,5 ;

b _экв – эффективная ширина уплотнения, м

,

b – ширина прокладки, м

.

.

.

МН.

Выбор линзового компенсатора.

Рисунок 5.7. Сварной линзовый компенсатор из двух штампованных полулинз

По давлению в межтрубном пространстве и наружному диаметру выбираем линзовый компенсатор по таблице 24.1 [6]. Таблица 2. Основные размеры линзового компенсатора.

Таблица 5.4

По конструктивным соображениям число линз принимаем

Расчет опор.

Установка аппаратов на фундаменты или на специальные несущие конструкции осуществляется в основном посредством опор. Непосредственно на фундаменты устанавливаются аппараты с плоским днищем. Опоры могут размещать или снизу аппарата или с боков и они жестко соединяются с аппаратом. Выбор типа опоры зависит от места установки теплообменного аппарата (в помещение или вне его), соотношения высоты к диаметру и массы аппарата.

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

Максимальное усилие G_max определяется по формуле:

Для расчета составляющих, входящих в формулу используем известное соотношение для i -го элемента:

где V_i – объем i -го элемента, м³

ρ _i – плотность материала элемента, кг/м³

Для стали: ρ_ст =7850 кг/м³

Для воды: ρ_в =1000 кг/м³

Для латуни: ρ_л =8430 кг/м³

Усилие от обечайки G_o определяем по следующей формуле:

где D _н – наружный диаметр обечайки, м

D _вн – внутренний диаметр обечайки, м

l – длина обечайки, м

g=9,81 м² /с

D_н =0,816 м

D_вн =0,8 м

l=3,7 м

Усилие, действующее на опору со стороны цилиндрических частей днищ:

где l_ц – высота цилиндрической части днища, м

конструктивно принимаем l_ц =0,2 м

Усилие от эллиптических днищ:

где m _дн – масса днища, кг

по таблице 16.2 [6] принимаем стандартную массу днища m _дн =49 кг

Усилие, действующее на опоры со стороны трубных решеток:

где D _н.р. – наружный диаметр трубной решетки, м

h – высота трубной решетки, м

d _н – наружный диаметр трубки в аппарате, м

D_н.р. =0,816 м

h=0,03 м

d_н =0,024 м

Усилие от фланцев на обечайке:

где h _ф – высота фланца, м

D _ф – наружный диаметр фланца, м

D _в – внутренний диаметр фланца, м

h _ф = 0,028 м

D _ф =0,93 м

D _в =0,8 м

Усилия, действующие со стороны теплообменных трубок, находятся из соотношения:

где n – количество трубок,

l _т – полная длина трубки, м

ρ_л – плотность латуни, кг/м³

d _н – наружный диаметр трубки, м

d _вн – внутренний диаметр трубки, м

Усилие, действующее со стороны патрубков греющего теплоносителя:

где l _п.г. – вылет патрубков греющего теплоносителя, м

D _н – наружный диаметр патрубка, м

D _вн – внутренний диаметр патрубка, м

Усилие, действующее со стороны патрубков греющего теплоносителя:

где l _п.г. – вылет патрубков нагреваемого теплоносителя, м

D _н – наружный диаметр патрубка, м

D _вн – внутренний диаметр патрубка, м

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

Усилие фланцев для патрубков:

Усилие, действующее на опоры со стороны фланцев для патрубков конденсата пара: так как их два, то

Н.

Усилие, действующее на опоры со стороны фланцев для патрубков воды на ХВО: так как их два, то

Н.

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

Н

Усилие, действующее на опоры со стороны всего аппарата:

Усилие, действующее со стороны воды, заполняющей весь аппарат:

Объем обечайки находится по формуле:

Емкость днищ определяется по формуле:

где h – высота цилиндрической части днища, м

V ^* – емкость эллиптического днища, м³

h =0,02 м

V ^* =0,0796м³

Объем цилиндрической камеры обечайки находится по формуле:

Суммарный объем воды в аппарате:

Вес воды, заполняющей аппарат:

Определяем усилие, действующее на опору со стороны всего аппарата, заполненного водой:

Количество опор составляет обычно от двух до четырех и нагрузка приходящаяся на одну опору равна:

где n – количество опор

Принимаем количество опор в аппарате равным п=2.

В качестве материала лапы назначаем ВСт3

В зависимости от величины нагрузки G по табличным данным выбираем стандартную опорную лапу.

Ориентируясь на таблицу 14.6 [7] принимаем опору по ГОСТ 26–1265–75

Таблица 5.5. Базовые размеры опоры аппарата

D н , мм	Q , кН	, мм	, мм	R , мм	L , мм	A , мм	, мм
800	80	8	14	442	740	500	400

6. Гидравлический расчет аппарата

Для скорости нагреваемого теплоносителя , равной 1,5 м/с:

Определим полную длину трубок:

,

где м толщина трубной решетки.

м высота выступа трубок.

м.

При турбулентном режиме движения воды коэффициент трения по трубному пространству находим по формуле Блазиуса:

,

.

Определим потери давления на трение по трубному пространству:

,

где количество ходов по трубному пространству.

Па.

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

,

где сумма коэффициентов местных сопротивлений, где

x_вх – коэффициент местного сопротивления при входе потока в камеру, принимаем x_вх =1,5;

x_вых – коэффициент местного сопротивления при выходе потока из камеры, принимаем x_вых =1,5;

x_п – коэффициент местного сопротивления при повороте потока на 180°, принимаем x_п =2,5;

x_вх.тр – коэффициент местного сопротивления при входе потока в трубки, принимаем x_вх =0,5;

x_вых.тр – коэффициент местного сопротивления при выходе потока из трубок, принимаем x_вых =1,0;

Па.

Определим потерю давления по трубному пространству:

,

Па.

Определим мощность, потребляемую насосом для перемещения воды по трубному пространству:

,

Вт.

При турбулентном режиме движения воды коэффициент трения по межтрубному пространству находим по формуле Блазиуса:

,

.

Определим потери давления на трение по межтрубному пространству:

,

где количество ходов по межтрубному пространству.

Па.

Определим потери давления на местные сопротивления в аппарате по межтрубному пространству:

,

где сумма коэффициентов местных сопротивлений

Па.

Определим потерю давления по межтрубному пространству:

,

Па.

Определим мощность, потребляемую насосом для перемещения конденсата по межтрубному пространству:

,

Вт.

7. Расчет тепловой изоляции

Тепловая изоляция необходима для предотвращения потерь тепла в окружающую среду и обеспечения безопасных условий эксплуатации теплообменных аппаратов. Изоляции подлежат все аппараты, температура стенки которых больше 50°С,а также аппараты с температурой меньше 0 °С. В качестве материала теплоизоляции выбираем спененный полиуретан.

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

где t _из – температура изоляции,

t _ос – температура окружающей среды, °С

t_из =50 °С

t_ос =20 °С

Тогда по таблице 9 [1] определяем

-число Прандтля

-коэффициент кинематической вязкости воздуха

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

Определяем число Нуссельта:

где Gr – число Гросгофа,

Pr – число Прандтля

где β – коэффициент объемного расширения воздуха, 1/К

Задаемся толщиной изоляции:

– наружный диаметр изоляции

ν – коэффициент кинематической вязкости воздуха, м/с²

Коэффициент объемного расширения воздуха находится из соотношения:

,

где t _пс – средняя температура пограничного слоя, °С

Число Нуссельта также определяется по формуле:

,

Для расчета изоляции воспользуемся равенством:

где t _т – температура наружной стенки обечайки, °С

Сосчитаем левую часть равенства:

Сосчитаем правую часть равенства:

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

Список литературы

1. Краснощеков Е.А. и Сукомел А.С. Задачник по теплопередаче. Изд. 2-е перераб. и доп. М. Энергия, 1969. – 264 с.

2. Калинин Э.К., Дрейцер Г.А., Копп И.З., Мякочин А.С. Эффективные поверхности теплообмена. – М.: Энергоиздат, 1998. – 408 с.: ил.

3. Кунтыш В.Б., Рощин С.П., Самородов А.В. Сравнение поверхностей теплообмена по энергетическим характеристикам: Методические указания к курсовому и дипломному проектированию. – Архангельск: Изд-во АГТУ, 2001. – 26 с.

4. Лебедев П.Д. и Щукин А.А. Теплоиспользующие установки промышленных предприятий. Учеб. Пособие для энергетических вузов и факультетов. – Москва: Изд-во «Энергия», 1970. – 408 с.

5. Водяные тепловые сети: Справочное пособие по проектированию/ Беляйкина И.В., Витальев В.П. и др.; Под ред. Н.И. Громова, Е.П. Шубина. – М.: Энергоатомиздат, 1988. - 376 с.

6. Лащинский А.А. Конструирование сварных химических аппаратов/ справочник. Л.: Машгиз, 1963. – 468 с.

7. Лащинский А.А. Конструирование сварных химических аппаратов / справочник. – Л.: Машиностроение, 1981. - 382 с.

8. Бакластов А.М. Проектирование, монтаж и эксплуатация теплоиспользующих установок. Учеб. Пособие для студентов специальности «Промышленная энергетика» высших учебных заведений. М., «Энергия», 1970.