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

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

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

 

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

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

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

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

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

В_с = В_б + С_п + В_ох,                                                                                                        (1)

где В_с — количество влаги в соломе, что попадает на прессование, кг

В_б — количество влаги в конечном продукте (брикете), кг;

С_п — потери влаги от нагрева в процессе прессования, кг

В_ох — потери влаги при охлаждении брикетов, кг.

Левая часть уравнения представляет количество влаги в сырье, попадает на прессование, а правая — распределение влаги после прессования (в брикете). Из общего баланса влаги количество, которое будет в образованном брикете

В_б = В_с - В_п - В_ох.                                                                                                        (2)

Общий баланс теплоты при брикетировании соломы записывается в виде

Q_т - Q_п = Q_б + Q_ох,                                                                                                       (3)

где Q_т — количество теплоты, выделяющейся при прессовании, Дж;

Q_п — потери теплоты в процессе прессования, Дж;

Q_б — количество теплоты, выносится образованными брикетами, Дж;

Q_ох — количество тепла, теряется при перемещении брикетов из пресса на сушильную установку, Дж.

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

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

На основе анализа литературных источников по тепломассообмену [1, 3], сушке [2, 4] и экспериментальных исследованиях были сделаны следующие расчеты и обобщения. Основной принцип предлагаемой схемы процесса сушки может быть представлен в следующем виде (рис. 1).

Рис. 1. Схема процесса сушки топливных брикетов

 

Принимаем, что входящий воздух соответствует характеристикам атмосферного воздуха: средняя температура t_o, ⁰С, относительная влажность φ_о, %.

Влагосодержание входящего атмосферного воздуха, кг/кг

Х₀ = 0,622· φ₀ · Р_нвп/(Р_атм — φ₀·Р_нвп),                                                                             (4)

где Р_нвп — давление насыщенного водяного пара при t_o, Па;

Р_атм — атмосферное давление, Па.

Удельная энтальпия входящего воздуха при принятых условиях, кДж/кг

h_взд0 = (С_взд + x · С_пв) · t_o+r₀ · Х,                                                                                   (5)

где С_взд — теплоемкость воздуха, кДж/кг×К;

С_пв — теплоемкость паров воды, кДж/кг×К;

r₀ — удельная теплота парообразования при t_o, кДж/кг.

Температура брикетов, которые выходят из пресса t_н, ^оС.

Определяем энтальпию сырого брикета (Н_брк0), равной количеству теплоты, что выносится образованными брикетами (Q_б) при t_н, ⁰С, кДж

Н_брк0 = Н_сс0 + Н_водж0 + Н_взд0,                                                                                         (6)

где Н_сс0, Н_водж0, Н_взд0 — соответственно энтальпия сухой соломы, свободной воды, влажного воздуха в сыром брикете, кДж.

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

Энтальпия сухой соломы Н_сс0 в брикете, кДж

Н_сс0= m_cc · C_cc · t_н,                                                                                                          (7)

где m_cc — масса сухой соломы в брикете, кг;

С_сс — принята теплоемкость сухой соломы, кДж/кг·К.

Энтальпия свободной воды Н_водж0 в брикете, кДж

Н_водж0 = m_вод · С_вж · t_н,                                                                                                   (8)

где m_вод — масса воды свободной в брикете, кг

С_вж — теплоемкость воды в брикете, кДж/кг ·К.

Энтальпия влажного воздуха в исходном брикете,

Н_взд0 = m_{. взд} · h_взд0,                                                                                                        (9)

где m_взд — масса сухого воздуха в порах брикета, кг

h_взд0 — удельная энтальпия влажного воздуха в брикете, кДж/кг.

Удельная энтальпия при принятых условиях влажного воздуха, находящегося в порах брикета при температуре брикета, кДж/кг

h_взд0 =(С_взд + Х_вздбр ·С_пв)·t_н + r_н ·Х_вздбр,                                                                         (10)

где С_взд — теплоемкость воздуха, кДж/кг ·К;

С_ПВ — теплоемкость паров воды, кДж/кг ·К;

