Авторефераты по всем темам  >>  Авторефераты по техническим специальностям  

На правах рукописи

ЧЕРВЯКОВ Виктор Михайлович

Теоретические основы методов расчета роторных аппаратов с учетом нестационарных гидродинамических течений

05.02.13 - Машины, агрегаты и процессы

(химическая промышленность)

Автореферат

диссертации на соискание ученой степени

доктора технических наук

Тамбов 2007

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

Научный консультант	доктор технических наук, профессор Воробьев Юрий Валентинович
Официальные оппоненты	Заслуженный деятель науки и техники РФ, доктор технических наук, профессор Холпанов Леонид Петрович; доктор технических наук, профессор Зимин Алексей Иванович; доктор технических наук, профессор Першин Владимир Федорович
Ведущая организация	ГОУ ВПО Московский государственный университет инженерной экологии (МГУИЭ), г. Москва

Защита диссертации состоится л24 октября 2007 г. в 14 часов 30 минут на заседании диссертационного совета по присуждению ученой степени доктора технических наук Д 212.260.02 в Тамбовском государственном техническом университете по адресу: г. Тамбов, ул. Ленинградская, 1, ауд. 60.

Отзывы на автореферат в двух экземплярах, скрепленные гербовой печатью, просим направлять по адресу: 392000, г. Тамбов, ул. Советская, 106, ТГТУ, ученому секретарю диссертационного совета Д 212.260.02.

С диссертацией можно ознакомиться в библиотеке университета.

Автореферат разослан л    2007 г.

Ученый секретарь диссертационного совета, доцент

В.М. Нечаев

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы. Роторные аппараты с модуляцией потока обрабатываемой среды с высокой эффективностью используют в различных гидромеханических, химических, тепло-массообменных процессах, например, в процессах диспергирования, растворения, экстракции и т.д. Они характеризуются низкой удельной энерго- и металлоёмкостью при высоком качестве получаемого готового продукта. Высокая степень воздействия на обрабатываемую среду объясняется развитой турбулентностью, интенсивной акустической импульсной кавитацией, большими сдвиговыми напряжениями, гидравлическими ударами и другими механическими воздействиями. Теоретическое и экспериментальное исследование, внедрение роторных аппаратов в промышленность провели отечественные учёные М.А. Балабудкин, А.А. Барам, А.М. Балабышко, В.И. Биглер, П.П. Дерко, А.И. Зимин, Г.Е. Иванец, Е.А. Мандрыка, В.А. Плотников, М.А. Промтов, О.А. Трошкин, В.М. Фридман, В.Ф. Юдаев и многие другие.

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

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

Диссертационная работа выполнялась в соответствии с координационным планом АН СССР по направлению ТОХТ код 2.27.1.4.14 1991-1995гг; МНТП Ресурсосберегающие технологии машиностроения 1996-2000 гг; НТП Научные исследования высшей школы в области химической технологии 2003-2005гг.

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

Для достижения этой цели были поставлены и решены следующие задачи:

- проведение комплексных теоретических и экспериментальных исследований особенностей нестационарного течения потоков несжимаемой и сжимаемой жидкостей в модуляторе роторного аппарата и их сопоставление;

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

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

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

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

Научная новизна

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

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

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

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

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

- математическая модель нестационарного течения в радиальном зазоре между цилиндрическими ротором и статором, позволяющая определить пик мощности в момент пуска роторного аппарата и таким образом обоснованно подобрать электродвигатель привода;

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

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

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

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

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

Правовая защищенность разработок обеспечивается 16 авторскими свидетельствами СССР и патентами РФ на изобретения.

Результаты исследований, предложенные методики расчета и разработанные на их основе конструкции роторных аппаратов использованы на Мичуринском ОАО Прогресс, ОАО Тамбовполимермаш, Борисоглебском заводе ОАО Патроны, ООО Грибановский машиностроительный завод, ОАО Тамбовский завод Комсомолец им. Н.С. Артемова, ОАО РЖД Локомотивное ДЕПО г. Тамбов для приготовления эффективной смазочно-охлаждающей жидкости, на Котовском заводе ОАО КЛКЗ при производстве эмали белой ПФ-115, на  НПФ Лионик г. Москва при производстве сухих концентратов натуральных напитков. Получен  реальный экономический эффект около 18 млн. руб.

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

