Математическое моделирование электротехнологического процесса озонирования яйцескладов птицефабрик

Аннотации

Санитарная обработка яиц — необходимая мера борьбы с бактериями, грибками и другими микроорганизмами. От качества санитарной обработки инкубационных яиц зависят такие показатели, как выводимость яиц и сохранность суточного молодняка. Ведущие ученые нашей страны доказали высокую эффективность применения озона для обработки инкубационных яиц. Для получения положительного результата от применения перспективного способа санитарной обработки инкубационных яиц озоном необходимо создать равномерную концентрацию озона по всему объему помещения яйцесклада, так как снижение концентрации озона от заданного значения в любой точке пространства может привести к недостаточному воздействию на болезнетворные организмы и как следствие? к снижению процента вывода цыплят. Значительное превышение концентрации от заданного значения может привести к губительному воздействию на зародыш. В результате математического моделирования электротехнологического процесса озонирования яйцескладов птицефабрик получена математическая модель электротехнологического процесса озонирования яйцескладов птицефабрик, учитывающая разложение озона на поверхности яиц и стен яйцесклада, а также обосновывающая требуемое управляющее воздействие при различных исходных данных, таких как: требуемая концентрация озона в яйцескладе, количество яиц, геометрические размеры яйцесклада, подача вентилятора, параметры конструкции разрядного устройства, расположение электроозонатора в помещении яйцесклада Ключевые слова: ОЗОНАТОР, МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ЭЛЕКТРОТЕХНОЛОГИЧЕСКОГО ПРОЦЕССА ОЗОНИРОВАНИЯ ЯЙЦЕСКЛАДОВ ПТИЦЕФАБРИК, ЯЙЦЕСКЛАД, КОНЦЕНТРАЦИЯ ОЗОНА.

UDC 62−932.2.

Technical sciences.

MATHEMATICAL MODELING OF ELECTRO TECHNOLOGICAL OZONIZATION OF EGG STORES OF POULTRY FARMS.

Voloshin Aleksandr Petrovich.

senior lecturer.

Kuban State Agrarian University, Krasnodar, Russia.

Sanitization of eggs is an essential way to fight bacteria, fungi and other microorganisms. Hatchability of eggs and the safety of day-old chicks are dependent on the quality of eggs processing. Leading scientists of our country have proved high efficacy of ozone application for processing of hatching eggs. To obtain a positive result by this method of sanitizing hatching eggs ozone, it is necessary to create a uniform concentration of ozone around the egg store volume. Decrease in ozone concentration from the set point may result in insufficient exposure to pathogens and because of this, may reduce hatchability. Significant excess of ozone concentration from the set point can kill the embryo. Because of mathematical modeling of electro eggs ozone treatment process, there was a mathematical model of the process of electroozonation of egg stores of poultry farms. This model takes into account decomposition of ozone on the surface of eggs and decomposition of ozone on the walls of an egg store. This mathematical model proves the desired control action at different initial data, such as: the required concentration of ozone in the egg store, the number of eggs, egg store geometric dimensions, the fan power, the design parameters of the discharge unit, the location of the ozone generator in the room.

Keywords: OZONE GENERATOR, MATHEMATICAL MODEL OF OZONIC TECHNOLOGICAL PROCESS IN THE EGG STORE, EGG STORE, CONCENTRATION OF OZONE.

Doi: 10.21 515/1990;4665−121−098

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

Возмущающими воздействиями являются: температура воздуха t_возд, площадь стен помещения яйцесклада S_ст, площадь поверхности яиц S_я. Такое количество параметров представляет большую сложность для равномерного распределения концентрации озоно-воздушной смеси в каждой точке пространства помещения яйцесклада. Снижение концентрации от заданного значения в любой точке пространства может привести к недостаточному воздействию на болезнетворные организмы и как следствие? к снижению процента вывода цыплят. Значительное превышение концентрации от заданного значения может привести к губительному воздействию на зародыш.

Рисунок 1 — Функциональная схема яйцесклада, как объекта управления.

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