Х_вздбр — влагосодержание при условиях формирования брикета, кг/кг

Х_вздбр = 0,622·φ_вздбр·Р_нвпбр /(Р_атм — φ_вздбр ·Р_нвпбр),                                                          (11)

где Р_нвпбр — давление насыщенного водяного пара в воздухе пор брикета при принятых условиях формирования брикета, Па;

r_н — удельная теплота парообразования при соответствующей температуре брикета, кДж/кг;

φ_вздбр — относительная влажность воздуха в порах сырого брикета, дол. ед.

На входе в 1 ступень аппарата.

Атмосферный воздух перед подачей в аппарат охлаждает брикеты, которые кулачковым транспортером направляются после сушки и охлаждения в накопитель (тару), подогревается до температуры t_п, ⁰С, но сохраняет влагосодержание атмосферного воздуха Х₀.

Удельная энтальпия входного подогретого атмосферного влажного воздуха, кДж/кг

h_вздп =(С_взд + x·С_пв) · t_п + r_п · X,                                                                                    (12)

где С_взд — теплоемкость воздуха, кДж/кг ·К;

С_ПВ — теплоемкость паров воды, кДж/г ·К;

r_п — удельная теплота парообразования, кДж/кг;

Х_п = Х₀ — влагосодержание атмосферного воздуха, кг/кг.

С формулы 4 для данных условий находим значение относительной влажности подогретого воздуха φ_п.

В высушенном брикете после 2 ступени аппарата содержится свободной воды, кг

m_вод2 = u_к × m_cс,                                                                                                              (13)

где U_к — влагосодержание, или абсолютная влажность высушенного брикета;

m_cc — масса сухой соломы в сыром брикете, кг.

Необходимо при сушке удалить из брикета воду, кг

ΔМ_вод = m_вод — m_вод2,                                                                                                   (14)

где m_вод — масса свободной воды в сыром брикете, кг.

Энтальпия сухой соломы Н_СС в брикете с конечной влажностью, кДж

Н_сс0= m_cc × C_cc × t_к2,                                                                                                         (15)

где m_cc — масса сухой соломы в брикете, кг

С_сс — теплоемкость сухой соломы, кДж / кг·К;

t_н2 — уточнена конечная температура готового брикета после 2 степени аппарата, ⁰С.

Энтальпия свободной воды Н_водж2 в высушенном брикете, кДж

Н_водк2 = М_водк2 × С_вж × t_к2,                                                                                               (16)

где М_водк2 — масса свободной воды в высушенном брикете, кг

С_ВЖ — теплоемкость жидкой воды, кДж/кг·К.

Удельная энтальпия влажного воздуха, находящегося в порах высушенного брикета, кДж/кг

h_вздк2 =(С_взд + Х_вздк2 × С_пв) × t_к2+r_к2 × Х_вздк2,                                                                    (17)

где С_взд — теплоемкость воздуха, кДж / кг·К;

С_ПВ — теплоемкость паров воды, кДж/кг·К;

Х_вздк2 — влагосодержание воздуха в порах высушенного брикета, кг/кг

Х_вздк2 =0,622× φ_вздк2×Р_вздк2/(Р_атм — φ_вздк2 × Р_вздк2),                                                          (18)

где Р_вздк2 — давление насыщенного водяного пара в воздухе пор высушенного брикета, Па;

r_к2 — удельная теплота парообразования при конечной температуре готового брикета после 2 ступени аппарата, кДж/кг;

φ_вздк2 — относительная влажность воздуха в порах высушенного брикета после второй ступени аппарата, дол. ед.

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

Н_вздк2 = m_взд × h_вздк2,                                                                                                       (19)

где m_взд — масса сухого воздуха в порах высушенного брикета, кг

h_взд0 — удельная энтальпия влажного воздуха при конечной температуре готового брикета после 2 степени аппарата, кДж/кг.

Полная энтальпия высушенного брикета, кДж

Н_брк2= Н_сск2 + Н_водк2 + Н_вздк2.                                                                                        (20)

Изменение энтальпии брикета в теплообменном аппарате, кДж

