# Макрокинетика сушки

Герман Иванович Ефремов

Изложены теоретические основы и экспериментальные данные по кинетике сушки химических, пищевых и текстильных материалов. Показана аналогия различных процессов переноса массы, тепла и импульса на примерах различного вида обработки ряда химических, пищевых и волокнистых материалов. Рассмотрены типы и свойства тепло- и влагоносителей и свойства высушиваемых материалов. Показано развитие теории статики и кинетики сушки и увлажнения (внешняя и внутренняя задачи) на основе аналитических решений уравнения диффузии и квазистационарного метода расчета коэффициентов диффузии в процессах сушки. Описаны конструкции и работа ряда сушильных установок. Изложены инженерные методы расчета макрокинетики процессов сушки. Монография предназначена для научных, инженерно-технических работников и аспирантов, занимающихся рассматриваемой проблемой. Она может быть использована также в качестве учебного пособия для студентов химико-технологических, текстильных, пищевых и других ВУЗов и колледжей.

## Предисловие

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

Макрокинетика процессов переноса массы, тепла и импульса, что характерно и для процессов сушки, как самостоятельный раздел науки впервые была выделена профессором Д. А. Франком-Каменецким. Им, по словам академика Н. Н. Семенова, показано, что вопросы, рассматриваемые в этом разделе науки, имеют не только частное и прикладное, но и общенаучное значение и представляют общий теоретический научный интерес.

Общеизвестна фундаментальная монография Д. А. Франка-Каменецкого “Диффузия и теплопередача в химической кинетике”, выдержала три издания на русском языке и переиздана за рубежом. В ней автор рассматривает макрокинетику химических реакций и теорию горения в сочетании с диффузией, теплопередачей и гидродинамикой. Такое сочетание позволило решить ряд прикладных нелинейных задач, разработать новые методы расчета кинетики химических реакций и ряда нестационарных технологических процессов.

В развитие исследований макрокинетики процессов переноса внесли свой вклад такие известные ученые как академик Н. М. Жаворонков (теория гидродинамики и процессов химической технологии), профессор О. Левеншпиль (расчет и инженерное оформление химических реакторов), академик В.В. Кафаров (математическое моделирование химико-технологических процессов), профессор Я. Циборовский (тепло массообменные процессы) и многие другие исследователи.

Вопросы макрокинетики процессов рассматриваются также в настоящее время и под другими названиями “Теория диссипативных структур” – изучение диссипативных процессов, “Синергетика”– наука, изучающая периодические химические реакции и пространственную их картину, “Теория катастроф” – описание быстропротекающих процессов от механики полимеров до биологических процессов. “Теория катастроф” берет свое начало от теории теплового взрыва, разработанной академиком Н. Н. Семеновым, профессорами Д. А. Франком-Каменецким и О.М. Тодесом.

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

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

В книге приведены экспериментальные данные, которые были получены на кафедрах процессов и аппаратов Московского государственного открытого университета им. В. С. Черномырдина и Московского государственного текстильного университета им. А.Н. Косыгина, а также и литературные данные по сушке различных материалов.

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

## Основные обозначения

### Латинские


### Греческие


## Индексы


## Введение

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

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

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

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

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

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

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

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

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

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

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

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

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

Следует отметить, как важную особенность книги, что ряд расчетов и примеров в данной монографии с необходимыми пояснениями и программами выполнен в компьютерных средах MathCAD, MathLab и FemLab.

## Глава I

### Общие закономерности процессов переноса


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

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

Механические – процессы механической обработки твердых материалов (основываются на законах механики твердого тела). Это процессы резания, деформации, дробления, смешения, сепарации и т. д.

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

Тепловые – процессы разного вида распространения тепла в гомогенных и гетерогенных системах, основаны на законах теплопроводности, теплообмена (перенос тепла). Эти процессы нагрева, охлаждения, конденсации, испарения т. д.

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

Химические – процессы получения новых веществ на основе протекающих химических реакций, описываются законами химической кинетики. Эти процессы получения аммиака, каучука, различных кислот, щелочей, солей, горения топлив т. д.

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

#### 1.1 Закон сохранения массы

Одним из главных законов при переносе массы является закон сохранения массы. Этот закон установлен М. В. Ломоносовым. Для элементарного объема он может быть получен следующим образом.


Рис. 1.1 К выводу закона сохранения массы.


Рассмотрим поток вещества через грани элементарного объема. Плотность ρ и скорость потока u в общем случае изменяются в пространстве и во времени:


Рассмотрим изменение массы вдоль оси х (Рис. 1.1). Если проекция скорости потока на входе в элементарный объем ux, то на выходе из него, с учетом изменения на длине dx она составит:


.

Тогда изменение массы вдоль оси х за счет изменения скорости составит:


.

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


Выражение в скобках в уравнении (1.2) называется дивергенцией вектора скорости и обозначается div u. С учетом его получим для (1.2):


Это выражение закона сохранения массы и оно известно в гидродинамике, как уравнение сплошности, неразрывности потока. В элементарной форме это уравнение для одномерного потока, движущегося со средней скоростью v примет вид:


где М – массовый расход потока, S – площадь его поперечного сечения.

Для несжимаемых жидкостей (ρ = Const) уравнение (1.3) упрощается:


Для описания химического процесса в уравнении (1.2) вместо плотности подставляют массовую концентрацию компонента С. С учетом скорости образования этого компонента по химической реакции r, если она имеет место, для уравнения (1.2) получим:


С учетом, что концентрация компонента изменяется в пространстве и во времени, получим:


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

#### 1.2 Закон сохранения количества движения

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


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


Рис. 1.2 К выводу закона сохранения количества движения.


Рассмотрим равновесие сил в движущемся потоке в проекциях на ось х (Рис. 1.2). На правую и левую грани действуют силы давления. Их проекция на ось х составит –


.

Проекция массовой силы Q на ось х запишется:


.

На верхнюю и нижнюю грани действуют силы вязкостного трения. Их проекция на ось х составит .

С учетом закона Ньютона для вязкостного трения:


имеем проекцию сил вязкостного трения на ось х:


.

Здесь выражение в скобках – оператор Лапласа от проекции скорости на ось х, он обозначается Δux или 2ux.

Так как сумма проекций всех сил равна проекции силы инерции:


,

относя все силы к единице объема, получим:


Последние два уравнения получены аналогично для осей у и z, а в целом система уравнений (1.10) в гидрогазодинамике называется уравнениями движения вязкой жидкости Навье-Стокса и выражает закон сохранения количества движения.

Система уравнений Навье-Стокса может быть записана более детально, если раскрыть полную производную проекции скорости. Так уравнение для оси х, например, при делении всех его членов на ρ , c учетом, что ν = μ/ρ, будет иметь вид:


.

Аналогично записываются уравнения для осей у и z.

#### 1.3 Закон сохранения энергии

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


.

Помножим эти уравнения соответственно на dx, dy, dz и сложим. Тогда, преобразуя, получим следующее уравнение:


В поле силы тяжести (Х = 0; У = 0; Z = – g) уравнение (1.12) примет вид:


Это уравнение определяет в дифференциальном виде закон сохранения энергии для движения идеальной жидкости и представляет собой, соответственно, сумму удельных (отнесенных к единице массы) потенциальных энергий положения и давления и кинетической энергии. При интегрировании уравнения (1.13) для потока несжимаемой жидкости (ρ = Const) получим уравнение Бернулли для одномерного потока, движущегося со средней скоростью v:


Уравнение Бернулли показывает, что для идеальной жидкости сумма потенциальной и кинетической энергий остается постоянной вдоль всего потока.

В более общей форме закон сохранения энергии описывает 1-й закон термодинамики: теплота, подводимая к системе, идет на производство работы и увеличение энергии системы:


Выражение для потока энергии в дифференциальном виде включает ее члены, входящие в уравнение (1.13) плюс, внутренняя энергия dU. С учетом этого запишем уравнение (1.15) в следующем виде:


Сумма второго и третьего членов правой части уравнения (1.16) представляет собой изменение энтальпии dh. С учетом этого получим другой вид уравнения (1.16):


#### 1.4 Микро- и макроперенос

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

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

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

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

Макроперенос – это перенос определенных объемов массы, перенос энергии этих объемов, перенос количества движения, которым обладают эти объемы. Макроперенос обусловлен наличием конвекции (свободной или вынужденной), вихреобразованиями. Кинетика процесса макропереноса называется макрокинетикой.

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

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

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

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

#### 1.5 Перенос массы

Рассмотрим поток вещества J в пространстве (Рис. 1.3). Концентрация вещества С изменяется в пространстве и во времени т. е. .

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


Рис. 1.3 Поток вещества в пространстве.


Поток вещества J, отнесенный к единице поверхности S в этом случае будет пропорционален градиенту концентрации:


Это выражение называется 1-м законом Фика. Здесь D – коэффициент диффузии.

Закон молекулярного переноса массы (вещества) для элементарного объема может быть получен следующим образом. Рассмотрим поток вещества через грани элементарного объема.

Расход вещества, обусловленный диффузией вдоль оси х через левую поверхность dy·dz (Рис. 1.4), с учетом закона Фика, составит.


На расстоянии dx, с учетом изменения концентрации, получим ее значение


.

Тогда расход вещества, обусловленный диффузией вдоль оси х через правую грань dy·dz, составит

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


Рис. 1.4 К выводу 2-го закона Фика.


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


Выражение в скобках в уравнении (1.19) является оператором Лапласа. С учетом его сокращенного обозначения и образования вещества за счет химической реакции при скорости образования r получим для выражения (1.19):


Уравнение (1.20) является следствием закона сохранения массы и 1-го закона Фика и называется 2-м законом Фика. Оно определяет поле концентраций для молекулярной диффузии в рассматриваемой среде.

Для изотропной диффузии уравнение (1.20) может быть записано для изменения вдоль любой оси, например х, в следующем виде:


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


Уравнение (1.22) определяет поле концентраций с учетом молекулярной и конвективной диффузии. Как следует из этого уравнения, для описания макрокинетики процесса переноса массы необходимо учитывать гидродинамику и поэтому уравнение переноса в движущейся среде (1.22) должно быть дополнено уравнениями движения Навье-Стокса (1.10), определяющими поле скоростей и перенос импульса в движущейся среде.

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

#### 1.6 Перенос тепла

Рассмотрим молекулярный перенос тепла аналогично переносу массы (Рис. 1.3). Температура t в общем случае изменяется как в пространстве, так и во времени т. е. .

Тогда поток тепла q, отнесенный к единице поверхности S будет пропорционален градиенту температуры в произвольном направлении n:


Это выражение называется законом теплопроводности Фурье. Здесь – коэффициент теплопроводности среды.

Закон молекулярного переноса теплоты для элементарного объема может быть получен аналогичным образом, как и для переноса массы. Тогда суммарное изменение количества тепла вдоль всех 3-х координат элементарного объема, проявляющееся в изменении температуры, равно изменению температуры, вовремя помноженному на произведение теплоемкости ср на плотность ρ [6]:


Выражение в скобках в уравнении (1.24) является оператором Лапласа, с учетом его сокращенного обозначения и дополнительного источника тепла qr, например, за счет химической реакции, с введением коэффициента температуропроводности а = / ср ρ, получим для выражения (1.24):


Уравнение (1.25) является следствием закона Фурье и также называется его именем. Оно определяет поле температур для молекулярного переноса тепла в рассматриваемой среде.

Для одномерной теплопроводности, например, вдоль оси х, уравнение (1.25) может быть записано в следующем виде:


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


Уравнение (1.27) называется уравнением Фурье-Кирхгофа и определяет поле температур с учетом молекулярного и конвективного переноса тепла. Как следует из этого уравнения, для описания макрокинетики процесса теплопереноса необходимо учитывать гидродинамику и уравнение переноса в движущейся среде (1.27) должно быть дополнено уравнениями движения Навье-Стокса (1.10), определяющими поле скоростей в движущейся среде.

#### 1.7 Перенос количества движения

По второму закону Ньютона изменение количества движения в единицу времени (импульс) численно равно силе – уравнение (1.8). В движущемся потоке газа или жидкости под действием массовых и поверхностных сил происходят соударения молекул, что обуславливает перенос количества движения. Баланс сил в движущемся потоке представляет собой закон сохранения количества движения (импульса). На основе баланса сил получена выше система уравнений Навье-Стокса (1.10).

Строго говоря, т. к. система уравнений Навье-Стокса получена на основе закона Ньютона для вязкостного трения (1.9), учитывающего молекулярный перенос количества движения (микроперенос), она применима только для струйчатого, ламинарного движения вязкой жидкости.

Если жидкости отклоняются от закона Ньютона, то их называют неньютоновскими. К ним относятся жидкие полимеры, растворы высокомолекулярных полимеров, суспензии и др.

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


Можно показать [6], что в этом случае система уравнений Навье-Стокса формально останется той же (1.10), но вместо ньютоновской вязкости μ необходимо подставить в нее сумму ньютоновской и турбулентной вязкости μт.

Поле скоростей при микропереносе количества движения определяется при решении системы уравнений Навье-Стокса (1.10), а при макропереносе количества движения определяется из решения системы уравнений Навье-Стокса с учетом турбулентной вязкости. Следует отметить, что если ньютоновская вязкость является величиной постоянной для данной среды, то турбулентная вязкость зависит от масштаба турбулентности и поэтому при ее расчете, возникают определенные трудности.

#### 1.8 Аналогия процессов переноса

Сравнивая выражения для трех видов молекулярного переноса массы (1.18) – закон Фика, переноса тепла (1.23) – закон Фурье и переноса количества движения (1.9) – закон Ньютона, нетрудно заметить, что по форме они абсолютно аналогичны.


.

Здесь эти уравнения дополнены еще одним видом переноса – переносом электричества. Где u – градиент напряжения или электрический потенциал, ρ – удельное сопротивление электрического тока.

Имеется очевидное подобие выражений для конвективного переноса массы (1.22), переноса тепла (1.27) – закон Фурье и переноса количества движения (1.10), записанных ниже для одномерного потока:


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

На аналогию трех видов переноса указывает и одинаковая размерность коэффициентов переноса а, D и ν равная [м2/с]. Для газов механизм молекулярного переноса одинаков для переноса массы, теплоты и количества движения, т. к. связан с тепловым движением молекул. Поэтому все три коэффициента имеют один порядок с произведением длины свободного пробега молекул L на скорость их теплового движения u*:


.

Для жидкостей соотношение этих коэффициентов следующее:

ν а D.

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

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

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

#### 1.9 Экспериментальный подход к исследованию процессов переноса

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

Так для расчета процессов теплопереноса используют эмпирический коэффициент теплоотдачи , равный отношению потока тепла q к разности температур Δt. Тогда тепловой поток равен:


Распределение температур от одной среды к другой в стационарной теплопередаче (постоянство во времени потока тепла q) через стенку толщиной δ с коэффициентом теплопроводности приведено на Рис. 1.5.

Применяя уравнение (1.29) для 1-й и 2-й среды, с учетом теплопроводности через стенку толщиной δ и общего коэффициента теплопередачи К, получим равенство выражений для стационарного теплового потока:


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


С учетом этих соотношений получим уравнение:


Физический смысл соотношения (1.31) заключается в том, что общее сопротивление теплопередачи через стенку 1/К равно сумме сопротивления переноса тепла от 1-й среды к стенке 1/α1, термосопротивления стенки δ/ и сопротивления переноса тепла от стенки ко 2-й среде 1/α2.

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

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

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


Значения концентраций на границе xs и уs трудноопределимы, поэтому записывают другое выражение диффузионного потока для коэффициентов массопереноса для первой среды К1 и для второй К2 через соответственно равновесные концентрации x* и y*.


Рис. 1.6. Распределение концентраций при массопередаче через поверхность раздела фаз.


Обычно принимают линейный закон (m – константа равновесия, тангенс угла наклона линии равновесия) для определения равновесных концентраций на границе [6]:


Из очевидного равенства:


находим с учетом (1.32) и (1.33):


К1 – коэффициент массопереноса по газовой фазе.

Физический смысл соотношения (1.36) заключается в том, что общее сопротивление процесса массопередачи через межфазную поверхность 1/К1 равно сумме сопротивления переноса вещества от 1-й среды к межфазной поверхности 1/1 и сопротивления переноса вещества от межфазной поверхности ко 2-й среде m/2.

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


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

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

#### 1.10 Основы теории подобия

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

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

Соотношения между сходственными величинами образца и модели называются константами (масштабами) подобия. Так для геометрического подобия – подобия геометрических размеров образца (l’ – длина, b’– ширина, h’– высота) и соответствующих размеров модели – l”, b”, h” получим константу геометрического подобия


Из этого выражения можно записать инварианты (симплексы) подобия, записав отношения параметров для образца и для модели, например:


Для подобия физических величин имеем константы подобия, например для кинематической вязкости и плотности:


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


Константы кинетического подобия включают отношение скоростей u и ускорений a в сходственных точках объекта и модели:


Следует отметить, что подобие кинетическое (подобие планов скоростей и ускорений) может иметь место только при наличии подобия геометрического.

Динамическое подобие описывается соотношениями сходственных сил, энергий и масс для модели и объекта, причем динамическое подобие может иметь место только при наличии подобия геометрического и кинетического.

#### 1.11 Подобное преобразование дифференциальных уравнений переноса

Перенос количества движения. Рассмотрим динамическое подобие на примере подобного преобразования уравнения Навье-Стокса (1.10) для одномерного потока (перенос количества движения). В этом случае скорость потока .

Запишем это уравнение для объекта:


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


Чтобы уравнения (1.42) и (1.43) были идентичны, все множители из констант подобия должны быть равны, тогда:


Домножив все члены уравнения (1.44) на kl /kv2, получим безразмерные соотношения:


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


Опуская индексы, получим аналогично и другие критерии:


.

Для обеспечения подобия критерии подобия для объекта и модели должны быть численно одинаковы. Это положение определяет содержание 1-й теоремы подобия.

Так как члены уравнения Навье-Стокса представляют собой сумму удельных сил тяжести, давления, вязкости и инерции, то сравнивая их с выражениями критериев, можно сделать вывод, что критерий Фруда представляет собой отношение силы инерции к силе тяжести, критерий Эйлера – отношение силы давления к силе инерции, критерий Рейнольдса – отношение силы инерции к силе вязкости (учитывает режим движения) и критерий гомохронности – отношение инерционных сил модели и объекта (учитывает нестационарность процесса).

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


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

Выше было сказано, что для обеспечения подобия критерии подобия для объекта и модели должны быть численно равны, но тогда согласно зависимости (1.46) и определяющие критерии должны быть равны. Таким образом, для обеспечения подобия достаточно равенства определяющих критериев модели и объекта. Это положение определяет содержание 3-й теоремы подобия.

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

Теплоперенос. Критерии теплопереноса выводятся аналогично выводу критериев переноса количества движения. Рассмотрим подобное преобразование уравнения Фурье-Кирхгофа (1.27) для одномерного теплопереноса, когда t = f (x, τ) в отсутствии источников тепла для объекта моделирования (индексы опустим):


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


Чтобы уравнения (1.47) и (1.48) были идентичны, все множители из констант подобия должны быть равны, тогда:


Из первого равенства выражения (1.49), подставив значения констант подобия, получим безразмерный комплекс – критерий Фурье, который характеризует изменение теплопереноса теплопроводностью во времени:


Из второго равенства выражения (1.49), подставив значения констант подобия, получим другой безразмерный комплекс – критерий Пекле, который представляет собой отношение теплопереноса за счет движения среды (конвективный) и за счет теплопроводности (молекулярный):


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

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

Так как в критерии Fr и Re входит скорость, используется их комбинация – критерий Галилея, не содержащий скорости:


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


Если разность плотностей вызвана термическим расширением среды (T – коэффициент термического расширения), то:


,

где Т – абсолютная температура.

После подстановки этого соотношения в критерий Архимеда получим критерий Грасгофа:


Для газов ,

тогда критерий Грасгофа примет вид:


При сочетании критерия Рет с критерием Re получим критерий Прандтля, характеризующий теплофизические свойства среды или соотношение полей скоростей и температур:


Если граничным условием теплопереноса является линейность теплового потока – уравнение (1.23) тогда в сочетании с уравнением (1.29) получим:


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


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


Для стационарного теплопереноса из зависимости (1.60) исключаются критерии. Но и Foт, содержащие время.

Критерий Nu является наиболее удобной величиной для расчета молекулярного переноса или переноса в ламинарном движении [10]. В случае развитой турбулентности более удобным оказывается использование критерия Стентона:


Тогда в зависимости (1.61) определяющим критерием вместо Nu будет критерий Stт.

Массоперенос. Критерии массопереноса и общая критериальная зависимость выводятся аналогично теплопереносу. Для этого рассмотрим подобное преобразование дифференциального уравнения массопереноса (1.22). Для одномерного потока, когда С = f (x, τ) в отсутствии источников для объекта моделирования запишем (индексы опущены):


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


Чтобы уравнения (1.62) и (1.63) были идентичны, все множители из констант подобия должны быть равны, тогда:


Из первого равенства выражения (1.64), подставив значения констант подобия, получим безразмерный комплекс – диффузионный критерий Фурье, который характеризует изменение молекулярного массопереноса во времени