Рисунок 2 — Схема предлагаемого электротехнологического процесса озонирования яйцескладов птицефабрик: 1 — разрядное устройство электроозонатора, 2 — вентилятор, 3 — датчик концентрации озона в непосредственной близости от электроозонатора, 4 — датчик концентрации озона на удалении от электроозонатора, 5 — контейнеры с яйцами.

Функциональная схема предлагаемого электротехнологического процесса озонирования яйцескладов птицефабрик представлена на рисунке 3. На схеме, электроозонатор является управляющим устройством УУ формирующим входное воздействие на объект управления ОУ (помещение яйцесклада).

ПР — это программный регулятор, необходимый для формирования сигнала управления, который подаётся на УУ. В зависимости от этого УУ выдаёт необходимое количество озоновоздушной смеси Q_вент с определённым значением концентрации озона на выходе С_оз.вых.. Два датчика концентрации озона представлены на функциональной схеме как измерительные органы ВО1 и ВО2.

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

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

Для того чтобы узнать какая концентрация озона будет в удалённой точке помещения яйцесклада С_оз.я. при различных значениях подачи вентилятора Q_вент., и подаваемого на электроозонатор напряжения U_оз. необходимо произвести математическое моделирование.

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

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

При электроозонировании яйцесклада птицефабрики концентрация озона C зависит от: потерь через дверные, оконные проёмы или ворота C_ПР, потерь через вентиляцию C_ВЕНТ, потерь, связанных с разложением озона на стенах яйцесклада C_СТ, потерь, связанных с разложением озона на инкубационных яйцах C_Я.

C = C_ОЗ — С_ПР — С_ВЕНТ — С_СТ — С_Я, (1).

где, С_ОЗ — концентрация озона на выходе из электроозонатора, мг/м³;

C_ПР — потери озона через дверные, оконные проёмы или ворота, мг/м³;

C_ВЕНТ — потери через вентиляцию, мг/м³;

C_СТ — потерь связанных с разложением озона на стенах яйцесклада, мг/м³;

C_Я — потерь связанных с разложением озона на инкубационных яйцах, мг/м³.

Изменение концентрации озона в помещении от времени описывается уравнением:

(2).

где, C — объемная концентрация озона, %;

q — расход озона, м³/ч;

L_ВЕНТ — вентиляционный расход, м³/ч;

D — коэффициент диффузии озона м²/ч;

S_ПР — площадь открытого проема (двери и т. д.), м²;

L_ПР — характерное расстояние от электроозонатора до открытого проема, м;

V — объем помещения, м³.

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

(3).

Сделаем замену переменной С. Пусть, тогда,, имеем:

(4).

(5).

(6).

(7).

(8).

Проинтегрируем обе части. Имеем:

(9).

Отсюда, выразим переменную v:

(10).

(11).

Имеем:

(12).

Отсюда, общее решение дифференциального уравнения имеет вид:

(13).

(14).

Данное уравнение позволяет построить график зависимости изменения концентрации озона C_ОЗ от времени t в воздухе при нормальных условиях (рисунок 4):

Рисунок 4 — График зависимости изменения концентрации озона CОЗ от времени t в воздухе при нормальных условиях.

Как видно, с увеличением времени концентрация озона в помещении яйцесклада C_ОЗ падает. Через час она уменьшится в 4 раза. Постоянная времени ф = 25 минут.

Данное уравнение не учитывает поглощение озона стенами помещения С_ОЗ.СТ. и числом яиц N, а также не учитывает распределение озона на расстоянии от электроозонатора.

Распределение концентрации озона C_ОЗ по объему помещения яйцесклада и ее изменение во времени t описывается уравнением диффузии:

инкубационный яйцесклад озонирование электротехнологический.

(15).

где,.

С_ОЗ — концентрация озона, мг/м³;

D — диффузии озона, м²/ч;

Q — объемный расход вещества, м³/ч;

v — скорость воздушного потока в помещении, м/ч;

x, y, z — пространственные координаты, м;