Апробация работы. Основные положения диссертации доложены и обсуждены на Всесоюзной научно-технической конференции Роль молодых исследователей и конструкторов химического машиностроения в реализации целевых комплексных программ и проблем (Москва, 1983); Всесоюзной научно-технической конференции Новые разработки в области ультразвуковой техники и технологии и опыт их применения в машиностроении (Новосибирск, 1989); Всесоюзном научном симпозиуме Акустическая кавитация и проблемы интенсификации технологических процессов (Одесса, 1989); Всесоюзной научно-технической конференции Биотехника-89 (Москва, 1989); Всесоюзной научно-технической конференции НТП в химмотологии топлив и смазочных материалов (Днепропетровск, 1990); 4 Всероссийской научной конференции Динамика ПАХТ (Ярославль, 1994); Всероссийской научно-технической конференции Математическое моделирование в научных исследованиях (Ставрополь, 2000); 4, 5 и 6 Международных конференциях Математическое моделирование физических, экономических, социальных систем и процессов (Ульяновск, 2001, 2003, 2005); 15 и 18 Международных научных конференциях Математические методы в техники и технологиях, ММТТ - 15 и 18 (Тамбов, 2000; Казань, 2005); Международной научно-технической интернет-конференции Фундаментальные и прикладные проблемы технологии машиностроения Технология-2002 (Орёл, 2002); 2 Международной научно-практической конференции Актуальные проблемы экологии (Караганда, 2003); 4 Международной конференции Актуальные проблемы современной науки (Самара, 2003); 5 Международной конференции Наука и образование (Белово, 2004); 1 и 2 Всероссийских научных конференциях Математическое моделирование и краевые задачи (Самара, 2004, 2005); 2 и 3 Международных научно-практических интернет-конференциях Энерго- и ресурсосбережение ЦXXI век (Орёл, 2004, 2005); Международной научной конференции Энерго-ресурсосберегающие технологии и оборудование, экологически безопасные производства (Иваново, 2004); 6 Международной научно-практической конференции Методы и алгоритмы прикладной математики в технике, медицине и экономике (Новочеркасск, 2006).

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

Структура и объём работы. Диссертация включает введение, семь глав, основные выводы и результаты, список литературы (503 наименования) и четыре приложения. Работа изложена на страницах основного текста, содержит 126 рисунков.

___________________________________________

Автор выражает искреннюю благодарность д.т.н., профессору Юдаеву В.Ф. за неоценимую помощь при проведении данного исследования.

Основное содержание работы

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

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

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

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

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

a) б)	в)
Рис.1. Конструктивная схема роторного аппарата: а) с цилиндрическим ротором, б) с коническим ротором, в) схема модулятора.

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

Уравнение нестационарного движения несжимаемой среды в канале ротора рассматривается на основании классических уравнений Навье-Стокса и неразрывности в цилиндрических координатах, при этом система координат связана с вращающимся каналом, вследствие чего в указанное уравнение входит в явном виде центробежная сила инерции. Для упрощения уравнений сделаны следующие допущения: скорость течения зависит только от радиальной координаты и времени. Кроме того, ввиду осевой симметрии имеем . Для приведения к безразмерному виду уравнения движения считаем, что , , . С учетом сделанных допущений скорость в канале как функция координат r и t выглядит следующим образом

,  (1)

где , . (2)

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

- выражения (2) приведены к виду

,.  (3)

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

Отметим, что при принятом масштабе радиальной скорости в нашем случае, критерий Eu=0,5.

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

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

,  (4)

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

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

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

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

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

а)	б)
Рис.2 Зависимость амплитуды модуля отрицательного ускорения: а) от критерия Kк, =0,1; б) от симплекса  ,  Kк =6,7.

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

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

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

Из анализа рис. 2а следует, что с уменьшением симплекса величина модуля ускорения увеличивается. Это можно объяснить тем, что при увеличении ширины канала ротора и уменьшении его длины инертность массы жидкости в канале уменьшается, т.е. она быстрее приводится в движение при открывании канала статора и набирает большую скорость за время открывания.

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

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

В работе определена нижняя граница применения модели течения несжимаемой жидкости по критерию в зависимости от и .

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

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

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

, (5)

а уравнение неразрывности для сжимаемой жидкости

.  (6)

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

.  (7)

После подстановки (7) в (6) и представления полученного уравнения и выражения (5) в безразмерном виде имеем систему уравнений

.  (8)

Проведя оценку малости величин системы (8) и  учитывая, что в данном конкретном случае , получаем уравнение гиперболического типа течения сжимаемой жидкости в виде

. (9)

Для аналитического решения введена новая безразмерная функция

.

Относительно этой функции получено уравнение

. (10)

Для решения уравнения (10) применен метод Фурье

.

Граничные условия для R(r) имеют вид

; .  (11)

Начальные условия для данного случая

.  (12)

С учетом начальных и граничных условий в результате решения уравнения (9) получено следующее выражение

  (13)

Для определения закономерностей течения сжимаемой среды в канале статора использовано уравнение неразрывности в виде (4), в предположении, что для малосжимаемой жидкости выполняются условия

.

С использованием компьютерного моделирования и пакета программ Maple 9.5 установлено влияние критерия и симплекса на особенности течения сжимаемой жидкости. Некоторые результаты представлены на рис.3.

Рис.3 Влияние Kк и на амплитуду модуля отрицательного ускорения: M=0,02; аp=0,002м; hp=0,01м.

При значении величина возрастает на большую величину (при=const), чем при . Таким образом, увеличивать нецелесообразно, кроме того, возрастает металлоемкость ротора. Изменение модуля амплитуды отрицательного ускорения проявляется значительно только, примерно, со значений . На основании полученных результатов можно рекомендовать определять параметры роторного аппарата, работающего в условиях гидравлического удара, при значениях и . При увеличении симплекса в 3 раза с 0,1 до 0,3 (при ), возрастает ~ в 3,4 раза, а при увеличении критерия - в 2 раза с 10 до 20 (при ) возрастает в 4 раза. Таким образом, при постоянном расходе через аппарат для повышения эффективности его работы выгоднее увеличивать угловую скорость его вращения в пределах исследуемых параметров.

