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

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

Рубрика Химия
Вид автореферат
Язык русский
Дата добавления 03.02.2018
Размер файла 1,0 M

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

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

АВТОРЕФЕРАТ
диссертации на соискание ученой степени доктора технических наук
Моделирование и оптимизация процессов смешивания сыпучих материалов

05.17.08 - Процессы и аппараты химических технологий

05.13.01 - Системный анализ, управление и обработка информации (в химических технологиях, нефтехимии)

Баранцева Елена Александровна

Иваново 2010

Работа выполнена в Государственном образовательном учреждении высшего профессионального образования «Ивановский государственный энергетический университет имени В.И. Ленина».

Научный консультант: доктор технических наук, профессор Мизонов Вадим Евгеньевич.

Официальные оппоненты: доктор технических наук, профессор Зайцев Анатолий Иванович;

доктор технических наук, профессор Першин Владимир Федорович;

доктор технических наук Волынский Владимир Юльевич.

Ведущая организация: Российский химико-технологический университет им. Д.И. Менделеева», г. Москва.

Защита состоится «1» марта 2010 г. в 10 часов на заседании совета по защите докторских и кандидатских диссертаций Д 212.063.05 при ГОУВПО «Ивановский государственный химико-технологический университет» по адресу: 153000, г. Иваново, проспект Ф. Энгельса, 7

С диссертацией можно ознакомиться в библиотеке ГОУВПО «Ивановский государственный химико-технологический университет» по адресу: 153000, г. Иваново, проспект Ф. Энгельса, 10.

Автореферат разослан « » 2010 г.

Ученый секретарь совета, доктор физико-математических наук, профессор Зуева Г.А.

Общая характеристика работы

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

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

Все отмеченное и определило цель настоящей работы, которая выполнялась в рамках ФЦП «Интеграция» (2.1 - А118 Математическое моделирование ресурсосберегающих и экологически безопасных технологий), НТП «Развитие инфраструктуры научно-технической и инновационной деятельности высшей школы и ее кадрового потенциала» №2/05 (37788) и международными договорами о научном сотрудничестве между ИГЭУ и Горным институтом г.Алби, Франция, Ченстоховским политехническим институтом, Польша, университетом г. Веспрем, Венгрия, и исследовательским центром Tel-Tek, Норвегия.

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

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

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

Задачи исследования

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

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

3. Разработать рациональные (оптимальные) способы воздействия на перемешиваемые компоненты, обеспечивающие повышение качества смесей и/или производительности смесителей.

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

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

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

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

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

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

4. Для эмпирического обеспечения модели лопастного смешивания выполнено экспериментальное исследование миграции частиц при прохождении плоской лопасти через слой сыпучего материала и показано, что вероятность частиц остаться в зоне прохождения лопасти зависит только от высоты слоя и не зависит от угла установки лопасти, а отношение вероятностей перейти назад и вперед зависит только от угла установки лопасти и не зависит от высоты слоя. Степень захвата материала плоской лопастью составляет 0,2…0,3, что ниже теоретически полученного оптимального по скорости перемешивания значения 0,4…0,5 и определяет направление совершенствования процесса путем создания профилированных лопастей.

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

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

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

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

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

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

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

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

6. Разработанные подходы к построению математических моделей процессов в дисперсных средах и программно-алгоритмическое обеспечение моделирования процессов смешивания нашли применение в практике выполнения промышленных и исследовательских проектов в горном институте г.Алби, Франция, Ченстоховском политехническом институте, Польша, университете г. Веспрем, Венгрия, и исследовательском центре Tel-Tek, Норвегия; новые конструкции смесителей и рекомендации по модернизации действующих смесителей внедрены в ООО «Полимепластбетон», ООО «Инженерный центр», г. Ярославль, ЗАО «Резинотехника», г. Ярославль, ООО «Спецтехника», г. Кострома с реальными техническими и экономическими эффектами.

Апробация результатов работы.

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