Из второго равенства выражения (1.64), подставив значения констант подобия, получим другой безразмерный комплекс – диффузионный критерий Пекле, который представляет собой отношение массопереноса за счет движения среды (конвективный) и за счет молекулярной диффузии (молекулярный)


Он аналогичен по форме критерию Рейнольдса, и тепловому критерию Пекле.

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

При сочетании критерия Ре с критерием Re получим критерий Шмидта (в отечественной литературе его часто называют диффузионным критерием Прандтля), характеризующий соотношение молекулярного переноса количества движения и молекулярной диффузии или соотношение полей скоростей и концентраций:


Если граничным условием массопереноса является линейность потока вещества – уравнение (1.18), тогда в соответствии с уравнением (1.32) получим равенство:


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


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


Для стационарного массопереноса из зависимости (1.70) исключаются критерии. Но и Fo, содержащие время.

Критерий Sh является наиболее удобной величиной для расчета молекулярного переноса или переноса в ламинарном движении [6]. В случае развитой турбулентности более удобным оказывается использование диффузионного критерия Стентона, т. к. он в этом случае близок к постоянному значению


Тогда в зависимости (1.70) определяющим критерием вместо критерия Sh будет критерий St. Общий вид зависимости при этом сохранится.

#### 1.12 Внешняя и внутренняя задачи

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

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

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

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

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

Необходимо отметить, что поскольку условия переноса (граничные условия) во внешней и внутренней задачах различны, решения дифференциальных уравнений и значения критериев для этих задач также будут различными. Так, например, при движении шара в сплошной среде ламинарный режим сохраняется до значения критерия Рейнольдса равном 2, а при движении среды в трубах – до значения 2100 – 2300.

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

Примером совместной внешней и внутренней задач при теплопереносе является теплопередача через стенку (Рис. 1.5). Здесь внешняя задача – конвективный теплоперенос от среды к стенке, а внутренняя – теплопроводность внутри стенки. Если лимитирует теплопроводность (материал стенки теплоизолирующий), то коэффициент теплопередачи по уравнению (1.31). К , т. е. конвективными сопротивлениями можно пренебречь. Теплопередачу в условиях совместной внешней и внутренней задач характеризует тепловой критерий Био (аналог критерия Нуссельта)


Если в критерии Нуссельта оба параметра α и относятся к одной среде, то в критерии Био α – определяет конвективный теплообмен от среды к стенке (или наоборот) – внешняя задача, а параметр / l – определяет теплопроводность стенки – внутренняя задача. Если Biт 0, то лимитирует внешняя задача, если Biт → ∞, то лимитирует внутренняя задача. Если влияние внешней и внутренней задач при теплопередаче одного порядка, то говорят о смешанной задаче передачи тепла.

Примером совместной внешней и внутренней задач при массопереносе является процесс конвективной сушки пористого материала. Изменение влажности пористого материала происходит при его сушке потоком горячего воздуха. Здесь внешняя задача – конвективный массоперенос от среды к материалу, а внутренняя – перенос влаги (массопроводность) внутри материала. Перенос влаги внутри материала может быть учтен коэффициентом диффузии D внутри материала. Если лимитирует массоопроводность, то процесс массопередачи определяет диффузия внутри материала, т. е. конвективным массопереносом можно пренебречь. Массопередачу в условиях совместной внешней и внутренней задач характеризует диффузионный критерий Био (аналог критерия Шервуда)


Однако, если в критерии Шервуда оба параметра и D относятся к одной среде, то в критерии Био – определяет конвективный массооперенос от материала к среде (или наоборот при сорбции) – внешняя задача, а параметр D / l – определяет массопроводность материала – внутренняя задача. Если Bi → 0 (практически при Bi 0,2), то лимитирует внешняя задача, если Bi → ∞ (практически при Bi 50), то лимитирует внутренняя задача. Если значения критерия Био лежат в интервале от 0,2 до 50, то имеет место смешанная задача и влияние обеих внешней и внутренней задач существенно, одного порядка [10].

#### 1.13 Критериальные зависимости стационарного переноса

Перенос количества движения. Рассмотрим примеры стационарного переноса количества движения. Внутренней задачей гидродинамики является описание движения жидкостей и газов в трубах. Для стационарного горизонтального движения в трубах (отсутствуют критерии Фруда и гомохронности) критериальную зависимость (1.46) записывают в виде


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


Подставив выражение критерия Эйлера в (1.74), получим уравнение Дарси:


Можно показать [6], что в ламинарном режиме (Re 2100):


а в развитом турбулентном режиме для гладких труб (Re 10000) применяется зависимость Блазиуса:


Для шероховатых труб (Δ – высота выступов шероховатости) и переходного режима [6] используют зависимости общего вида:


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


Равномерное движение частиц обусловлено равновесием сил, действующих на частицу – тяжести, архимедовой и сопротивления среды [6]:


где Сх – коэффициент лобового сопротивления частицы.

С учетом, что потери давления при обтекании частицы равны отношению силы сопротивления к сечению частицы


получим из (1.79):


Таким образом, движение частицы сводится к зависимости коэффициента лобового сопротивления Сх от числа Рейнольдса. В ламинарном режиме (Re 2) движение частицы описывается законом Стокса,

в переходном (2 Re 500), –

а в турбулентном (500 Re 2105) Сх практически не зависит от Re и составляет Сх = 0,44.

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


Конкретный вид зависимости (1.81) для ламинарного режима найден Левеком [10]:


Для труб большой длины в установившемся ламинарном режиме, при (Реd/L) 20 величина Nu стремится к постоянному значению 3,695 теоретически полученному Нуссельтом.

Для турбулентного течения в трубах (Re 10000) при отношении L/d 50 в литературе используется уравнение:


Для газов последний множитель в уравнении (1.83) равен единице и Pr зависит только от атомности газа. Так для воздуха в этом случае получим:


Для стационарного теплопереноса при обтекании тел (внешняя задача) вид зависимости (1.83) сохраняется. Так при перпендикулярном обтекании коридорных и шахматных пучков труб при Re 1000 используется уравнение:


а для обтекания шахматных пучков труб при Re 1000 применяется уравнение:


Для газов последний множитель в уравнении (1.86) равен единице и Pr зависит только от атомности газа. Для воздуха в этом случае получим:


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


Конкретный вид зависимости (1.88), например, для адсорбции на насыпном активированном угле с размером частиц 1,7-2,2 мм при средней скорости потока 0,3-2 м/с имеет вид


Аналогичный вид зависимости с близкими численными значениями коэффициентов приведен [6] для адсорбции в кипящем слое активированным углем паров четыреххлористого углерода.

Если одновременно имеет место тепло- и массоперенос (например, при сушке твердых материалов), то для расчета процесса необходимо использовать критериальные уравнения для тепло- и массопереноса.

#### 1.14 Сушка как процесс тепломассопереноса

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


Рис. 1.7 Виды сушки материалов различного агрегатного состояния.


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

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

Влага из материалов может быть удалена различными способами: 1) механическим: 2) физико-химическим; 3) собственно сушкой.

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

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