Граница разграничения применимости моделей течения несжимаемой и сжимаемой жидкости определена из условия возникновения гидравлического удара, а именно

.  (14)

Следует отметить, что в (14) фазовая скорость звука в двухфазной среде определяется в зависимости от газосодержания (27). Из этого условия определены значения граничных параметров и .

На рис.4 показана форма графиков ускорения течения среды при наличии и отсутствии гидравлического удара.

График представленный на рис. 4а соответствует модели течения несжимаемой среды, а на рис.4б и 4в - модели сжимаемой среды при выполнении условия (14). Экспериментальное подтверждение этих результатов представлено в 5 главе.

a)	б)
в)	Рис.4. Зависимость ускорения течения жидкости в модуляторе роторного аппарата от времени: а=0,05;  zp=zc=20; =0,1; l/ap=5 и значениях Kк: a) 0,3; б) 1,1; в) 1,3.

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

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

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

,  (15)

.

Граничные условия в общем виде, когда оба цилиндра вращаются, будут

, , .

Начальное условие .

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

. (16)

Для решения использован метод Фурье и в результате получено решение в виде

, (17)

где ; ;

.

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

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

В связи с этим рассмотрена трехмерная осесимметричная относительно оси оz задача течения вязкой жидкости, т.е. . Дифференциальные уравнения движения Навье-Стокса и неразрывности представлены в конических координатах . При этом начало отсчета оси совпадает с осью z. Данная система координат ортогональна. Коэффициенты Ламе для нее соответственно равны: ; ; . Полученные уравнения движения и неразрывности с учетом сделанных допущений ввиду их громоздкости приведены в диссертации.

Уравнения приведены к безразмерному виду с помощью подстановок

; ;

; . (18)

Относительная координата определяется в виде

.

Используя (18) и проведя оценку значимости их членов с учетом численных величин параметров, реально осуществимых в роторных аппаратах, получена система уравнений движения жидкости в зазоре

, , .  (19)

Условие неразрывности удовлетворяется с большой точностью.

Граничные условия имеют вид: а) на наружной поверхности ротора при ()

, , ;

б) на внутренней поверхности статора при ()

, , .

Система уравнений (19) не имеет решения в замкнутой форме. Для решения функции H и G представлены разложенными в ряды Тейлора вблизи координаты . После использования граничных условий и преобразований получена система уравнений относительно , ,. Использовав численный метод решения, получаем результаты по значениям , и С а также графические зависимости меридиональной , азимутальной G и нормальной к поверхности конуса H, составляющих скорости от безразмерного зазора и нормальной к поверхности конуса составляющей скорости на входе и выходе из зазора W. Зависимость требуется при вычислении энергопотерь в радиальном зазоре.

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

Баланс энергопотерь имеет вид:

;

Рис.5 Зависимость нормальной H, меридиональной и окружной G составляющих скорости от

Мощность, затрачиваемая для сообщения кинетической энергии жидкости в роторе:

- в случае цилиндрического ротора

;  (20)

- в случае конического ротора

. (21)

Формулы (20) и (21) учитывают, что в реальных аппаратах . Силы внутреннего трения в радиальном зазоре и в осевом зазоре между торцом ротора и корпусом роторного аппарата определяются согласно закону Ньютона. Затем через момент трения находится затрачиваемая мощность в радиальных зазорах:

- в случае цилиндрического ротора

;  (22)

- в случае конического ротора

.  (23)

Величина определяется по предложенной методике.

Рис.6 Зависимость диссипируемой в зазоре мощности от и времени.

Трение в осевом зазоре определяется по выражению

.  (24)

Механические потери обычно не превышают (5Е10%) от суммы .

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

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

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

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

(25)

Переменное содержание свободного газа в жидкостной среде определяется по зависимости

.  (26)

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

. (27)

При

,  (28)

где ; .

Форма импульса отрицательного давления определяется для стадии закрывания канала статора из решения уравнений нестационарного течения: а) несжимаемой среды выражений (1)-(4), б) сжимаемой среды (13), (4), после взятия производной по времени.

Начальные условия для уравнения (25) имеют вид

, . (29)

Уравнения (25)-(28) с учетом (29) решены численным методом при ; , м/с; м/с, МПа; кг/м3; ; ; ; .

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

, , , где .

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

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

Определено изменение величины импульса кавитационного давления при постоянном газосодержании от акустического числа кавитации по выражению

.

Рис.7 Зависимость кавитационного импульса давления от : ; ; ; .

Рис.8 Зависимость от начального газосодержания : .

На всех полученных графиках имеется резко выраженный максимум (рис.7). Экспериментальное подтверждение изложено в 5 главе.

Основным результатом проведенного исследования является впервые полученная зависимостьот начального газосодержания (рис.8), которую можно аппроксимировать выражением

. (30)

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

Определены зависимости от начального газосодержания. Показано, что с увеличением с до 0,5 максимальная амплитуда кавитационного давления падает почти в 500 раз.