ΔН_бр= Н_брк0 — Н_брк2.                                                                                                     (21)

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

-          на поверхности брикета идет нагрев потока воздуха конвективным путем с отбором части тепла от брикета;

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

Расчет конвективного теплообмена между поверхностью брикета и потоком воздуха осуществляли, используя приемы теории подобия [2]. Учитывая малое число рядов брикетов на транспортере, для определения среднего по поверхности брикетов коэффициента теплоотдачи, используем критериальное уравнение Жукаускаса, предлагаемое для расчета теплообмена при омовения потоком теплоносителя одиночного цилиндра

Nu_жd= 0,25Re_жd ^0,6 × Pr_ж ^0,38 (Pr_ж / Pr_с) ^0,25                                                                      (22)

где Nu_жd — число Нуссельта,

Re_жd — число Рейнольдса,

Pr_ж — число Прандтля.

Учитывая малую зависимость Pr для газов от температуры, формула может записываться

Nu_жd= 0,25Re_жd ^0,6.

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

Выполнен поливариантный расчет условий теплообмена в первом ряду брикетов. При расчетах варьировали скоростью потока воздуха, принимая ее дискретные значения 2; 1; 0,5 и 0,25 м/с. Указанные скорости, как показали расчеты, обеспечивали турбулентный режим течения воздуха и пригодность уравнения Жукаускаса.

Расчет числа Рейнольдса выполняется по формуле

Re_жd=w × d_бр / ν_ж,                                                                                                           (23)

где w — скорость линейная потока воздуха при температуре, м/с;

ν_ж — кинематическая вязкость воздуха при температуре, м²/с.

Значение чисел Рейнольдса для принятого интервала скоростей удовлетворяет требованию по использованию уравнения Жукаускаса

1 × 10³ < Re_жd < 2 × 10⁵.                                                                                                 (24)

Расчет размерного коэффициента теплоотдачи α вели исходя из формулы

Nu_жd = α × d_бр / λ,                                                                                                           (25)

λ — коэффициент теплопроводности;

при скорости потока воздуха w, м/с коэффициент теплоотдачи α находим как

α = Nu_жd × λ/d_бр.                                                                                                                                                 (26)

Расчет первой ступени аппарата

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

Поверхность теплоотдачи одиночного брикета F, м² определяется как

F = π × d_бр × L_бр + π × d_бр²/2,                                                                                           (27)

где L_бр — длина брикета, м.

Сушка ведется на подвижном кулачковом транспортере.

Уравнения передачи тепла Q, кДж одиночным брикетом

Q = F × α × Δt_ср × τ,                                                                                                          (28)

где τ — экспериментальный время пребывания брикета в зоне первичной сушки, с;

Δt_ср — средняя разность температур между потоком воздуха и стенкой брикета,

Δt_ср = [(t_с1 — t₀) + (t_с1 — t_k1)]/2,                                                                                    (29)

где t_с1 — средняя температура поверхности брикетов на 1 степени аппарата, ⁰С;

t₀ — температура воздуха на входе в первую ступень аппарата, ⁰С;

t_k1 — температура воздуха, которая подлежит определению на выходе с 1 ступени аппарата, ⁰С.

Уравнение расчета количества тепла полученного потоком воздуха за время τ от одного брикета на первой ступени аппарата имеет вид, кДж

Q = M_взд × C_взд (t_k1 — t₀).                                                                                               (30)

Масса воздуха, который подогревается омывая один брикет на первой ступени аппарата за время τ, кг

М_взд = s × w × ρ_взд × τ,                                                                                                       (31)

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

ρ_взд — плотность воздуха подогретого за счет охлаждения брикетов ⁰С, кг/м³;

w — скорость потока воздуха в узком месте между брикетами, м/с.

Можно записать

Q = Q_пвзд,                                                                                                                       (32)

где Q_пвзд — количество тепла, затрачиваемое на подогрев конвекцией внешнего потока воздуха за время прохождения брикета через первую ступень аппарата, кДж;

Q_пвзд = ΔН_бр — Н_пвод,                                                                                                    (33)