На Рис. 1.8 показаны разные способы удаления влаги из материалов с их детализацией: механический способ удаления влаги, физико-химический и удаление влаги способом сушки. Обычно при сушке из материала удаляется поверхностная влага в начале процесса, а затем и внутренняя. В ряде производственных процессов из материала могут удаляться и органические растворители (ацетон, бензин и др.), например, в производствах химических волокон при сухом и мокром формовании, в производстве искусственной кожи, нефтепереработке. Многие виды сушилок и методы их расчета могут быть использованы как для удаления влаги, так и для процессов удаления органических растворителей.


Рис. 1.8 Способы удаления влаги из материалов.


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


где определяющим размером при вычислении критериев Sh и Re является длина поверхности испарения в направлении движения сушильного агента. В уравнение (1.90) входит параметрический критерий Гухмана:


где tc и tм – температуры сухого и мокрого термометров.

Величины коэффициентов А и n в уравнении (1.90) зависят от режима движения сушильного агента, определяемого критерием Рейнольдса. Эта зависимость в форме таблицы представлена ниже.

Таблица 1.1 Коэффициенты критериального уравнения сушки


При тепломассопереносе в процессе сушки уравнения тепло- и массопереноса взаимосвязаны. Поток тепла q, подводимый для удаления влаги из материала, может быть найден как произведение потока массы удаляемой влаги J на теплоту парообразования r:


или с учетом уравнений (1.29) и (1.32):


В периоде постоянной скорости сушки скорость сушки с поверхности материала F, с учетом соотношения (1.93) может быть найдена как:


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

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

Как следует из уравнения (1.94), скорость сушки в первом периоде можно определить как из уравнения тепло-, так и из уравнения массопереноса. Коэффициенты α и находят из соответствующих критериальных уравнений.

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

## Глава II

### Типы и свойства тепло- и влагоносителей

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


Рис. 2.1 Типы сушилок по виду теплоносителя (способу подвода тепла).

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

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

#### 2.1 Водяной пар

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

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

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

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

В таблице 2.1 приведены для ряда значений абсолютного давления параметры насыщенного водяного пара: температура конденсации, плотность, теплота парообразования, его энтальпия и энтропия. Более подробная таблица приведена в Приложении. Аналогичные таблицы приведены также в литературе по тепломассообменным процессам [1, 4, 9].

Таблица 2.1 Параметры сухого насыщенного пара

Р – давление; t к —температура; ρ – плотность; r —теплота парообразования; h —энтальпия; S – энтропия.


Теплота парообразования r – это количество тепла, необходимое для превращения 1 кг кипящей воды (при данном давлении) в сухой насыщенный пар. Энтальпия (теплосодержание) сухого насыщенного пара h равна энтальпии жидкости при температуре кипения tк плюс теплота парообразования r.


Аналогично энтропия сухого насыщенного пара S равна энтропии жидкости при температуре кипения tк плюс теплота парообразования r деленная на температуру кипения в градусах Кельвина.


Значение коэффициента теплоотдачи α, характеризующего теплоперенос от насыщенного пара к нагреваемому им материалу, в среднем составляет 1500 Вт/(м2·К). Коэффициент теплоотдачи значительно уменьшается в случае присутствия в нем инертных примесей, обычно воздуха и воды, поэтому их количество должно быть минимальным.

Расход греющего пара в процессе сушки Gг определяют по уравнению:


где Q – количество теплоты, отдаваемое паром (рассчитывается по уравнению теплового баланса – см. ниже); r – теплота парообразования (табл. 1.1); х – степень сухости пара.

Степень сухости пара х – это массовая доля сухого насыщенного пара во влажном паре. Степень сухости насыщенного пара, часто в среднем принимают равной 0,95, это означает, что на 100 кг влажного пара приходится 95 кг сухого насыщенного пара и 5 кг воды.

Пример 1. Определить температуру и расход насыщенного водяного пара, поступающего в калорифер при избыточном (манометрическом) давлении 0,5 МПа, если его степень сухости 95%, а количество передаваемого паром в калорифере тепла 10000 Вт.

Решение: По табл. 2.1 для абсолютного давления 0,1+ 0,5= 0,6 (МПа) находим температуру пара 158,1 oС.

Расход греющего пара определяем по уравнению (2.3):


Здесь r = 2095 кДж/кг, также находится по таблице 2.1.

#### 2.2 Горячий воздух

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

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

Влажный воздух характеризуется следующими параметрами: температурой t, точкой росы tр, влагосодержанием d, абсолютной и относительной влажностью φ, энтальпией h и энтропией S.

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

Абсолютная влажность – это масса водяного пара (кг), содержащаяся в 1 м3 влажного воздуха. Она равна плотности водяного пара ρп при температуре воздуха t и парциальном давлении пара рп.

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


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


Относительную влажность воздуха определяют специальным прибором-психрометром, либо по диаграмме Рамзина h-d для воздуха (см. ниже). Влагосодержание воздуха определяется отношением плотности пара к плотности абсолютно сухого воздуха при той же температуре:


Ряд значений d для температур от 0 до 100°С и относительной влажности φ от 5 до 100% приведены в табл. 2.2 (для среднегодового барометрического давления 99,3 кПа, т. е. 745 мм рт. ст. [1]).

Для установления аналитической связи между влагосодержанием и относительной влажностью φ запишем выражение (2.6) с учетом уравнения Менделеева-Клайперона:


Таблица 2.2 Влагосодержание воздуха d · 103 кг/кг при различных φ и t для давления 745 мм рт. ст. (99,3 кПа)


Подставляя мольную массу воды Мп = 18 и мольную массу абсолютно сухого воздуха Мсв = 29, с учетом (2.5) и закона Дальтона, согласно которому давление рсв равно разности общего давления влажного воздуха Р и парциального давления водяного пара в нем рп получим выражение:


Точка росы достигается при температуре tр, которую будет иметь воздух с влагосодержанием d в результате его охлаждения до состояния насыщения (φ = 1, в процентах – 100%, т. е. при температуре насыщения в данной точке tн).

Энтальпия h влажного воздуха – это количество содержащейся в нем теплоты, отнесенное к 1 кг сухого воздуха при той же температуре. Отсчет ведут от 0°С. Величина h равна сумме энтальпий 1 кг сухого воздуха и энтальпии для d кг водяного пара:


где hв = сcвt – энтальпия сухого воздуха, t – его температура и сcв – теплоемкость (обычно принимают среднюю, не зависящую от температуры 1,01 кДж/кг·К); hn —энтальпия водяного пара.

Энтальпия водяного пара определяется как сумма


где: hн – энтальпия насыщенного водяного пара, сп – теплоемкость перегретого водяного пара и tн – температура насыщения. Применяя предложенную Рамзиным линейную зависимость от температуры, имеем [1].


Эта зависимость применима и для перегретого пара и равна сумме теплосодержаний при 0 оС (2493 кДж/кг) и тепла перегрева от 0 до t оС. Аналогично определяется энтропия влажного воздуха:


Величину энтропии воздуха удобнее определять по Т-S диаграмме для воздуха [1].

Для определения параметров воздуха удобно пользоваться также диаграммой Рамзина [1] в координатах энтальпия h – влагосодержание d. Угол между координатными осями на такой диаграмме равен 135°, но для расчета пользуются вспомогательной осью d, проведенной под углом 90° к ординате. Диаграмма h-d для воздуха представлена на Рис. 2.2. В более крупном масштабе в единицах системы СИ эта диаграмма дана также в приложении. На ней нанесены изотермы – линии постоянных температур (они идут с наклоном относительно горизонтали) и пучок расходящихся кривых, построенных по уравнению (2.8), выходящих из точки t = – 273 °С и d = 0, каждая из которых соответствует постоянному значению относительной влажности φ.