5-я Международная НК «Теоретические и экспериментальные основы создания новых высокоэффективных химикотехнологических процессов и оборудования», Иваново, 2001, 10-й Российско-польский семинар Теоретические основы строительства, Москва-Иваново, 2001, Международные НТК „Состояние и перспективы развития электротехнологии (Бенардосовские чтения). Иваново, 2001, 2003, 2007, Международная школа молодых ученых «Методы кибернетики в технологиях, экономике и управлении производством», Иваново, 2002 г., Международные конференции «Математические методы в технике и технологиях»: ММТТ-16, Санкт-Петербург, 2003; ММТТ-18, Казань, 2005; ММТТ-19, Воронеж, 2006; ММТТ-20, Ярославль, 2007; ММТТ21, Саратов, 2008, Международная НК «Теоретические основы создания, оптимизации и управления энерго- и ресурсосберегающими процессами и оборудованием», Иваново, 2007, 13-я Международная Плесская конференции по нанодисперсным магнитным жидкостям, Плес, 2008, The 1-st Int. Congr. (the 2-nd French Congr.) on Tracers and Tracing Mehtods. 29-31 May 2001, Nancy, France, The 4th International Conference for Conveying and Handling of Particulate Solids, Budapest, Hungary, 2003, The 4th Powder Science and Technology French Forum, France, Compiegne, 2004, The World Congress on Particle Technology, Orlando, USA, 2005, The Fifth International Conference for Conveying and Handling of Particulate Solids CHoPS-05, Sorrento, Italy, 2006, The Int. Conf. “Science and Technology of Particles”, Albi, France, 2007, The International Symposium on Reliable Flow of Particulate Solids IV (RELPOWFLO IV), Tromso, Norway, 2008.

Публикации. По теме диссертации опубликована 61 печатная работа, в том числе, 1 монография, 17 работ в изданиях, предусмотренных перечнем ВАК, 2 патента на полезные модели, 3 положительных решения о выдаче патентов на полезные модели, 1 свидетельство о государственной регистрации программы для ЭВМ.

Структура и объем работы. Диссертация состоит из введения, 6-х глав, основных выводов, списка использованных источников (236 наименований) и приложения.

Содержание работы

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

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

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

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

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

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

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

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

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

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

В заключение главы приведены детализированные задачи исследования.

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

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

(1)

где cj - относительное его содержание в j-ом объеме. При полностью однородной по выделенным локальным объемам смеси все cj=1/N и ?=0. В диссертации выполнено исследование влияния объемов, по которым определяется качество смеси, и ограниченного числа проб на отклонение характеристики смеси, определяемой по уравнению (1). Подчеркнуто, что технологически оправданным является объем пробы, равный объему смеси, отпускаемой потребителю (например, объем многокомпонентной лекарственной таблетки, для которой в большинстве случаев безразлично, как компоненты распределены по ее объему, если их содержание удовлетворяет стандарту). Однако и в этом случае может возникнуть неопределенность. Например, стиральный порошок отпускается в объеме пакета, но в пакете смесь должна быть устойчиво равномерно распределена по его объему, так как содержимое пакета используется по частям.

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

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

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

Пусть состояние системы характеризуется набором вероятностей каждого из состояний Si , i= 1, 2, 3, …, m. Очевидно, что если рассматривается эволюция фиксированной порции частиц, то пространство состояний является полным и

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

(2)

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

, (3)

где ' обозначает транспонирование вектора или матрицы.

Будем наблюдать за процессом, то есть регистрировать эти вероятности состояния, через конечные промежутки времени ?t. Тогда моменты фиксации состояния процесса будут tk=(k-1)?t, k=1,2 …, где величина k - номер перехода - является дискретным целочисленным аналогом времени. Таким образом, обе переменные процесса j и k становятся целочисленными.

Очевидно, что в любой момент времени вектор состояния полностью характеризует смесь. Если цепь является цепью Маркова, то два последовательных состояния Sk и Sk+1 связаны рекуррентной матричной формулой:

Sk+1 = PSk, (4)

