Математическая модель системы вентиляции с переменным расходом воздуха вдоль вертикальной оси реактора в процессе компостирования

Размещено на http://www.allbest.ru/

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

При построении модели были приняты следующие допущения.

1. За время моделирования перемешивание субстрата отсутствует.

2. Теплообмен через стенки реактора отсутствует.

3. Направление движения воздуха - от основания к верхушке вдоль оси биореактора.

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

5. Теплопроводность и влагопроводность субстрата не рассматриваются вследствие их малого влияния на рассматриваемый процесс.

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

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

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

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

, (1)

где и - масса соответственно в точках и , кг;

- температура, С;

- теплоемкость, Дж/(кг*0С);

- плотность, кг/м3.

В элементарном объеме субстрата происходят процессы массопереноса между субстратом и газовой фазой :

- испарение воды с поверхности субстрата в газовую фазу;

- конденсация воды из газовой фазы на поверхность субстрата.

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

- конвективный теплообмен между поверхностью субстрата и газовой фазой ;

- изменение энтальпии субстрата (газовой фазы) при испарении (конденсации) воды .

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

, (2)

где - изменение массы влаги в субстрате или газовой среде (с соответствующим знаком), кг;

- коэффициент массопереноса влаги между субстратом и газовой фазой, кг воды/(с*м3);

- влагосодержание субстрата, кг воды/кг субстрата;

- равновесное влагосодержание субстрата при известных температуре и влагосодержании газовой фазы, кг воды/кг субстрата.

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

, (3)

где - удельная теплота испарения воды, Дж/кг воды.

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

, (4)

где - коэффициент теплообмена между субстратом и газовой фазой, Дж/(с*м3*К);

- температура газовой фазы, 0С;

- температура субстрата, 0С.

Выражение для ,, определено экспериментально в виде эмпирического выражения:

, (5)

где - экспериментально определенные коэффициенты.

, (6)

где - экспериментально определенные коэффициенты.

, (7)

где - экспериментально определенные коэффициенты.

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

, (8)

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

- плотность субстрата, кг субстрата/м3;

- влагосодержание газовой фазы, кг воды/кг воздуха.

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

Тогда поток воздуха системы вентиляции в точке х можно представить в виде G(x), причем для непрерывности потока необходимо выполнение условия:

, (9)

где G(0) - поток воздуха в основании реактора.

Таким образом, выражение (9) описывает поток воздуха в точке вертикальной оси биореактора для традиционной системы вентиляции при G(0)=G, ; для системы вентиляции с боковой подачей воздуха при , .

Из (8) имеем:

, (10)

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

Запишем массовый поток в точках и :

, (11)

. (12)

Подставляя (11), (12) в (10), получаем:

(13)

Производя предельный переход при , , получаем:

. (14)

Из (8) имеем:

. (15)

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

. (16)

Производя предельный переход при , получаем:

. (17)

Из (8) имеем:

, (18)

где - температура подаваемого воздуха системы вентиляции, oС.

Запишем массовый поток в точках и :

, (19)

. (20)

Подставляя (19), (20) в (18), получаем:

. (21)

Производя предельный переход при , и учитывая, что

, , где - удельная теплоемкость сухого газа, Дж/(кг *К); - удельная теплоемкость водяного пара, Дж/(кг*К), получаем:

. (22)

Из (8) имеем:

. (23)

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

. (24)

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

. (25)

Из (14), (17), (22), (25) получаем математическую модель разработанной системы вентиляции биореактора в виде системы уравнений, описывающих тепловой и массовый баланс двухфазной среды:

. (26)

Учитывая, что традиционная система вентиляции является частным случаем системы с боковой подачей воздуха, для которого [4], получаем математическую модель традиционной системы вентиляции биореактора в виде системы уравнений, описывающих тепловой и массовый баланс двухфазной среды:

(27)

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

Данная статья подготовлена в рамках выполнения научного проекта 16-48-230441 а(р) «Математическое моделирование процессов, протекающих в автоматизированной установке для круглогодичного производства органических удобрений в условиях Краснодарского края», финансируемого РФФИ и администрацией Краснодарского края.

Литература

1. MacDonald, L. Physical and Mathematical Modelling of the Composting Process / L. MacDonald - University of Guelph, 1995.

2. Московец А.Л. Постановка задачи и анализ оптимального управления процессом компостирования/ А.Л.Московец, С.В.Усатиков, Д.Л.Пиотровский //депонированная рукопись № 1926-B2004. - 06.12.2004

3. Пиотровский Д.Л. Автоматизация производства органических удобрений/ Д.Л. Пиотровский, Т.Г.Шарапкина// Автоматизация. Современные технологии. 2004. - № 7. - С. 9-11

4.Пиотровский Д.Л. Теоретические основы построения автоматических систем управления процессами производства органических компостов: диссертация … доктора технических наук: 05.13.06 . - Краснодар, 2007

References

вентиляция биореактор компостирование

1. MacDonald, L. Physical and Mathematical Modelling of the Composting Process / L. MacDonald - University of Guelph, 1995.

2. Moskovec A.L. Postanovka zadachi i analiz optimal'nogo upravlenija processom kompostirovanija/ A.L.Moskovec, S.V.Usatikov, D.L.Piotrovskij //deponirovannaja rukopis' № 1926-B2004. - 06.12.2004

3. Piotrovskij D.L. Avtomatizacija proizvodstva organicheskih udobrenij/ D.L. Piotrovskij, T.G.Sharapkina// Avtomatizacija. Sovremennye tehnologii. 2004. - № 7. - S. 9-11

4. Piotrovskij D.L. Teoreticheskie osnovy postroenija avtomaticheskih sistem upravlenija processami proizvodstva organicheskih kompostov: dissertacija … doktora tehnicheskih nauk: 05.13.06 . - Krasnodar, 2007

Размещено на Allbest.ru