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

бесплатно 0
4.5 83
Математическое моделирование как метод оптимизации процессов. Расчет сушилок, баланс влаги. Моделирование процесса радиационно-конвективной сушки. Уравнение переноса массы. Период условно-постоянной скорости. Градиент влагосодержания и температуры.


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

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

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

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

1. Математическое моделирование как метод оптимизации процессов

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

Традиционные методы расчета ХТП, основанные на учете при вычислениях упрощенных механизмов их протекания. Абсолютно не удовлетворяют современным требованиям [3].

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

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

При разработке математического описания объекта в настоящее время применяется два основных подхода: системный и эмпирический [1].

2. Расчет сушилок

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

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

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

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

Или же как отношение общего количества влаги W в материале к количеству абсолютно сухого ВЕЩЕСТВАGСУХ.

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

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

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

(4)

Ведем следующие обозначения: G1 - количество влажного продукта, входящего в сушилку, кг/ч;

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

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

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

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

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

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

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

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

Введем следующие обозначения: V - объем сушильной камеры, м3;

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

Найдем объем сушильной камеры по влагосодержанию, кг/(м3?ч): (10) где W - количество удаляемой влаги, кг.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Баланс тепла выражают так: Приход: Расход: С сушильным агентом С сушильным агентом

С обрабатываемым материалом С готовой продукцией

С влагой, которая удаляется из продукта С удаленной из продукта влагой

Физическое тепло транспортных механизмов Физическое тепло нагретых транспортных механизмов

Подводимое тепло Потери тепла в атмосферу

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

(14)

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

(15)

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

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

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

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

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

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

(16)

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

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

(17)

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

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

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

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

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

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

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

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

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

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

(21)

Данное уравнение является уравнением теплового баланса сушильных камер конвективного типа [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], движущей силой процесса массоуноса является градиент влагосодержания и температуры.

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

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

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

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

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

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

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

Заключение

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

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

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

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

Список литературы
1. Касаткин А.Г. Основные процессы и аппараты химической технологии. - М.: Химия, 1971. - 784 с.

2. Кафаров В.В., Глебов М.Б. Математическое моделирование основных процессов химических производств. - М.: Высшая школа, 1991. - 400 с.

3. Гартман Т.Н., Клушин Д.В. Основы компьютерного моделирования химико-технологических процессов. - М.: Академкнига, 2006. - 416 с.

Размещено на
Заказать написание новой работы



Дисциплины научных работ



Хотите, перезвоним вам?