где ΔН_бр — общее количество тепла, отданное брикетом на испарение внутренней воды и на подогрев конвекцией внешнего воздушного потока за время пребывания брикета на первой ступени аппарата, кДж.

ΔН_бр = Н_брк0 — Н_брк1.                                                                                                                                     (34)

Находим энтальпию подсушенного брикета на выходе из первой ступени аппарата, кДж

Н_брк1 = [m_cc × С_сс + m_вод1 × С_вж + m_взд1 × С_взд] × t_брk1,                                                       (35)

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

Энтальпия паров воды диффундирует из брикета в воздушный поток, кДж

Н_пвод = h_tс × d_mв1,                                                                                                            (36)

где h_tс — удельная энтальпия паров воды диффундирует из подсушенного брикета при температуре стенки брикета после первой степени аппарата, кДж/кг;

d_mв1 — масса воды, удаляемой из брикета на первой ступени аппарата за счет тепловой энергии, которая вносится горячим брикетов, кг.

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

Тепловая мощность вентилятора Q_эд, кВт, передана потока воздуха определяется

Q_эд = 0,8 × W_эд,                                                                                                              (37)

где W_эд — установочная мощность электродвигателя вентилятора, кВт;

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

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

Тепловая мощность дополнительного источника тепла Q_пв1, кДж/с, что действует в брикете за счет изменения параметров потока воздуха и которая влияет на процесс испарения,

Q_пв1 = 0,7 × Q_ед.                                                                                                              (38)

Дополнительное количество энергии q_пвбр, кДж, подводимой каждому брикету за счет изменения параметров потока воздуха

q_пвбр = Q_пв1 × τ/n,                                                                                                            (39)

где n — количество брикетов, которые высушиваются за 1 час.

Оценка количества влаги d_mпв1, которая удаляется из брикета за счет действия дополнительного источника теплоты проводится по соотношению

d_mпв1 = q_пвбр / h_tс,                                                                                                            (40)

Общее количество влаги удаляется из брикета на 1 ступени аппарата, кг

d_mобщ1 = d_mв1 + d_mпв1.                                                                                                     (41)

Повышение влагосодержания воздуха Δх₁ на выходе из первой ступени аппарата, кг/кг сухого воздуха

Δх₁ = d_mобщ1/М_взд.                                                                                                          (42)

Влагосодержание воздуха Х_вздк1 на выходе из первой ступени аппарата, кг/кг

Х_вздк1 = Х + Δх₁.                                                                                                            (43)

Удельная энтальпия влажного воздуха на выходе из 1 ступени аппарата, кДж/кг

h_вздк1 = (С_взд + Х_вздк1 × С_пв) × t_к1 + r_к1 × Х_вздк1,                                                                 (44)

где R_к1 — удельная теплота парообразования при конечной температуре потока, кДж/кг.

Расчет второй ступени аппарата

Среднее содержание воды в брикетах m_вод1 на входе во вторую ступень аппарата, кг

m_вод1 = m_вод — d_mобщ1.                                                                                                   (45)

Количество воды Δ_mвод2, которую необходимо удалить из каждого брикета на второй ступени аппарата

Δm_вод2 = m_вод1 — m_водк2.                                                                                                (46)

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

Δх₂ = Δm_вод2 / М_взд.                                                                                                        (47)

Конечное влагосодержание воздуха на выходе из второй ступени аппарата, кг/кг

Х_вздк2 = Х_вздк1 + Δх₂.                                                                                                      (48)

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

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

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

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

Удельная энтальпия влажного воздуха на выходе из второй ступени аппарата, кДж/кг

h_вздк2 = (С_взд + Х_вздк2 × С_пв) × t_к2 + r_К2 × Х_вздк2,                                                                (49)

где t_к2 — температура отработанного воздуха на выходе из 2 ступени аппарата, ⁰С;

r_К2 — удельная теплота парообразования при соответствующей температуре (⁰С) брикетов на выходе из второй ступени, кДж/кг.