где Р - матрица переходных вероятностей или матрица переходов, имеющая следующий вид:

. (5)

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

, , j=1,2,…,m. (6)

Свойства матрицы Р в зависимости от моделируемого процесса и модельного представления о нем могут быть весьма разнообразными. В большинстве случаев вероятности покинуть ячейку вправо и влево не равны друг другу. Предположим, что в случайной миграции частиц преобладают переходы вправо, то есть pj+1,j> pj-1,j. Во многих случаях удобно выделить в этих вероятностях симметричную составляющую, определяющую так называемый диффузионный (чисто случайный) перенос частиц, и направленную только вправо и «неуравновешенную» составляющую, определяющую конвективный перенос. На рис.1б показан пример такого выделения для второй ячейки цепи. Симметричная диффузионная составляющая p12=d, а конвективная __ p32-p12=v. Поскольку сумма всех переходных вероятностей равна единице, то вероятность остаться рассчитывается как p22=1-p12-p32=1 - 2d - v.

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

(7)

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

(8)

(9)

где D - коэффициент макродиффузии, V - скорость конвективного переноса.

Для многих моделей смешивания особую роль играет размещаемая в конце цепи ячейка, моделирующая коллектор выходящих из рабочей зоны смесителя частиц (абсорбирующая ячейка). Частицы могут попасть в эту ячейку из предыдущей ячейки, но возвратиться обратно из нее не могут, то есть вероятность остаться в ней равна единице. В ряде случаев (в частности, для непрерывного смешивания) коллектор удобно вынести за пределы цепи, оставив в ней только «рабочие» ячейки. В этом случае пространство состояний оказывается неполным, и сумма элементов в последнем столбце матрицы (5) будет меньше единицы. Важной характеристикой процесса является кинетика выхода в коллектор единичной порции материала, первоначально помещенного в первую ячейку. Если последняя m-ая ячейка в цепи на рис.1а относится к коллектору, то накопление в ней материала описывается величиной Smk, k=1,2, …, но более распространенной характеристикой является доля материала, выходящая в нее на каждом переходе, которая может быть рассчитана как

qk= Smk+1 - Smk = Sm-1k pm,m-1 (10)

Получающуюся гистограмму называют распределением времени пребывания (РВП) частиц в цепи, имея в виду, что номер перехода k есть целочисленный аналог времени. Очевидно, что при k>? ?qk >1. РВП характеризуется моментами распределения, среди которых наиболее важную роль играют момент первого порядка - среднее время (число переходов) пребывания частиц в рабочих ячейках цепи

(11)

и центральный момент второго порядка - дисперсия РВП в ненормированной и нормированной форме

, (12)

Моменты более высоких порядков редко используются в практике моделирования.

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

Sk+1=РSk + Sfk (13)

где Sfk - вектор-столбец подачи такого же размера, что и вектор состояния Sk. В практике моделирования чаще всего встречаются следующие виды Sfk:

- однократная импульсная подача в 1-ю ячейку (ввод трассера)

Sfk =[1 0 … 0] при k=0, Sfk =[0 0 … 0] при k=1,2 (14)

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

непрерывная постоянная подача в 1-ю ячейку

Sfk =[1, 0, …,0], при k=1,2,…,z (15)

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

- непрерывная постоянная подача в 1-ю ячейку с наложенной на нее переменной составляющей

Sfk =[S01(k),0,…,0] при k=1,2,…,z. (16)

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

В ряде модельных случаев возможен непосредственный расчет асимптотического состояния, минуя переходный процесс. В цепи с полным пространством состояний, в которой все состояния достижимы (периодический смеситель с постоянной загрузкой) асимптотический вектор состояния S? есть собственный вектор матрицы Р, соответствующий собственному числу ?=1. Для того, чтобы это асимптотическое состояние было равномерным, необходимо, чтобы матрица Р удовлетворяла условию

, i=1,2,…,m (17)