Выявлено, что с увеличением начального радиуса пузыря с 0,1 до 4, возрастает примерно в 300 раз.

Установлено, что с увеличением с 0,1 до 4, критерий возрастает примерно в 1,5 раза, а влияние критериев Вебера и Рейнольдса на незначительно и для упрощения расчетов эти критерии с достаточной точностью можно не учитывать т.к. основным фактором воздействия является акустическая кавитация,  соответствующая слагаемому с в (25).

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

Влияние симплексов и на критерий не выявлено. Увеличение критерия влияет на незначительно. Таким образом, можно утверждать, что уравнение (30) можно использовать при изменении конструктивных и режимных параметров в исследованных границах.

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

.  (31)

Если канал закрыт промежутком между каналами ротора, основная частота определяется выражением

.  (32)

Кроме того, в канале статора возможен нелинейный резонанс, определяемый формулами (31) и (32) с удвоением знаменателя. Для камеры озвучивания, в первом приближении, резонансные частоты определяются выражением (32), а также возможностью возникновения нелинейного резонанса. Условием совместного резонанса в канале статора и в камере озвучивания является равенство .

Основная частота колебаний, генерируемая модулятором роторного аппарата, определяется при

.  (33)

Частота резонансных колебаний определяется задающей системой - каналом или камерой озвучивания, независимо от механизма возбуждения (уравнение (33)). Очевидно, наилучший случай - равенство выражения (33) выражению (31) или (32), конкретный результат должен быть подтвержден экспериментально.

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

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

Условие возникновения резонансного процесса определяется соотношением

,  (34)

где время, когда канал статора закрыт, равно

,  (35)

а время пробега импульса давления по длине канала

.  (36)

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

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

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

Полученные осциллограммы подтверждают необходимость использования 2х моделей течения - для несжимаемой и сжимаемой жидкости. Граница применимости моделей отвечает условию (14). Осциллограмма на рис.9 качественно соответствует аналитическим зависимостям, полученным в 3 главе (рис.4).

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

.  (37)

Критерий кавитации Стайлса удобен для практического использования, т.к. его легко регулировать, изменяя давление в камере озвучивания и перепад

а)	б)
в)	Рис.9 Осциллограммы пульсационного давления в камере озвучивания роторного аппарата при а=0,05;  zp=zc=20; =0,1; l/ap=5 и значениях Kк: а) 0,3 б) 1,1 в) 1,3.

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

.

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

Подтверждено возникновение в длинноканальном роторном аппарате явления резонанса, условие его существования в виде (34)-(36) и влияние на величину перепада давлений между ротором и статором (рис.10)

Рис.10 Зависимость от угловой частоты вращения ротора: а=0,01, =0,1, l/ap=2, Q=0,0041 м3/с.

Рис.11 Зависимость отрицательного импульса давления от критерия Kк: 1-=2,75105Па; 2-=2,75105Па; 3-=2,75105Па; а=0,05, =0,1, l/ap=5.

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

Рис.12 Зависимость кавитационнного импульса давления от статиченского давления: а=0,05; l/aр = 5.

Рис.13 Зависимость коэффициента скорости растворения от угловой скорости вращения ротора: а=0,05; l/aр = 2.

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

Показано, что согласно условию (14) при 0,4<<0,8 для описания течения служит модель течения несжимаемой жидкости, что соответствует рис. 2а (левая часть графика). Для значений >0,8 используется модель течения с учетом сжимаемости среды, что подтверждает теоретический график на рис.3.

Экспериментально подтверждена теоретическая зависимость кавитационного импульса давления от числа акустической кавитации (рис.7). Увеличение статического давления в камере при постоянном возбуждающем импульсе отрицательного давления равнозначно увеличению числа акустической кавитации (рис.12).

Подтверждено влияние массовых сил (величины критерия ) на интенсивность кавитации. График изменения аналогичен графику на рис.11, за исключением участка 0,76<<0,8, где кавитация отсутствует, т.к. для ее возбуждения недостаточна величина .

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

. (38)

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

. (39)

Необходимое условие работы роторного аппарата в кавитационном режиме соответствует неравенству

. (40)

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

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

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

В заводских условиях проведены испытания по выяснению влияния качества полученной СОЖ на стойкость режущего инструмента изготовленного из P18. Показано, что стойкость резцов, оцениваемая по износу их задней поверхности, гарантированно возрастает минимум на 30% без добавления дополнительных присадок. Показано, что использование высокодисперсной СОЖ позволило повысить чистоту обработанной поверхности на один - два класса.

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

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

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

В методике использованы экспериментальные зависимости (36), (38), (39) и условие (40), характеризующие работу роторного аппарата в кавитационном режиме. Для возникновения резонансного режима необходимо использовать выражения (31), (33)-(36). По зависимости (30) определяется критерий акустической кавитации. При выполнении условия (14) для определения модуля отрицательного ускорения используется модель течения для сжимаемой жидкости - выражения (4) и (13), при не выполнении условия (14) используется модель течения для несжимаемой жидкости - выражения (4), (1), (2). Сравниваются  выражение для и полученное, с использованием критерия Стайлса. В случае несовпадения производят пересчет, изменяя параметры в заданных границах. В заключение определяется потребляемая мощность: для цилиндрических ротора и статора - по зависимостям (20), (21), (24); для конических - по выражениям (21), (23), (24). Расчеты проводятся для каждого значения угловой скорости вращения ротора. Учитывается, что режим наиболее развитой кавитации соответствует максимуму.