Энтальпия высушенного брикета при соответствующей температуре (⁰С) брикетов на выходе из 2 степени аппарата, кДж

Н_брк2 = [m_cc × С_сс + m_вод2 × С_вж + m_взд2 × Св_зд] × t_брk2.                                                      (50)

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

Приход тепла

Q_прих = Q_вздк1 + Q_брк1 + Q_пвм2 + Q_подог,                                                                           (51)

где Q_вздк1 — теплота, поступающая с воздухом, омывает брикет на выходе с 1 ступени аппарата, кДж

Q_вздк1 = М_взд × h_вздк1,                                                                                                       (52)

где Q_брк1 — теплота, поступающая из брикетов с 1 ступени при соответствующей его температуре на выходе с 1 ступени (⁰С), кДж

Q_брк1 = Н_брк1,                                                                                                                  (53)

где Q_пвм2 — теплота превращения механической энергии потока воздуха в тепловую энергию на 2 ступени аппарата

Q_пвм2 = Q_пв2 × τ/120,                                                                                                       (54)

где Q_подог — теплота, которую необходимо подвести от внешних подогревателей воздуха перед подачей его на 2 ступень аппарата для реализации процесса досушки одиночного брикета, кДж;

τ — экспериментальное время пребывания брикета в зоне вторичной сушки (равное времени пребывания в зоне первичной сушки), с;

Расходы тепла

Q_расх = Q_вздк2 + Q_брк2,                                                                                                     (55)

где Q_вздк2 — тепло, выводится из 2 степени аппарата с отработанным воздухом, кДж

Q_вздк2= М_взд × h_вздк2,                                                                                                        (56)

где Q_брк2 — тепло выведено сухими брикетами из 2 ступени аппарата при соответствующей температуре брикетов, ^оС

Q_брк2 = Н_брк2.                                                                                                                  (57)

Удельная энтальпия подогретого воздуха на входе во 2 ступень аппарата

h_вздн2 = Q_вздк2 /М_взд.                                                                                                        (58)

Расчет удельной энтальпии воздуха после подогревателя, кДж/кг

h_кал = (С_взд + Х_вздк1 × С_пв) × t_кал + r_н2 × Х_кал,                                                                     (59)

где t_кал — температура подогретого воздуха на входе во вторую ступень аппарата, ⁰С;

r_н2 — удельная теплота парообразования при t_кал, кДж/кг;

Х_кал = Х_вздк1 — влагосодержание воздуха Х_вздк1 на выходе из первой ступени аппарата.

Определение общей мощности дополнительного нагревателя, кВт

W_под = Q_подог × n/3600,                                                                                                   (60)

где n — число брикетов, что высушиваются за 1час.

Выводы

В результате расчетов при максимально возможной входной влажности брикетов 20 % для достижения исходной влажности 14 %, которая обеспечивает эффективность сжигания и долговременное хранение брикетов, было определено, что на выходе из первой ступени сушки и охлаждения температура брикетов будет составлять 60 °С, влажность 16 %. Для обеспечения досушивания брикетов была установлена необходимость подогрева воздуха на входе во вторую ступень до 73 ⁰С. Установленная мощность подогревателя (электрокалорифера) 2,6 кВт, при фактической производительности пресса 90 кг/ч. За счет этого обеспечивается досушивание брикетов до необходимой влажности — 14 %, при этом их температура на выходе из второй ступени составит 55 °С, а температура воздуха, который оставляет аппарат будет составлять 58 °С.

 

Литература:

 

1.         Кутателадзе С. С. Основы теории теплообмена / С. С. Кутателадзе. — 5-е изд., перераб. и доп. — М.: Атомиздат, 1979. — 416 с.

2.         Лыков М. В. Сушка в химической промышленности / М. В. Лыков. - М.: Химия, 1970. - 432 с.

3.         Лыков А. В. Тепломассообмен: справочник / А. В. Лыков. — 2-е изд., перераб. и доп. - М.: Энергия, 1978. - 480 с.

4.         Сажин Б. С. Основы техники сушки / Б. С. Сажин. — М.: Химия, 1984. — 320 с.