то есть чтобы сумма переходных вероятностей равнялась единице в каждой строке матрицы. Такая матрица может быть названа «несегрегационной»; любая другая матрица в принципе не обеспечит равномерного распределения даже при неограниченной продолжительности процесса.

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

Sk= (I-P)-1Sf (18)

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

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

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

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

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

На рис.2,3 показано влияние параметров микромасштабного перемешивания на кинетику формирования качества смеси. Из графиков видно, что диффузионная компонента благоприятно влияет на кинетику и предельно достижимое качество смеси, а конвективная (сегрегационная) компонента резко их ухудшает, причем равномерное распределение достижимо только при v=0; во всех остальных случаях неоднородность смеси проходит через минимум, а затем она вновь начинает расслаиваться.

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

(19)

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

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

Sk+i= PDSk+i-1, i=1,2,…,L, Sk+L+1= PBSk+L (20)

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

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

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

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

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

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

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

На рис.8 показано влияние параметра ? лопасти и структуры лопастного аппарата (рис.7) на формирование качества смеси для ключевого компонента, не склонного к сегрегации. Из рис.8,а следует, что доля захватываемого лопастью материала имеет при прочих равных условиях оптимальное значение, разное для различных комбинаций лопастей, но группирующееся около величины 0,5…0,6. При этом наилучшее качество смеси обеспечивает комбинация лопастей по схеме С рис.7. Это подтверждается и кинетикой процесса при ?=0,5, показанной на рис.8,б, из которого следует, что схема С обеспечивает не только наилучшее качество смеси, но и скорость достижения его достаточно высоких показателей (малых ?).

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

Ситуация меняется практически на противоположную, если ключевой компонент имеет склонность к сегрегации в смеси. Об этом свидетельствуют данные рис.9, полученные при v=0,1 в матрице PD. Для сравнения на рисунке нанесена жирной линией кинетика процесса при полном отсутствии лопастного перемешивания. Качество смеси, пройдя через максимум, соответствующий минимуму ?, вновь начинает ухудшаться, и неоднородность смеси асимптотически может стать даже больше исходной из-за скопления ключевого компонента внизу рабочего объема. При действии лопастей качество смеси стабилизируется на определенном достаточно высоком уровне, хотя и пульсирует относительно среднего значения с большей или меньшей амплитудой.

Наилучшее качество обеспечивается лопастями, установленными по схеме А рис.7; схема С дает наихудшее качество и большие его пульсации около среднего установившегося значения.

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

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

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

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

(21)

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

Было установлено, что определяющим фактором, влияющим на VRR, является отношение времени пребывания Тf ?Тm к периоду пульсаций Тin, где среднее по расходу время пребывания определяется как

Tf=M/Q (22)

где М - масса загрузки смесителя материалом, Q - массовый расход смеси через него.

Некоторые обобщенные результаты расчетов показаны на рис. 10. Требуемое технологией значение VRR достигается при Тin/Тm ?1,5, что позволяет по Тf =M/Q и заданной производительности найти требуемую массу материала в смесителе и его объем. При дисперсии РВП больше 0,75 ее влияние на VRR незначительно.

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

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

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

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

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

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

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

Глава завершается описанием нового подхода к совершенствованию процесса смешивания сегрегирующих компонентов - оптимального управления потоком сегрегирующего компонента. При начальной загрузке всего ключевого компонента в верхнюю ячейку смесителя периодического действия начальный вектор его состояния имеет вид Sfi =[1 0 0 … 0] при i=0 и Sfi =[0 0 0 … 0] при i>0. Предложено реализовывать распределенную во времени подачу ключевого компонента - его непрерывную засыпку в смеситель в течение k первых переходов. Такой подаче соответствует вектор Sfi =[1/j 0 0 … 0] при 0?i?j и Sfi =[0 0 0 … 0] при i>j, то есть в течение j первых переходов в первую ячейку на каждом переходе добавляется доля 1/j общей единичной загрузки ключевого компонента. На рис.13 показана кинетика формирования качества смеси при различных значениях j, причем рассматривать изменение величины ? имеет смысл только при i>j, когда в смесителе находится уже весь ключевой компонент. Таким образом, показано, что существует оптимальная программа загрузки ключевого компонента в работающий смеситель непрерывного действия, обеспечивающая наилучшее качество смеси склонных к сегрегации компонентов.