Разработана методика оптимизационного расчета аппарата. Выбрана модификация метода градиента с постоянным шагом.

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

Целевая функция имеет вид

.

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

Используются соотношения математической модели течения сжимаемой среды (4), (13-14) соотношения модели радиально-сферических колебаний кавитационного пузыря (25-30), (37-40) соотношения условий возникновения резонанса (32-36).

Разработана программа на языке С++ Builder 6.0. Программа позволяет получить основные конструктивные и режимные параметры: ; ; ; ;и(каналы наклонные); ;;;;;; ;; ; М; ;;N.

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

Из условий возникновения резонансного режима и стоячей волны разработана конструкция роторного аппарата (SU 1187858).

Для достижения наиболее интенсивной кавитации изменением объемной концентрации содержания свободного газа в жидкости, влияющая на , запатентована конструкция (RU 2287360).

Всего запатентовано 16 конструкций роторных аппаратов.

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

Роторные аппараты нашли эффективное применение в процессах получения СОЖ. Реализованная на заводе Тамбовполимермаш конструкция защищена авторским свидетельством SU 1187858. Установки для приготовления СОЖ были использованы в трех цехах завода. За время эксплуатации роторные аппараты работали надежно, отказов в работе из-за конструктивных недостатков не выявлено. Стойкость режущего инструмента (резцы, сверла, фрезы) возросла не менее чем на 30%. Время расслаивания возросло почти в 8 раз. На Мичуринском заводе Прогресс изготовлены малой серией (5 шт.) роторные аппараты, защищенные а.с. SU 1389830. В получаемой СОЖ 80-90% частиц эмульсола имели размер 1Е2 мкм. Время расслаивания выросло в 10 раз, и стойкость режущего инструмента из Р18 возросла в среднем на 30%. После проведения успешных испытаний на заводе в корпусе, объединяющем 4 цеха, создан участок для централизованной раздачи СОЖ к металлообрабатывающим станкам. В результате были улучшены условия труда, снизилось количество ручного труда. Примерно такие же данные по повышению стойкости режущего инструмента и времени расслаивания получены при эксплуатации установок на базе роторного аппарата (RU 2225250, RU 2230616) на Борисоглебском заводе ОАО Патроны, ООО Грибановсий машиностроительный завод и ОАО РЖД Локомотивное ДЕПО г. Тамбов (RU 2287360).

На ОАО Котовский лакокрасочный завод внедрен роторный аппарат для диспергирования двуокиси титана (марка Р-02), применяемого для производства эмали ПФ-115 белая. Результаты показали следующее: снизились время приготовления продукта и удельный расход электроэнергии.

На НПФ Лионик г. Москва внедрение роторного аппарата при проведении совмещенных процессов смешения-диспергирования-гомогенизации при производстве сухих концентратов натуральных напитков позволило в два раза увеличить производительность при улучшении качества готовой продукции.

Основные результаты и выводы.

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

- выполнены с использованием зонного подхода теоретические исследования нестационарного течения среды в каналах роторного аппарата с учетом сжимаемости в условиях гидравлического удара, подтвержденные экспериментальными данными, позволяющие расширить область полученных режимных и конструктивных параметров роторных аппаратов и повысить точность расчетов; на основании поученных данных установлено, что для эффективной работы аппарата значение критерия 10, симплекса 0,3, величина должна приниматься минимально возможной, с учетом прочности ротора;

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

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

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

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

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

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

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

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

- результаты исследований, предложенные методики расчета и разработанные на их основе конструкции роторных аппаратов использованы на Мичуринском ОАО Прогресс, ОАО Тамбовполимермаш для приготовления эффективной смазочно-охлаждающей жидкости с реальным годовым эффектом 125000 рублей (в ценах до 1990 г.), а также на Борисоглебском заводе ОАО Патроны, ООО Грибановский машиностроительный завод, ОАО Тамбовский завод Комсомолец им. Н.С. Артемова, ОАО РЖД Локомотивное ДЕПО г. Тамбов, на Котовском заводе ОАО КЛКЗ при производстве эмали белой ПФ-115 с экономическим эффектом 410 руб. на 1 тонну эмали при годовом производстве до 3000 т. (в ценах 2005 г.). На НПФ Лионик г. Москва при использовании роторного аппарата для проведения совмещенных процессов смешивания - диспергирования - гомогенизации при производстве сухих концентратов натуральных напитков получен в 2005-2006 гг. реальный экономический эффект 16,0 млн. руб.

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

Основные обозначения

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

Индексы: - безразмерная величина;- первая и вторая производная; - наибольшая величина.

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

