Математическое моделирование процесса сушки

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

2.1 Материальный баланс сушилок

Материальный баланс сушилки играет большую роль в расчётах процесса сушки. Сначала необходимо определить начальное щ1 и конечное щ2 влагосодержание продукта.

Влажность можно представить как отношение общего количества влаги W в материале к сумме W+Gсух.

Или же как отношение общего количества влаги W в материале к количеству абсолютно сухого веществаGсух.

Следовательно, влажность в % можно определить как:

Абсолютная влажность в % выражается так:

При необходимости связать общую и абсолютную влажность используют следующие формулы:

Ведём следующие обозначения:

G1 - количество влажного продукта, входящего в сушилку, кг/ч;
G2 - количество высушенного материала, выходящего из сушилки, кг/ч;
W - количество влаги, удаляемой из продукта.

Количество абсолютно сухого вещества можно определить по формуле (5):

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

Таким же образом определим количество продукта, поступающего на сушку:

По формуле (8) определим количество влаги, которое удаляется из продукта:

Так же можно воспользоваться формулой (9):

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

Введём следующие обозначения:

V - объём сушильной камеры, м3;

- время сушки, ч.

Найдём объём сушильной камеры по влагосодержанию, кг/(м3Чч):

где W - количество удаляемой влаги, кг.

2.2 Баланс влаги в сушилке

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

Где L - количество абсолютно сухого воздуха, который необходим для процесса сушки, кг/ч;

d1 и d2 - влагосодержание сушильного агента (воздуха) на входе и выходе из сушильной камеры, г/кг сухого воздуха.

Удельный расход сухого воздуха на 1 кг испаряемой влаги равен l=L/W. Отсюда определим удельный расход сухого воздуха, кг/кг:

Калорифер нагревает воздух от температуры t0 до температуры t1. При этом количество влаги в сушильной камере остаётся неизменным, т.е. d0=d1. Исходя из этого запишем формулу (12) в виде:

2.3 Тепловой баланс сушилки

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

Рис. 1 - Принципиальная схема сушилки [1]

Чтобы составить тепловой баланс необходимо ввести обозначения:

- средняя удельная теплоемкость в (Дж/кг?К): сушильного агента; влаги; которая удаляется в процессе сушки; самого продукта и транспортировочных устройств сушилки;

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

- температура продукта перед входом в сушилку и на выходе из нее, °С;

- вес транспортных механизмов, кг;

- температура транспортного механизма на входе и выходе из сушильной установки, °С;

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

Баланс тепла выражают так:

Тепловой баланс выражается следующим равенством:

Будем решать уравнение (14) по подводимому теплу Q.

Исходя из этого уравнения можно сказать, на что тратится тепло, которое подводится в сушилку:

1) расход тепла с уходящим сушильным агентом:

2) на процесс испарения влаги из продукта:

3) на нагрев высушенного материала:

4) на нагрев транспортных механизмов:

5) в атмосферу .

Исходя из этого уравнение (14) может принять вид:

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

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

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

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

Теперь известное значение  подставляем в уравнение (17), получаем:

здесь  - сумма удельных расходов тепла в подогревающем воздух устройстве и сушильной установке.

Введем обозначение разницы величин:

Исходя из этого, уравнение (18) приводится к виду:

Рассчитаем количество тепла , которое вносится сушильным агентом:

здесь  - тепло воздуха из атмосферы,

- тепло, которое получил воздух в воздухоподогревателе.

Из (19) , а значит

Вынесем  за скобки и, подставив , поучим уравнение вида:

Учитываем, что , окончательно выражение принимает вид:

Данное уравнение является уравнением теплового баланса сушильных камер конвективного типа [1].

3. Моделирование процесса радиационно-конвективной сушки материала

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

При испарении жидкости с поверхности капиллярно-пористого тела имеют место три случая [3]:

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

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

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

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

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

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

Такой подход наиболее приемлем, так как в массе находящегося в аппарате материала присутствуют объекты, обрабатывающиеся в режиме как первого, так и второго периода одновременно [2].

В качестве модели материала принят плоский диск диаметром D и толщиной 2R. Ось X ориентирована вдоль поверхности, а Z перпендикулярна поверхности материала. Начало координатной системы выбрано на оси симметрии его поперечного сечения.

Период условно постоянной скорости сушки

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

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

Дифференциальные уравнения переноса для теплоносителя могут быть записаны в следующем виде [1].

Уравнение переноса массы:

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

Показано весьма слабое влияние фактора поперечного потока массы на процессы тепло- и массообмена в процессах испарения. Таким образом, в уравнениях (22) и (23) можно положить:

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

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

Общее выражение для плотности потока влаги в капиллярно-пористом теле (в направлении оси Z) записывается [2] в виде соотношения:

В таком случае уравнение сохранения массы для образца может быть записано в форме уравнения (22), а уравнение сохранения энергии для осушаемого плоского двумерного объекта, в условиях пренебрежения термическим сопротивлением тонкой пленки влаги, может быть записано в классическом виде [3]:

Для решения уравнений (22), (23), (24) тепломассопереноса в первом периоде сушки необходимо сформировать условия однозначности - краевые условия.

С учетом фазового перехода при испарении воды со свободной поверхности, граничные условия третьего рода для уравнений (22) и (23) примут вид:

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

±R - координаты Z для верхней и нижней поверхности материала при толщине 2R;

п - индекс для поверхности объекта.

Следует отметить, что граничные условия 3 рода для уравнения (24) при указанном способе осушки материала записываются в форме уравнения (26).

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

Период убывающей скорости сушки

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

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

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

математический моделирование сушка уравнение

В начальный момент времени температуру tн и влагосодержание Uо(t) материала принимаем постоянными:

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

Граничные условия могут быть записаны в виде:

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

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

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

Заключение

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

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

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

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

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