В четвертой главе описаны двумерные модели процесса смешивания, которые необходимо привлекать, когда одномерные модели в принципе не могут описать интересующие эффекты. Главным образом это относится к учету поперечной неоднородности потока и сегрегации материала в поперечном направлении. Принципы построения двумерной модели легко проиллюстрировать на примере периодического лопастного смесителя, показанного на рис.14. Сформировав двумерную цепь ячеек размером mxn, ячейки необходимо последовательно пронумеровать по столбцам и построить вектор-столбец состояния S размером (mxn)x1. На пронумерованной таким образом сетке матрицы переходных вероятностей размером (mxn)x(mxn) строятся по такой же схеме, что и одномерной модели: в столбце, соответствующем выбранной ячейке, в строках с номерами ячеек, куда разрешены переходы, следует разместить вероятности этих переходов. В частности, для схемы рис.14 матрица циркуляции РС и матрица конвективной диффузии РD имеют вид

, (23)

Для визуализации эволюции состояния вектор состояния удобно преобразовывать в матрицу, соответствующую сетке ячеек. Расчетные исследования данного смесителя показали, что оптимальное по скорости перемешивания значение степени захвата лопастью материала ? составляет 0,5…0,7, а оптимальная высота лопасти зависит от склонности ключевого компонента к сегрегации v и увеличивается с ее ростом.

Значительное внимание в работе уделено двумерным моделям непрерывного лопастного перемешивания. Расчетная схема процесса и его ячеечная модель показаны на рис.15. Рабочий объем смесителя, заполненный материалом и пересекаемый расположенными в шахматном порядке вращающимися на валу лопастями, разбит на n слоев высотой ?y в вертикальном направлении и m секций длиной ?x в горизонтальном направлении, то есть представлен двухмерной сеткой nxm ячеек. Высота ячейки равна высоте лопасти, а ее длина - длине лопасти. Средний уровень материала занимает (n - 1) нижних слоев.

На рис.15 показано состояние процесса, когда с материалом взаимодействуют лопасти с нечетными номерами. При моделировании лопастного перемешивания величина ?t связана со скоростью вращения лопастей ? и определяется соотношением ?t=?/?.


Подобные документы

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

    автореферат [1,2 M], добавлен 31.07.2009

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

    курсовая работа [211,9 K], добавлен 04.01.2009

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

    реферат [51,9 K], добавлен 21.10.2010

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

    дипломная работа [798,0 K], добавлен 16.04.2014

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

    контрольная работа [1,1 M], добавлен 18.01.2014

  • Процессы химической технологии. Разработка схемы химико-технологического процесса. Критерии оптимизации. Топологический метод и ХТС. Понятия и определения теории графов. Параметры технологического режима элементов ХТС. Изучение стохастических процессов.

    лекция [46,2 K], добавлен 18.02.2009

  • Обзор и анализ существующих методов оптимизации химико-технологических процессов. Определение параметров уравнения Аррениуса. Определение оптимальной температуры. Расчёт зависимости оптимальной скорости химической реакции от степени превращения.

    курсовая работа [498,1 K], добавлен 18.06.2015

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

    контрольная работа [1,8 M], добавлен 04.09.2009

  • Определение понятия и свойств полимеров. Рассмотрение основных видов полимерных композиционных материалов. Характеристика пожарной опасности материалов и изделий. Исследование особенностей снижения их горючести. Проблема токсичности продуктов горения.

    презентация [2,6 M], добавлен 25.06.2015

  • Особенности технологии изготовления полимерных материалов, основные параметры процессов переработки. Методы формования изделий из ненаполненных и наполненных полимерных материалов. Методы переработки армированных полимеров. Аспекты их применения.

    реферат [36,4 K], добавлен 04.01.2011

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.