1. Червяков, В.М. Гидродинамические и кавитационные явления в роторных аппаратах : монография/ В.М. Червяков, В.Ф. Юдаев. - М.: Машиностроение -1, 2007. - 128 с.
2. Червяков, В.М. Течение ньютоновской жидкости в зазоре между коническими проницаемыми поверхностями/ В.М, Червяков, А.А. Коптев// Инженерно-физический журнал. - 2006. - Т.79, №2. - С. 92-98.
3. Исследование процесса растворения серы в смеси масел в роторном аппарате/ В.М. Червяков [и др.]// Изв.ВУЗ. Химия и химическая технология. - 2006. - №4 - С.95-97.
4. Червяков, В.М. Нестационарное течение идеальной сжимаемой среды в каналах роторного аппарата/ В.М. Червяков, Ю.В. Воробьев// Теоретические основы химической технологии. - 2005. - Т.39, №1. - С. 65-71.
5. Червяков, В.М. Кавитационные явления в газожидкостной среде/ В.М. Червяков, В.Ф. Юдаев// Проблемы машиностроения и автоматизации. - 2004. - №4. - С.73-77.
6. Червяков, В.М. Определение энергозатрат в роторных аппаратах/ В.М. Червяков, А.А. Коптев// Химическое и нефтегазовое машиностроение. - 2005.- №4. - С.10-12.
7. Червяков, В.М. Обобщенная методика расчета роторного аппарата с учетом импульсной акустической кавитации/ В.М. Червяков, Ю.В. Воробьев, В.Ф. Юдаев// Вестник ТГТУ. - 2005. - Т.11, №3. - С. 683-689.
8. Режимы работы технологического оборудования с возбуждением кавитации/ В.М. Червяков [и др.]//Вестник ТГТУ. - 2005. - Т.11, №3. - С. 399-402.
9. Червяков, В.М. Нестационарное течение жидкости в зазоре между ротором и статором/ В.М. Червяков, В.И. Галаев, А.А. Коптев// Вестник ТГТУ. - 2003. - Т.9, №4. - С. 646-652.
10. Червяков, В.М. Подобие процессов нестационарного течения жидкости в модуляторе роторного аппарата/ В.М. Червяков, Ю.В. Воробьев// Вестник ТГТУ. - 2002. - Т.8, №4. - С.618-622.
11. Промтов, М.А. Кинетика растворения NaCl в воде при обработке в роторно-импульсном аппарате/ М.А. Промтов, В.М. Червяков// ИзвВУЗ. Химия и химическая технология. - 2000. - №6. - С.133-135.
12. Червяков, В.М. Нестационарное течение жидкости во вращающихся каналах роторного аппарата/ В.М. Червяков, В.И. Галаев, А.А. Коптев// Вестник ТГТУ. - 2000. - Т.6, №4. - С.611-616.
13. Воробьев, Ю.В. Исследование взаимного влияния процессов при работе роторного аппарата с модуляцией потока и вспомогательного оборудования/ Ю.В. Воробьев, М.А. Промтов, В.М. Червяков// Вестник ТГТУ. - 1996. - Т.2, №3. - С. 266-270.
14. Промтов, М.А. О формировании кавитационных образований в роторном аппарате с модуляцией потока/ М.А. Промтов, В.М. Червяков// Вестник ТГТУ. - 1995. - Т.1, №3-4. - С. 311-315.
15. Chervyakov, V.M. Determination of Power Consumption in Rotary Equipment/ V.M. Chervyakov, A.A. Koptev// Chemical and Petroleum Engineering. - 2005. - V.41, № 3 - 4. - P.180 - 184.
16. Chervyakov, V.M. Flow of a Newtonian fluid in the gap between conical permeable surfaces// V.M. Chervyakov, A.A. Koptev// Journal of engineering Physics and Thermophysics. - 2006. - V.79, №2. - P.301 - 308.
17. Червяков, В.М. Математическая модель течения жидкости во вращающихся каналах прямоугольного сечения роторного аппарата/ В.М. Червяков, В.И. Галаев, А.А. Коптев// Математическое моделирование физических, экономических, технических, социальных систем и процессов: тр. 4 междунар. конф. - Ульяновск: УГУ, 2001. - С. 158-160.
18. Нестационарное течение сжимаемой жидкости во вращающихся каналах роторного аппарата/ В.М. Червяков [и др.]// Математические методы в технике и технологиях. ММТТ-15: тр. 15 междунар. науч. конф. - Тамбов: ТГТУ, 2002. - Т.10. - С. 57-58.
19. Неустойчивость течения в каналах роторного аппарата/ В.М. Червяков [и др.]// Математические методы в технике и технологиях. ММТТ-15: тр. 15 междунар. науч. конф. - Тамбов: ТГТУ, 2002. - Т.10. - С. 58-59.
20. Червяков, В.М. Использование роторного аппарата для приготовления высококачественной смазочно-охлаждающей жидкости/ В.М. Червяков, Ю.В. Воробьев, Ю.В. Родионов// Фундаментальные и прикладные проблемы технологии машиностроения. Технология -2002: материалы Междунар. науч.-техн. интернет конф. - Орел: ОреГТУ, 2002. - С. 261-264.
21. Применение роторного аппарата для ультразвуковой гидроабразивной обработки пересекающихся отверстий/ В.М. Червяков [и др.] фундаментальные и прикладные проблемы технологии машиностроения. Технология-2002: материалы Междунар. науч. -техн. интернет конф. - Орел: ОреГТУ, 2002. - С. 265-267.
22. Математическая модель нестационарного течения слабосжимаемой жидкости в модуляторе роторного аппарата/ В.М. Червяков [и др.]// Математическое моделирование физических, экономических, технических, социальных систем и процессов: тр. 5 Междунар. конф. - Ульяновск: УГУ, 2003. - С. 205-207.
23. Червяков, В.М. Разгонное течение вязкой среды в зазоре между ротором и статором роторного аппарата/ В.М. Червяков, А.В. Кожуховский// Актуальные проблемы современной науки. Ч. 1-3 Математика. Механика. Машиностроение: тр. 4 Междунар. конф. - Самара: СамГТУ, 2003. - С. 120.
24. Использование роторного аппарата в ресурсосберегающих технологиях/ В.М. Червяков [и др.]// Наука и образование: материалы 5 Междунар. конф./ Кемеровский гос. ун-т. Беловский институт(филиал). - Белово: Беловский полиграфист, 2004. - Ч.1. - С.456-459.
25. Червяков, В.М. Математическая модель нестационарного течения несжимаемой жидкости в модуляторе роторного аппарата/ В.М. Червяков, А.В. Кожуховский, Л.В. Чичева-Филатова// Математическое моделирование и краевые задачи: тр. Всерос. науч. конф. - Самара: СамГТУ, 2004. - Ч.2.- С. 263-265.
26. Червяков, В.М. Определение энергетических затрат при работе роторного аппарата в энергосберегающих технологиях/ В.М. Червяков, Ю.М. Радько, А.И. Четырин// Энерго- и ресурсосбережение - XXI век: сб. материалов 2 Междунар. науч.-практ. интернет конф. - Орел: ОреГТУ, 2004. - С. 188-190.
27. Повышение энергетической эффективности роторного аппарата/ В.М. Червяков [и др.]// Теплофизические измерения при контроле и управлении качеством: материалы 5 Междунар. теплофиз. шк. - Тамбов: ТГТУ, 2004. - С. 290-291.
28. Червяков, В.М. Математическая модель течения среды в зазоре между коническими ротором и статором/ В.М. Червяков, Ю.В. Воробьев, А.И. Четырин// Математические методы в технике и технологиях ММТТ-18: сб. тр. 18 Междунар. науч. конф. - Казань: КГТУ, 2005. - Т.5. - С. 134-136.
29. Определение мощности роторных аппаратов/ В.М. Червяков [и др.]//Глобальный научный потенциал: сб. материалов Междунар. науч.-практ. конф. - Тамбов: Першина, 2005. - С. 60-61.
30. Перспективные направления в конструировании роторных аппаратов для снижения удельных энергозатрат/ В.М. Червяков [и др.]// Энерго- и ресурсосбережение - XXI век: сб. материалов 3 Междунар. науч.-практ. интернет конф. - Орел: ОреГТУ, 2005. - С. 226-227.
31. Обобщенная модель течения среды в модуляторе роторного аппарата/ В.М. Червяков [и др.]// Математическое моделирование физических, экономических, социальных систем и процессов: тр. Шестой Междунар. науч. конф. - Ульяновск: УГТУ, 2005. - С.136-137.
32. Червяков, В.М. Закономерности энергопотребления роторным аппаратом/ В.М. Червяков, А.И. Четырин// Методы и алгоритмы прикладной математики в технике, медицине и экономике: материалы 6 Междунар. науч.-практ. конф. - Новочеркасск: ЮРГТУ(НПИ), 2006. - Ч.1. - С.41.
33. Червяков, В.М. Влияние массовых сил на течение среды в каналах роторного аппарата/ В.М. Червяков, А.И. Четырин// Моделирование, теория, методы и средства: материалы 6 Междунар. науч. - практ. конф. - Новочеркасск : ЮРГТУ(НПИ), 2006. - Ч.4. - С.61-63.