Рис. 2.2. Диаграмма h-d для воздуха.

При температуре 99,4 °С давление насыщенного пара становится равным среднему барометрическому, т. е. 745 мм рт. ст., для которого и составлена диаграмма. Если температура воздуха выше температуры насыщения, то максимальное давление водяного пара будет равно барометрическому В и относительная влажность воздуха, согласно (2.5) определится как:


Тогда с учетом отношения (2.13) уравнение (2.8) примет вид:


Из уравнения (2.14) следует, что при температурах выше 99,4 °С кривые для φ не зависят от температуры. Поэтому эти кривые на диаграмме при t = 99,4 °С (температура кипения воды) имеют перелом и далее идут почти вертикально. Отклонение их от вертикали связано с тем, что плотность перегретого пара в некоторой степени зависит от температуры.

Линия насыщенного пара φ = 100 % является граничной между областью влажного воздуха (выше неё) и двухфазной областью насыщенного пара с капельками воды (ниже неё), где использование полученных выше формул для смеси сухого воздуха и водяного пара неправомерно. Пересечение изотерм (линий постоянных температур) с линией насыщенного пара φ = 100 % дает точку конденсации (точка росы) tр.

По уравнению парциального давления водяного пара [1]:


из центра координат на h-d диаграмме проведена наклонная прямая. Определяют парциальное давление пара по заданному его влагосодержанию d, проводя соответствующую вертикаль до пересечения с прямой линией, построенной по уравнению (2.15) и затем из точки пересечения проводят горизонталь вправо до шкалы рп (1 мм рт.ст. = 133,3 Па).

Каждая точка диаграммы определяется двумя параметрами по пересечению 2-х линий. Остальные параметры точки находят по диаграмме. Рассмотрим примеры определения параметров влажного воздуха с помощью h-d диаграммы.

Пример 2. Определить энтальпию, влагосодержание, точку росы влажного воздуха и парциальное давление водяного пара при t = 60 °С и φ = 40%.

Решение. По диаграмме h-d (Рис. 2.2) находим точку пересечения изотермы t = 60 °С с линией φ = 40%. Этой точке соответствуют: энтальпия h = 205 кДж/кг сухого воздуха, влагосодержание d = 0,055 кг/кг сухого воздуха, точка росы tр= 43 °С. Схема решения примера дана на Рис. 2.3, а. Влагосодержание для заданных условий по Таблице 2.2 составит d = 0,05448 кг/кг, а расчет энтальпии по уравнению (2.9) дает (кДж/кг):


По диаграмме h-d (Рис. 2.2) находим также точку пересечения влагосодержания d = 0,055 кг/кг с линией парциального давления водяного пара и по шкале справа находим парциальное давление пара рп = 60 мм рт.ст. Расчет парциального давления водяного пара по уравнению (2.15) для общего давления Р = 745 мм рт. ст. дает близкое значение


Пример 3. Определить относительную влажность воздуха, точку росы и парциальное давление водяного пара при t = 90°С и влагосодержании d = 0,07 кг/кг сухого воздуха.


Рис. 2.3. Схема решений примеров по диаграмме h-d для воздуха: а—к примеру 2, б—к примеру 3.


Решение. По диаграмме h-d (Рис. 2.2) находим точку пересечения изотермы t = 90 оС с линией постоянного влагосодержания d = 0,07 кг/кг сухого воздуха. Этой точке соответствует относительная влажность воздуха φ = 15% и точка росы tр= 47 оС. Схема решения дана на Рис. 2.3, б. Затем по диаграмме h-d (Рис. 2.2) находим точку пересечения влагосодержания d = 0,07 кг/кг с линией парциального давления водяного пара и по шкале справа находим парциальное давление пара рп = 75 мм рт.ст. Расчет парциального давления водяного пара по уравнению (2.15) при общем давлении Р = 745 мм рт. ст. дает близкий результат:


#### 2.3 Изменение состояний воздуха при сушке на h-d диаграмме

При конвективной сушке сушильный агент – воздух предварительно нагревается в калорифере от температуры t0 до температуры t1, необходимой для сушки (Рис. 2.4, а). Соответствующие значения энтальпии составят h0 и h1. При этом влагосодержание воздуха d остается постоянным, а относительная влажность φ – уменьшается. На диаграмме h-d этот процесс изображается вертикальным отрезком АВ (от точки А к точке В). Расход тепла в калорифере:


где L – расход сухого воздуха.

Процесс охлаждения воздуха на диаграмме также изображается вертикально, но в противоположном направлении. Если охлаждение воздуха идет до температуры насыщения tн, линия идет вниз до пересечения с линией φ =100% (отрезок ВС). Точка пересечения линий d = const и φ == 100% (точка С на Рис. 2.4, а) характеризует состояние воздуха в результате его охлаждения при d = const (точка росы). Изотерма, проходящая через эту точку, определяет температуру точки росы tр. Дальнейшее охлаждение воздуха ниже температуры точки росы (например, до температуры tп) приводит к конденсации из него части влаги и соответственно – к уменьшению его влагосодержания от d0 до dп. На диаграмме процесс охлаждения насыщенного воздуха совпадает с линией φ = 100% (кривая СЕ).

После нагрева воздуха в калорифере (линия АВ, Рис. 2.4, а) он поступает в сушилку. Если влага из материала будет испаряться только за счет тепла, передаваемого материалу воздухом, и отсутствуют теплопотери или дополнительный нагрев в сушилке (политропная сушка), то энтальпия воздуха после сушки h2 будет равна его энтальпии перед сушкой h1, так как все тепло, отданное воздухом на испарение влаги, возвращается обратно в воздух с удаляющимися из материала парами (адиабатическая сушка). Одновременно в сушилке понижается температура, увеличиваются влагосодержание и относительная влажность воздуха. Такой процесс носит название теоретического процесса сушки (h2 = h1 = h = const, линия ВD на Рис. 2.4, а).


Рис. 2.4 Изменение состояний воздуха на h-d диаграмме:

а – при нагреве (АВ) охлажденнии (ВАСЕ) и сушке (ВD)

б – линии сушки: адиабатическая (ВD), с теплопотерями (ВD′′) и с подводом тепла в сушилке (ВD′).


Для сушильной практики большое значение имеет понятие о теоретическом процессе адиабатического испарения в системе поверхность испарения – воздух. В этом процессе воздух только испаряет, но не нагревает влагу. Непосредственно над поверхностью испарения воды (а в равной степени и над поверхностью влажного материала в начальный период сушки) образуется слой насыщенного пара (φ = 100%), находящегося в равновесии с водой. Температура влаги при этом имеет постоянное значение, равное температуре мокрого термометра tм. Данная температура в процессе испарения не меняется, в то время как температура воздуха по мере его насыщения все время понижается, приближаясь в пределе к температуре мокрого термометра tм (при φ = 100%). Эту температуру, которую примет воздух в конце процесса насыщения, называют также температурой адиабатического насыщения. Поступающая в воздух испаренная влага W вносит в него некоторое количество тепла W∙c∙tм, поэтому адиабатический процесс охлаждения воздуха в этом случае происходит с повышением его энтальпии (h2>h1). С учетом расхода сухого воздуха L на испарение получим:


или


где с – теплоемкость воды.

Величина W/L – отражает увеличение влагосодержания воздуха в процессе его адиабатического охлаждения, оно равно (dм-d1), где dм—влагосодержание воздуха при его полном насыщении влагой при температуре мокрого термометра tм. Тогда при d1= 0 (уравнение оси ординат) получим


Уравнение (2.19) служит для нанесения на h-d диаграмму линий адиабатического насыщения воздуха.

Если, например, из точки М (см. Рис. 2.4, а) провести линию h1 = h2 = соnst до пересечения с осью ординат и отложить от точки пересечения R в масштабе, выбранном для энтальпий, отрезок RS, равный dм∙с∙tм , то полученная точка S и будет искомой, а прямая SМ будет линией адиабатического насыщения воздуха. Эту линию иногда называют также линией постоянной температуры мокрого термометра (tм = соnst) потому, что мокрый термометр, помещенный в воздух, насыщение которого происходит по линии адиабатического насыщения, будет показывать постоянную температуру tм. На h-d диаграммах для воздуха часто наносят пунктиром линии tм = соnst (см. на Рис. 2.4, а и приложение).

Приборы, состоящие из сухого и мокрого термометров (психрометры), широко используются в сушильной практике. Значение tм определяют с помощью термометра, нижний конец которого обернут влажной тканью («мокрый» термометр). Для получения надежных данных приток тепла излучением к нему должен быть минимален (экранирование при скорости воздуха около 5 м/с). По разности показаний «сухого» термометра (температура которого равна температуре воздуха tв) и мокрого термометра tм, пользуясь h-d диаграммой находят относительную влажность воздуха φ и другие его параметры (влагосодержание, энтальпия и парциальное давление водяного пара). Для этого по изотерме для температуры tм находят ее пересечение с линией φ = 100% (Рис.2.4, а). Далее двигаясь по линии постоянной температуры мокрого термометра (tм = соnst) доходят до пересечения с изотермой для температуры сухого термометра tв. В точке пересечения находят искомое значение φ и другие параметры.

Для более точного определения характеристик воздуха необходимо ввести поправку А на скорость движения воздуха u, м/с (обычно u>0,5 м/с) в месте установки психрометра [1] по уравнению:


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


где рн’– давление насыщенного водяного пара при температуре мокрого термометра; В – барометрическое давление, атм.

На Рис. 2.4, б показаны варианты процесса сушки от температуры t1 на входе до температуры t2 на выходе из сушилки. Линия ВD соответствует процессу в теоретической сушилке (h = соnst). Линии действительной сушки проходят либо выше нее ВD’, с повышением энтальпии (при подводе дополнительного тепла в сушилку), либо ниже ВD” с понижением энтальпии (при наличии теплопотерь в окружающую среду, на нагрев материала, на нагрев транспортных устройств). По линии теоретической сушки (адиабатического насыщения воздуха влагой) на h-d диаграмме для воздуха происходит изменение его параметров (температуры, влагосодержания и относительной влажности) за счет адиабатического испарения свободной поверхностной влаги в начальный первый период сушки.

Разность между температурой воздуха (сухого термометра) tс и температурой мокрого термометра tм, характеризует способность воздуха поглощать влагу из материала и носит название потенциала (движущей силы) сушки ε:


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

Более удобно определять потенциал сушки по ε-d диаграмме (Рис. 2.5), предложенной Г. К. Филоненко [2], дающей связь между ε, d и tм. Она построена в более простых, прямоугольных координатах.

Поскольку теплота испарения, теплоемкость влажного воздуха и диффузионный критерий Прандтля, в сущности, не зависят от давления, то линии адиабаты (h = соnst) и температуры мокрого термометра (tм = соnst) на h-d диаграмме для воздуха одни и те же для разных давлений. Это дает возможность ее использовать при давлениях, отличных от атмосферного, с корректировкой других параметров диаграммы.

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


Рис. 2.5 Диаграмма ε-d Г. К. Филоненко [2] для определения потенциала сушки.


Ниже рассмотрены примеры применения h-d диаграммы для воздуха при нахождении потенциала сушки и относительной влажности воздуха в сушилке.

Пример 4. Найти потенциал сушки для теоретической сушилки на входе ε1, выходе ε2 и среднее значение εср по следующим данным:

На входе в калорифер температура воздуха t0 = 22 °С, относительная влажность φ0 = 0,75; на выходе из сушилки температура воздуха t2 = 50 оС, относительная влажность φ2 = 0,5.

Решение. По известным параметрам t0 и φ0 находим на h-d диаграмме точку А (Рис 2.4, а) и влагосодержание d0 = 0,0125 кг/кг. По известным параметрам t2 и φ2 находим на диаграмме точку D и энтальпию h2 = 162 кДж/кг. Проводим через точку D линию, соответствующую процессу в теоретической сушилке (h = соnst) до пересечения с линией насыщения φ = 100% (точка М) и получаем tм = 38 оС. Продолжая линию h = соnst вверх до пересечения с вертикальной линией d0 = d1 = соnst (точка В, Рис. 2.4, а) и получаем t1 = 125 оС.

Находим потенциал сушки в точке В:

ε1 = t1 – tм = 125-38 = 87 (оС).

Находим потенциал сушки в точке D:

ε2 = t2 – tм = 50-38 = 12 (оС).

Среднее значение εср находим по уравнению:


Пример 5. Определить по показаниям психрометра tм = 70 °С и tв = 85 °С относительную влажность воздуха в сушилке, если давление в ней 700 мм рт.ст., а скорость воздуха 2 м/с.

Решение. Т. к. h-d диаграмма построена для давления 745 мм рт.ст., то расчет проводим аналитически. Находим поправку А по уравнению (2.20) на скорость движения воздуха в сушилке u, в месте установки психрометра:


По таблице свойств насыщенного водяного пара (Приложение) находим давление пара рн’= 0,3177 при температуре мокрого термометра tм = 70 оС. Для температуры tв = 85 оС давление насыщенного водяного пара рн находим по формуле линейной интерполяции табличных данных:


Определяем парциальное давление водяного пара в движущемся воздухе по уравнению (2.21):


Окончательно находим относительную влажность воздуха в сушилке по уравнению (2.5):


#### 2.4 Топочные газы

Топочные (дымовые) газы в смеси с атмосферным воздухом широко используются при сушке различных материалов, в том числе органических веществ и пищевых продуктов. Многие материалы, например, песок, глина, твердое топливо, неорганические соли и т. д. высушивают при довольно высоких температурах – от 300 до 800 °С и выше. Для сушки этих материалов наиболее рационально использовать топочные газы, разбавляя их до нужной температуры атмосферным воздухом.

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

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

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

Топочные (дымовые) газы получают при сжигании газообразного (природный газ), жидкого (мазут, нефть) или твердого (уголь, торф, сланцы, дрова) топлива в топках и смешивают их в специальных камерах смешения с атмосферным воздухом для получения смеси заданной температуры. Горючая часть твердого и жидкого топлива состоит из углерода С, водорода Н и серы S. Элементарный состав рабочей массы топлива составляет 100 %. В общем случае:


Сера в топливе может быть в трех видах: органическая So (в составе сложных органических соединений), колчеданная Sк (в составе соединений с металлами типа FeS2) и сульфатная (в процессах горения не участвует). Зола топлива Ар – это смесь негорючих минеральных соединений. Золу и влагу топлива Ар + Wр называют внешним балластом, а кислород и азот Ор + Nр – внутренним балластом [5], т. к. они не участвуют в горении топлива. При отсутствии балластных соединений получаем горючую массу топлива, обозначая (г), а при полном отсутствии влаги в топливе пользуются понятием сухой массы топлива – (с).

Характеристики основных видов топлива. Угли марки А (антрацитные): с высокой степенью углефикации (С = 90-93%) и малым выходом летучих (Vг=2… 9%) относят к антрацитам. Угли марки Б (бурые): с неспекающимся коксом, выходом летучих Vг>40 % и QнР = 10,5… 15,9 МДж/кг.

Торф – наиболее молодое ископаемое топливо с большим выходом летучих порядка 70%, высокой влажностью Wp =40… 50% и низкой теплотой сгорания QнР= 8,38… 10,47 МДж/кг.

Сланцы – топливо с большой зольностью Ар = 50… 60%, повышенной влажностью Wp = 15…20%, низкой QнР =5,87…10 МДж/кг при высокой Qнг = 27,2… 33,5 МДж/кг. Сланцы относятся к легко воспламенимым топливам вследствие высокого содержания водорода Нг = 7,5… 9,5% и летучих Vг=80… 90%. Сланцы являются местным топливом, так как вследствие большой зольности и влажности транспортировка его на дальние расстояния неэкономична.

Большое значение в топливном балансе России имеют природные газы. Основной горючей составляющей их является метан СН4 (80-98%). Теплота сгорания сухого природного газа Qнг = 33,52… 35,61 МДж/м3. Состав газообразного топлива задается в процентах по объему горючих составляющих СН4, СО, Н2, углеводородов СnН2n+2 и негорючих газов СО2, О2, N2, Н2О. Характеристики различных видов топлив [5] приведены в приложении.

Образующиеся при сгорании топочные газы по химическому составу состоят из кислорода О2, азота N2, окиси CО и двуокиси CO2 углерода, сернистого газа SO2 и водяных паров H2O. Состав топочных газов зависит от вида топлива и количества воздуха, подводимого в топку для его сжигания (первичный воздух) и подмешиваемого к продуктам сгорания для понижения их температуры до заданной (вторичный воздух). Особенно вредны соединения серы, содержащиеся в топливе, т. к. они в соединении с влагой коррозируют металлические конструкции сушилок и ухудшают качество высушиваемых материалов. Поэтому применение сернистых топлив для сушки очень ограничено.

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

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

При обычных условиях сжигания топлива в котельных и других промышленных установках теплота конденсации водяных паров продуктов сгорания не может быть использована, так как их температура значительно превышает точку росы. Поэтому при характеристике твердого и жидкого топлива эту теплоту исключают и получают значение его низшей теплоты его сгорания QнР. В значении высшей теплоты сгорания топлива QнР эта теплота учитывается. Связь между ними выражается формулой [5]:


где: коэффициент 0,025 учитывает теплоту парообразования, МДж/кг; величина 9Нр – количество воды, образующейся при сгорании водородной составляющей топлива (2: 9) по уравнению:


Низшая теплота сгорания углей QнР в зависимости от их месторождения, марки и сорта колеблется в пределах от 10000 до 28500 кДж/кг, жидких топлив (мазут, гудрон) – от 30000 до 40000 кДж/кг, газообразных топлив (природный газ) – от 38000 до 50 000 кДж/кг. Для сравнения различных топлив используется понятие условного топлива, имеющего теплоту сгорания 29330 кДж/кг (7000 ккал/кг). Таким образом, наибольшей теплотой сгорания обладает наиболее дешевое газообразное топливо, удобное в применении. Оно и наиболее широко используется.

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

Наиболее точная эмпирическая формула для расчета теплоты сгорания твердых и жидких топлив предложена Д. И. Менделеевым:


Низшая теплота сгорания сухого газообразного топлива (МДж/м3) может быть определена [5] по формуле:


Так как в расчетах газовых сушилок расход газа берется в кг, то теплоту сгорания сухого газообразного топлива необходимо пересчитать на МДж/кг. Для этого все члены уравнения (2.26) необходимо поделить на плотность смеси газов ρсм, находимую по правилу аддитивности:


где у – весовое содержание газа в смеси, ρ – его плотность.

С учетом этого для уравнения (2.27) получим в МДж/кг:


На основании уравнения Менделеева-Клайперона для идеальных газов плотность любого газа при заданной температуре Т и давлении Р рассчитывают по формуле:


где ρ0 – плотность газа при нормальных условиях (Р0 = 760 мм рт.ст и Т0 = 273,15 К). Расчеты упрощаются для природного газа, т. к. он состоит в основном из метана (более 90 %), а для него ρ0 = 0,72 кг/м3.

Высшая теплота сгорания сухого газообразного топлива (МДж/кг) может быть определена по формуле:


#### 2.5 Расчет процесса горения топлива

Горение топлива – сложный физико-химический процесс взаимодействия топлива с окислителем (кислородом воздуха), сопровождающийся интенсивным выделением тепла и быстрым подъемом температуры.

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

Скорость химического взаимодействия в процессе гомогенного горения может быть выражена через изменение концентраций реагирующих веществ C в единицу времени τ c учетом константы скорости реакции k. Для простейшего случая – реакция первого порядка выражается уравнением:


Интегрирование уравнения (2.31) позволяет получить кинетическое уравнение процесса горения вида:


где С0 – начальная концентрация вещества.

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

Скорость реакции в сильной степени зависит от температуры. Как известно, эта зависимость выражается законом Аррениуса:


где k0 – постоянная, определяемая экспериментально; R – газовая постоянная, кДжДмоль ∙К); T– термодинамическая температура, К; Е – энергия активации.

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


Реакция для водорода (2.24) записана выше. Из этих уравнений количество кислорода для сгорания 1 кг элементов топлива составит для С – 32/12 =2,67; для S – 1; для Н – 32/4=8.

Состав воздуха в объемных % обычно принимается [1] 21 % кислорода и 79 % азота. Мольная масса воздуха составит:

М = 0,21∙32 + 0,79∙28 = 28,8 кг/моль.

Состав воздуха в массовых % составит 0,21∙32/28,8=0,233 % кислорода и 0,79∙28/28,8 = 0,768 % азота. С учетом доли кислорода необходимое количество воздуха для сгорания 1 кг топлива составит для С —2,67/0,233 = 11,5; для S (как и для О) – 1/0,233=4,3; для Н – 8/0,233 =34,5. С учетом выражения элементарного состава топлива в процентах по весу, получим теоретическое количество абсолютно сухого воздуха для сжигания 1 кг твердого или жидкого топлива составит:


Теоретическое количество абсолютно сухого воздуха Lо, необходимое для сжигания 1 кг твердого или жидкого топлива, может быть рассчитано также по приближенной формуле [5]:


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