t — время, ч;

V_см — объем озоновоздушной смеси, м³/ч.

Q = q-L_ВЕНТ· С_ОЗ, (16).

где, q — подача воздуха, м³/ч;

L_ВЕНТ — вентиляционный расход, м³/ч.

Рассмотрим одномерный случай распространения озона. Фундаментальное решение однородного уравнения с постоянными, не зависящими от x и t, D.

(17).

(18), где:

(19).

В случае произвольного распределения, общее решение данного уравнения:

(20).

(21).

В итоге получаем уравнение диффузии в направлении оси х:

(22).

Данное уравнение позволяет построить график зависимости концентрации озона от координаты x. Уравнение не учитывает поглощение озона стенами помещения и числом яиц N.

Найдём массу озона m_ОЗ в помещении яйцесклада. Плотность озона с составляет [17]:

Масса озона в этом случае составит:

(23).

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

(2.24).

где, а, b, с — полуоси эллипсоида (рисунок 5), м;

k_Я — коэффициент поглощения озона яйцами (k = 0,0001);

N — количество яиц, шт.

Рисунок 5 — Представление яйца, как эллипсоида с тремя полуосями.

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

В случае, когда электроозонатор установлен сверху, при высоте потолка в 3 м, расстояние между электроозонатором и контейнером составляет 1,34 м (0,3 м — высота корпуса электроозонатора, высота контейнера составляет 1,36 м).

Поэтому количество обрабатываемых озоном яиц будет находиться по формуле:

(25).

где, х_ОЗ-К — расстояние от электроозонатора до контейнера, м;

N_ЯР — количество ярусов (для контейнера КТПЯ-05.00.00.00−01 N_ЯР = 3 шт.), шт.;

N_K — количество контейнеров, шт.;

N_Я — количество яиц в одном ярусе (для контейнера КТПЯ-05.00.00.00−01 N_Я = 1080 шт.), шт.

Подставляя (25) в (24) получим:

(26).

Если электроозонатор расположен на одной из стен яйцесклада, то формула (25) будет иметь вид:

(27).

Подставляя (27) в (24) получим:

(28).

Помещение яйцесклада имеет форму параллелепипеда, то есть площадь поверхности стен яйцесклада найдём по формуле:

S_СТ = 2((a₁· c₁) + (b₁· c₁) + (a₁· b₁)), (29).

где, a₁, b₁, c₁ — стороны параллелепипеда (длина, ширина и высота помещения соответственно), м.

Масса озона, поглощенная стенами яйцесклада будет равна:

(30).

где, k_СТ — коэффициент поглощения озона стенами (k_СТ = 0,042).

Таким образом, остаточная масса озона m_ОЗ.ОСТ. с учетом поглощения части озона яйцами и стенами равна:

(31).

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

(32).

Остаточная концентрация озона за время t определяется функцией:

(33).

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

(34).

Заключение

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