Новизна предложенных технических решений подтверждена авторскими свидетельствами и патентами на изобретение:

34. А. с. 1187858 SU B01F 7128. Роторный аппарат/ В.В. Белик, В.А. Колдин, М.М. Свиридов, В.М. Червяков, Е.С. Шитиков; Тамб. ин-т. хим. машиностроения. - №3685908/23-26; заявл. 15.11.83; опубл. 30.10.85, Бюл. №40.
35. А. с. 1389830 SU B01F 7/28. Роторный аппарат/ В.М. Червяков, Ю.В. Воробьев, В.И. Токарев, А.Г. Ткачев, В.И. Жеребятьев; Тамб. ин-т. хим. машиностроения. - №4073350/31-26; заявл. 02.06.86; опубл. 23.04.88, Бюл. №15.
36. А. с. 1428402 SU B01D 3/30. Роторный аппарат/ В.М. Червяков, Ю.В. Воробьев, Б.С. Спиридонов, А.Г. Ткачев, Э.И. Приходько; Тамб. ин-т. хим. машиностроения. - №4072175/31-26; заявл.28.05.86; опубл. 07.10.88, Бюл. №37.
37. А. с. 1558654 SU B24C 3/06, B24B 1/04. Способ гидроабразивной обработки сквозных каналов/ В.М. Червяков, М.А. Промтов, Ю.В. Воробьев, Н.Н. Чернецкий, А.Г. Ткачев, А.В. Суворов; Тамб. ин-т. хим. машиностроения. - №4366902/25-08; заявл. 15.01.88; опубл. 23.04.90, Бюл. №15.
38. А. с. 1674942 SU B01F 7/28. Роторный аппарат/ В.М. Червяков, М.А. Промтов, Ю.В. Воробьев, А.Г. Ткачев; Тамб. ин-т. хим. машиностроения. - №4440434/26; заявл. 03.06.88; опубл. 07.09.91, Бюл. №33.
39. А. с. 1719045 SU B01F 7/28. Роторный аппарат/ В.М. Червяков, Е.С. Шитиков, Ю.В. Воробьев, М.А. Промтов, В.А. Колдин; Тамб. ин-т. хим. машиностроения. - №4770235/26; заявл.18.12.89; опубл. 15.03.92, Бюл. №10.
40. А. с. 1722802 SU B24C 3/06. Способ гидроабразивной обработки сквозных каналов/ В.М. Червяков, М.А. Промтов, Ю.В. Воробьев, В.Н. Поляков, В.И. Жеребятьев, Н.Н. Чернецкий; Тамб. ин-т. хим. машиностроения. - №4771218/08; заявл.18.12.89; опубл. 30.03.92, Бюл. №12.
41. А. с. 1768269 SU B01F 7/28. Роторный аппарат/ М.А. Промтов, В.М. Червяков, Ю.В. Воробьев; Тамб. ин-т. хим. машиностроения. - №4786904/26; заявл.30.01.90; опубл. 15.10.92, Бюл. №38.
42. А. с. 1773469 SU B01F 7/28. Роторный аппарат/ М.А. Промтов, В.М. Червяков, Ю.В. Воробьев, В.И. Тебекин; Тамб. ин-т. хим. машиностроения. - №4820501/26; заявл.26.04.90; опубл. 07.11.92, Бюл. №41.
43. Пат. 2150318 RU B01F 7/00. Роторный аппарат/ А.А. Коптев, В.М. Червяков, М.А. Промтов; Тамб. гос. техн. ун-т. - №98120226/12; заявл.10.11.98; опубл. 10.06.2000, Бюл. №16.
44. Пат. 2155634 RU B01F 7/00. Роторный аппарат/ М.А. Промтов, В.М. Червяков, Ю.В. Воробьев, М.В. Монастырский; Тамб. гос. техн. ун-т. - №98116768/12; заявл.08.09.98; опубл. 10.09.2000, Бюл. №25.
45. Пат. 2165292 RU B01F 7/00. Роторный аппарат/ В.М. Червяков, М.А. Промтов, А.А. Коптев; Тамб. гос. техн. ун-т. - №99119141/12; заявл.06.09.1999; опубл. 20.04.2001, Бюл. №11.
46. Пат 2225250 RU B01F 7/28, 7/00, 3/08. Роторный аппарат/ В.М. Червяков, Ю.В. Воробьев, В.Ф. Юдаев, Е.С. Шитиков; Тамб. гос. техн. ун-т. - №2002107488/15; заявл.25.03.2002; опубл. 10.03.2004, Бюл. №7.
47. Пат. 2230616 RU B06B 1/20. Роторный аппарат/ В.М. Червяков, Е.С. Шитиков, А.А. Коптев, В.И. Галаев; Тамб. гос. техн. ун-т. - №2002107487/28; заявл.25.03.2002; опубл. 20.06.2004, Бюл. №17.
48. Пат. 2294236 RU B01F 7/26 Роторный аппарат/ В.М. Червяков, Ю.В. Воробьев, А.А. Коптев, В.Ф. Юдаев, Ю.В. Родионов, Л.В. Чичева-Филатова, В.А. Алексеев; Тамб. гос. техн. ун-т. - №2004133695/15; заявл.18.11.2004; опубл.27.02.2007, Бюл. №6.
49. Пат 2287360 RU B01F 7/00. Устройство для физико-химической обработки жидкой среды/ В.М. Червяков, В.Ф. Юдаев, В.И. Биглер, Л.В. Чичева-Филатова, В.А. Алексеев, Н.И. Акулов; Тамб. гос. техн. ун-т. - №2004133696/15; заявл.18.11.2004; опубл. 27.04.2006, Бюл. №32.
50. Свидетельство об официальной регистрации программы для ЭВМ. - №2005610721. Программное обеспечение для определения параметров, характеризующих режим наиболее интенсивной кавитации в газожидкостной среде/ В.М. Червяков, В.Ф. Юдаев, О.А. Ступников, А.И. Четырин; Тамб. гос. техн. ун-т. - №2005610163; заявл. 31.01.2005. Зарегистрировано 25.03.2005.

Авторефераты по всем темам  >>  Авторефераты по техническим специальностям