• 1. Лунин В. В. Физическая химия озона. / В. В. Лунин, М. П. Попович, С. Н. Ткаченко // - М: Изд-во МГУ, 1998. — 480 с.
• 2. Лытнев А. С. Результаты экспериментальных исследований модернизированной конструкции разрядного устройства пластинчатого типа / А. С. Лытнев, А. П. Волошин // Материалы VI международной научно-практической конференции «Актуальные проблемы энергетики АПК»: / Под общ. ред. Трушкина В. А. — Саратов: ООО «ЦеСАин», 2015. — 327 с. С. 30−33.
• 3. Нормов Д. А. Математическое моделирование физических процессов в электроозонаторах барьерного типа / Д. А. Нормов // КубГАУ.: «Труды Кубанского государственного аграрного университета» № 3(12), 2008. — С 231−235.
• 4. Овсянников Д. А. Математическая модель нагрева разрядного устройства электроозонатора / Овсянников Д. А., Николаенко С. А. — Материалы пятой Всероссийской конференции «Энергои ресурсосберегающие технологии и установки». — Краснодар, 2007. — 4 с.
• 5. Оськин С. В. Электротехнологии в сельском хозяйстве: учебник для студентов вузов / С. В. Оськин. — Краснодар: КубГАУ, 2016. — 501с.
• 6. Пат. РФ № 2 417 159, МПК С2 С01В13/11 (2006.01) Электроозонатор / Д. А. Овсянников, С. А. Николаенко, С. С. Зубович, А. П. Волошин, Д. С. Цокур; заявитель и патентообладатель КГАУ. — № 2 009 126 863 заявл. 13.07.2009; опубл. 27.04.2011. Бюл. № 2. — 5 с.
• 7. Пат. РФ № 2 429 192, МПК С2 С01В13/11 (2006.01) Электроозонатор / Д. А. Овсянников, С. А. Николаенко, С. С. Зубович, А. П. Волошин, Д. С. Цокур; заявитель и патентообладатель КГАУ. — № 20 091 330 067 заявл. 2.09.2009; опубл. 20.09.2011. Бюл. № 26. — 6 с.
• 8. Сторчевой В. Ф. Математическое моделирование стационарных процессов ионизатора-озонатора. Природообустройство № 2. — М.: РГАУ-МСХА им. К. А. Тимирязева. 2012. С. 78−82.

References

• 1. Lunin V.V. Fizicheskaja himija ozona. / V.V. Lunin, M.P. Popovich, S.N. Tkachenko // - M: Izd-vo MGU, 1998. — 480 s.
• 2. Lytnev A.S. Rezul’taty jeksperimental’nyh issledovanij modernizirovannoj konstrukcii razrjadnogo ustrojstva plastinchatogo tipa / A.S. Lytnev, A.P. Voloshin // Materialy VI mezhdunarodnoj nauchno-prakticheskoj konferencii «Aktual'nye problemy jenergetiki APK»: / Pod obshh. red. Trushkina V.A. — Saratov: OOO «CeSAin», 2015. — 327 s. S. 30−33.
• 3. Normov D.A. Matematicheskoe modelirovanie fizicheskih processov v jelektro-ozonatorah bar’ernogo tipa / D.A. Normov // KubGAU.: «Trudy Kubanskogo gos-udarstvennogo agrarnogo universiteta» № 3(12), 2008. — S 231−235.
• 4. Ovsjannikov D.A. Matematicheskaja model' nagreva razrjadnogo ustrojstva jelek-troozonatora / Ovsjannikov D.A., Nikolaenko S.A. — Materialy pjatoj Vserossijskoj konferencii «Jenergoi resursosberegajushhie tehnologii i ustanovki». — Krasnodar, 2007. — 4 s.
• 5. Os’kin S.V. Jelektrotehnologii v sel’skom hozjajstve: uchebnik dlja studentov vuzov / S.V. Os’kin. — Krasnodar: KubGAU, 2016. — 501s.
• 6. Pat. RF № 2 417 159, MPK S2 S01V13/11 (2006.01) Jelektroozonator / D.A. Ov-sjannikov, S.A. Nikolaenko, S.S. Zubovich, A.P. Voloshin, D.S. Cokur; zajavitel' i patentoobladatel' KGAU. — № 2 009 126 863 zajavl. 13.07.2009; opubl. 27.04.2011. Bjul. № 2. — 5 s.
• 7. Pat. RF № 2 429 192, MPK S2 S01V13/11 (2006.01) Jelektroozonator / D.A. Ov-sjannikov, S.A. Nikolaenko, S.S. Zubovich, A.P. Voloshin, D.S. Cokur; zajavitel' i patentoobladatel' KGAU. — № 20 091 330 067 zajavl. 2.09.2009; opubl. 20.09.2011. Bjul. № 26. — 6 s.
• 8. Storchevoj V.F. Matematicheskoe modelirovanie stacionarnyh processov ionizatora-ozonatora. Prirodoobustrojstvo № 2. — M.: RGAU-MSHA im. K.A. Timirjazeva. 2012. S. 78−82.