WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!


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

Белоусов Александр Сергеевич

ГИДРОДИНАМИКА ПРОЦЕССОВ С НЕОДНОРОДНЫМИ СТРУКТУРАМИ ЗАКРУЧЕННЫХ ГЕТЕРОГЕННЫХ ПОТОКОВ В ВИХРЕВЫХ АППАРАТАХ

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

АВТОРЕФЕРАТ

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

Москва 2010

Работа выполнена в Государственном образовательном учреждении высшего профессионального образования «Московский государственный текстильный университет им. А.Н. Косыгина» на кафедре «Процессы, аппараты химической технологии и безопасность жизнедеятельности»

Научный консультант: доктор технических наук, профессор Сажин Борис Степанович

Официальные оппоненты: доктор технических наук, профессор Авдюнин Евгений Геннадиевич доктор технических наук, профессор Кабанов Александр Николаевич доктор технических наук, профессор Горшенин Павел Александрович

Ведущая организация: Государственное образовательное учреждение высшего профессионального образования «Российский заочный институт текстильной и легкой промышленности» (ГОУ ВПО «РосЗИТЛП»)

Защита диссертации состоится 28 декабря 2010 г. в 10 часов на заседании заседании диссертационного совета Д 212.139.03 при Московском государственном текстильном университете им. А.Н. Косыгина по адресу:

119071, г. Москва, ул. Малая Калужская, д.

С диссертацией можно ознакомиться в библиотеке Государственного образовательного учреждения высшего профессионального образования «Московский государственный текстильный университет имени А.Н.Косыгина».

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

Ученый секретарь диссертационного совета, доктор технических наук А.В. Фирсов

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



Актуальность темы.

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

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

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

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

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

Диссертационная работа выполнялась в соответствии с планами Проблемного совета РАН по теоретическим основам химической технологии (ТОХТ) 2003-2008 гг.; планами важнейших НИР ГКНТ АН СССР по ТОХТ на 1986-1990 гг.; хоздоговорными и госбюджетными планами МТИ им. А.Н. Косыгина 1981-1991 гг.; планами НИР МГТУ им. А.Н. Косыгина 2003-2008 гг.

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

Основные задачи

исследования:

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

– Разработана группа многопоточных моделей неоднородных структур и методы идентификации распределения времени пребывания (РВП) твердой фазы в вихревых аппаратах. Методика идентификации неоднородных структур РВП, основана на нелинейной оптимизации совместно с анализом приближенных функций интенсивности.

– Выявлены закономерности РВП для аппаратов с неоднородными потоками: найдено, что для режима винтового течения в вихревых аппаратах адекватна ячеечная модель с числом секций N = 8 50 ; для циклонов при скоростях * больше критической V0 установлено существование эффекта распада РВП на две ветви; для многофункциональных аппаратов ВЗП-ВК разделение РВП на две ветви происходит при определенной высоте перелива; установлено существование двух циркуляционных контуров в аппарате фонтанирующего слоя; выявлено изменение структуры перемешивания в вихревой камере.

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

Практическая значимость и реализация результатов работы.

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

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

– Предложенные методы идентификации позволяют с высокой достоверностью определять фракционные характеристики и параметры РВП; методика и предложенные модели позволяют эффективно выявлять и описывать гидродинамические неоднородности потоков для широкого класса аппаратов.

– Разработан и исследован новый вид вихревых пылеуловителей типа ВП, которые имеют номинальное соотношение входных расходов 1:1, широкий рабочий диапазон по соотношению расходов, малочувствительны к колебаниям нагрузок; разработаны типоразмеры аппаратов ВП диаметром от 0,2 до 1,2 м и установки ВП-700х2 и ВП-800х3, которые обеспечивают очистку газов в диапазоне от 500 до 22000 м3 /час.

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

– По результатам выполненных исследований разработаны и внедрены в различных производствах аппарат ВП-300, пылеуловители ВП-700, установка ВП-700х2, комбинированная установка ВП-800х3. Внедрены пылеуловители типа ВЗП в составе установок ВЗП-600х3, ВП-1, ВП-2, ВП-3. Результаты работы были использованы при разработке и внедрении и эксплуатации пылеуловителей ВЗП-2000, ВЗП-2000У, ВЗП–ЦРП, при проектировании сушилки с одновременным улавливанием СВЗП-600, микрогранулятора ВЗП-350-МГ.

– Результаты работы используются в учебных курсах МГТУ им. А.Н.

Косыгина: «Моделирование химико-технологических процессов», «Оптимизация химико-технологических процессов», «Прикладная математика и кинетика процессов переноса».

Апробация работы. Результаты работы докладывались на семинарах Совета РАН по технологическим процессам с твёрдой фазой (2005г.), по теоретическим основам химической технологии (1987г.); на международном совещании РАН, посвящённом столетию со дня рождения академика Н.М. Жаворонкова (2007г.); at The Second Nordic Drying Conference (NDC–03), Copenhagen, Denmark, 2003, а также более чем на 30 международных, всероссийских и всесоюзных конференциях, в том числе на: международных конференциях «Экологические и ресурсосберегающие технологии промышленного производства» (Витебск, 2006г.); «Современные наукоемкие технологии и перспективные материалы текстильной и легкой промышленности» (Иваново, 2004, 2006, 2007гг.); «Современные технологии и оборудование текстильной промышленности» (МГТУ им. А. Н. Косыгина, 2003–2008гг.); на международных конгрессах по химии и химической технологии «МКХТ» (РХТУ им. Д.И. Менделеева 2003–2008гг).; всесоюзных конференциях «Очистка газовых выбросов промышленных предприятий» (Тольятти, 1990г.); «Технология сыпучих материалов» (Ярославль, 1989г.); «Новые процессы, оборудование и гибкие производственные системы для многономенклатурных химических производств» (Черкассы, 1989г.); «Современные машины и аппараты химических производств» (Чимкент, 1988г.); «Аэрозоли и их применение в народном хозяйстве» (Юрмала, 1987г.); «Создание и внедрение современных аппаратов с активными гидродинамическими режимами для текстильной промышленности и произ водства химических волокон» (Москва, 1981,1983,1985,1989 г.); «Повышение эффективности тепломассообменных и гидродинамических процессов в текстильной промышленности и производстве химических волокон» (Москва, 1985,1988 г.); «Методы кибернетики в химии и химической технологии» (Грозный, 1984г.); на «Всесоюзном совещании по математическому моделированию и управлению высокотемпературными процессами в циклонных и вихревых аппаратах» (Одесса, 1980г.); научных конференциях МТИ им. А. Н.

Косыгина, 1982–1987гг. и др.

Публикации.

По теме диссертации опубликовано 105 научных работ, из них 50 работ приведены в списке автореферата, 18 статей в журналах из перечня ВАК РФ.

Структура и объем диссертации.

Диссертация состоит из введения, 8 глав, приложений, основных выводов и списка литературы. Общий объем диссертации 372 страниц, таблиц – 14, рисунков – 93, библиография включает 315 наименований.

СОДЕРЖАНИЕ РАБОТЫ

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

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

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

Вместе с тем основные прикладные методы моделирования гидродинамики вихревых аппаратов исходят из однородной и фиксированной структуры потоков, и, следовательно, не могут прогнозировать возможные неоднородности. Другим направлением является применение фундаментальных методов вычислительной гидродинамики (CFD– Computational Fluid Dynamics), основанных на так называемых стандартных (базовых) полуэмпирических моделях турбулентности. Однако из анализа применения наиболее популярной группы двухпараметрических базовых моделей (2P) следует, что они в основном работоспособны только при малых степенях крутки потока.

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

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

На основе анализа сформулированы основные задачи исследования.

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

Выполнен анализ системы физических эффектов турбулентных потоков в вихревых технологических аппаратах. В соответствие с результатами М.А.

Гольдштика применен аналог вариационного принципа о минимизации диссипации кинетической энергии. Предложена вихревая концепция: приосевая генерация турбулентности, пропорциональная крутке потока и частичное подавления турбулентности в основном ядре течения. В этом случае универсальной постоянной задачи является турбулентное число Рейнольдса RT = Q /(2 µt H ), где Q – расход воздуха через аппарат, H – высота аппарата, µt – турбулентная вязкость. Исходя из аналога вариационного принципа получено RT = kt , где – приведенный параметр закрутки, kt 26.

Турбулентное течение газа в аппарате описано уравнениями Рейнольдса в цилиндрической системе координат ( z,r, ). Для замыкания уравнений Рейнольдса применена модель анизотропной турбулентной вязкости.

Предложена гидродинамическая схема течения в технологических вихревых аппаратах, состоящая из ядра потока, обладающего свойствами свободной турбулентности, и тонкого пограничного слоя. Таким образом, для условий течения в аппарате получена векторная турбулентная вязкость µt = {µ0, µr}, где µr - турбулентная вязкость, соответствующая напряжению r, а µ0 - турбулентная вязкость остальных компонент тензора турбулентных напряжений. Для среды с постоянной плотностью и постоянными коэффициентами µ, µ0 и µr получены следующие уравнения Г Г 2Г Г - = (µ - µr )r - + r(µ - µ0) Г, (1-а) r z z r r2 r z 1 1 1 - + = r z r r dz r r2 z, (1-б) 2 2 1 2Г Г = (µ - µ0) + + - + z2 r2 r r r2 r3 z 2 2 1 + - = -r, (1-в) z2 r2 r r где – функция тока, – циркуляция, – азимутальная компонента вихря.





По аналогии с корреляцией П. Брэдшоу принято µr = µ0 /, где S S - интегральный параметр анизотропии, зависящий от эффективной закрутки.

Используя результаты исследований по струйным и отрывным течениям для турбулентной вязкости в ядре потока получено 0, 2 (1- K)3(1+ ФM ) (1+ ФM ) Q K1 µ0 = +, (2) 2 kT D H 2 2 1- RT H R2 - R1 где ФMi, i = 1,2 - модельные степени закрутки во входных сечениях первого (нижнего приосевого) и второго (верхнего периферийного) потоков, D – диаметр аппарата; K - соотношение расходов потоков: K = Q2 (Q1 + Q2). Для интегрального параметра анизотропии предложено уравнение S 0, ={ DT [(1- K)ФM1 (2R2)+ K ФM2 (2Rвх)]} 0,32 +1, (3) S где Rвх - средний относительный радиус входа потока в тангенциальном завихрителе, R1– радиус втулки, R2 – радиус нижнего завихрителя.

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

В качестве универсального критерия, характеризующего влияние вращения на гидродинамические процессы в аппарате, принят интегральный параметр крутки S, определяемый через тангенциальнуюV и осевую Vz скорости R R S = Vr2dr R )rdr. (4) z z V (V 0 Краевые условия на входах различных типов завихрителей определяются через эквивалентные параметры крутки, которые рассчитаны по известным экспериментам для закрученных течений в трубах. Расчет модельной степени крутки ФMi, соответствующей эквивалентной, проводился по уравнению (4).

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

M = Mi /(1 + kMiX / R), (5) где X – координата вдоль оси в зоне однородного течения; Мi – момент импульса в начальном сечении; k – коэффициент трения.

Для лопаточного завихрителя при допущениях о равномерном распределении осевых скоростей и задании профиля V в виде V rn =const получаем:

Sc = 2tg0 N, (6) где: Sc – расчетный параметр крутки; 0 – угол закрутки при r=R; R - радиус 2 3-n канала; для n=3: N = ln(1 d0) (1- d ) ; для n3: N = (1- d0 ) [(3 - n)(1- d )];

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

Sc = (R2 - RT ) /(Fin) Rin / R sin , (7) где Fin – площадь входа, Rin – радиус центра тяжести площади входа до оси аппарата, - угол между осями входного патрубка и аппарата.

Для определения действительных значений S0 на срезе завихрителей, данные, полученные в литературе для X D 3 экстраполировались с помощью (5) до кромок завихрителей.

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

K(Sc ) = 1.28 Sc0,25 exp(-0.075 Sc ) Sc, (8) где Sc рассчитывается по уравнению (7). В этом случае погрешность составляет не более 8 %. Применяя уравнения (7) и (8) для циклонов ЦН-24, ЦН-15, ЦН-11 и СКЦН-34 получены значения крутки, хорошо коррелирующие с параметрами эффективности указанных циклонов.

Для расчета течения в пристеночной области предложена схема тонкого пограничного слоя в котором местные коэффициенты трения зависят от интегральных свойств потока. Для касательных напряжений на стенке i получено i = 0,158(1+ Ki 3 ФA V Vi 4 Re, i =Z, , (9) где V - характерная скорость, Vi - компонента скорости на стенке. Значения параметров в (9) приняты: KZ=10,2; = KZ / K = 4,15. Проектируя напряжения, действующие со стороны ядра потока на границе - слоя на плоскость, касательную к стенке, получены необходимые условия сопряжения.

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

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

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

а) б) в) Рис.1 Структура потоков газовой фазы в аппарате.

H =2,15; S1=0,6; S2 =5,2; а) K =0,8; б) K =0,62; в) K =0,39.

При K =0,39 структура движения газа наиболее сложная (рис. 1, в). Из представленных результатов следует, что в приосевой области возникает замкнутая циркуляционная зона III, содержащая сильное возвратное течение на оси (явление пузыревидного распада вихря – ПРВ). Выявлены особенности распада вихря применительно к аппаратам ВЗП: крупномасштабная неоднородность по всему сечению аппарата; устойчивая форма; хорошая управляемость формой распада при изменении режимно – конструктивных параметров. Таким образом установленное явление ПРВ можно классифицировать как специальную управляемую форму распада в отличие от классического неуправляемого распада вихря, возникающего в авиации при срыве потока с дельтовидных крыльев.

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

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

Между точками 1 и 2 осевая скорость в пристеночной зоне направлена вверх, что дает удержание в этой зоне материала.

Для экспериментального исследования гидродинамики полей течений был разработан стенд, на котором пневмометрическим методом исследовались поля скоростей для двух лабораторных моделей аппаратов, в шести сечениях, при трех соотношениях расходов. Сравнение результатов расчетов и экспериментов по полям скоростей при числе Рейнольдса Re0 =15000 показано на рис.2. Как видно из приведенных данных экспериментальные измерения достаточно хорошо подтверждают качественные и количественные результаты расчетов полей скоростей. Подтверждено существование неоднородных структур потоков, зоны распада вихря, значительных градиентов тангенциальной и осевой скоростей. Экспериментально полученная ширина зоны распада даже больше расчетной. Расчетные и экспериментальные поля турбулентного смешения также хорошо согласуются между собой.

Рис.2 Расчетные и экспериментальные значения скорости в сечениях: Z = Z / D.

Vz – осевая скорость: – расчет, – эксперимент; – тангенциальная скорость:

V – расчет, – эксперимент. а) Z =1,25, К=0,8; при К=0,39:, – расчет и эксперимент для ;, – соответственно расчет и эксперимент для ; б) Z =0,7, Vz V К=0,39; в) Z =1,65, К=0,62; г) Z =1,65, К=0,62.

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

Оценка смещения была проведена путем обработки литературных данных для радиальной скорости в циклоне, полученных с помощью специальной измерительной системы. Установлено, что для циклонов смещение оси по высоте близко к спиралевидному распаду вихря. Предложено исследовать оценку радиального стока по осесимметричной модели течения, а корректировку самой модели проводить по результатам замеров осевых и тангенциальных скоростей Из расчетов получено, что радиальный сток существенно изменяется по высоте аппарата и зависит от режима его работы. В основной части аппарата радиальный сток мал и близок к приближенной оценке Vr /V0 = -0,2. Вблизи выхлопа ( Z =2,0), увеличиваются радиальные токи, направленные к оси.

В целом, в результате расчетов установлена сильная чувствительность структуры потоков и полей скоростей к значениям параметров крутки завихрителей S1, S2 и к высоте рабочей части аппарата H. В режиме разделения потоков вблизи нижнего завихрителя появляется область с радиальной скоростью Vr > 0 (частичный срез туннельного эффекта на рис. 1, в). При этом прилегающая периферийная область потока имеет значение радиального стока примерно в 2 раза меньшее, чем при K = 0,8. Эти явления способствуют процессу улавливания частиц и, следовательно, повышают эффективность пылеуловителя.

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

В наиболее общем виде модель вихревых течений получена из уравнений Громеки-Ламба и имеет вид 2 1 1 0 + + c2 + c1 + [0 + 1( + 0)]= -c20 - . (10) 2 Здесь: {,,}– безразмерные цилиндрические координаты; – отклонение функции тока от граничного условия на нижнем завихрителе; 0,ci,i – параметры задачи. Аналитическое решение получено в виде ряда (,) = ){Dn 2 exp(- C 2 2)GnF(1- n (4C),2,c 2 Sh(n) Sh(nH )- n=- - 2n AnCh-1(n H 2)Sh[n ( - H ) 2] Sh(n 2)}, (11) где F – функция Куммера, C = i 1, n = n - c2, коэффициенты An, Dn,Gn определяются из граничных условий на торцах аппарата, n – корень уравнения F(1- n (4C),2,c)= 0.

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

Предложены также структурные модели на основе (10), определение параметров в которых выполнено по методу М.А. Гольдштика. Решения получены в виде рядов (11) и описывают только структуры течений в цилиндрических или кольцевых аппаратах. Установлено, что уменьшение кольцевого зазора или введение фильтрации устраняет неоднородности течений.

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

В работе предлагается подход к исследованию закрученных гетерогенных потоков в вихревых устройствах с единых позиций Эйлеро-Лагранжева описания гидромеханики течений. Уравнения движения плотной фазы приведены к безразмерному виду. В качестве масштаба скорости принята безразмерная скорость в плане аппарата V0, масштаба размера - радиус аппарата R, масштаба времени - TA = R V0. Уравнения движения частицы получены в виде:

dWz St = B f1 (Vz -Wz )-WB (12) d W dWr St = B f1 (Vr -Wr )+ St (13) d r dW W Wr St = B f1 (V -W )+ St (14) d r где: W и V - безразмерные скорости частицы и газа; - безразмерное время; St - число Стокса St = TA ; - время релаксации частицы = m d (18 µ); WB - безразмерная скорость витания; B = Cd Cd St - относительный коэффициент сопротивления частицы; Cd St - коэффициент сопротивления частицы, рассчитанный по закону Стокса Cd St = 24 Red ; m,d - плотность и диаметр частиц; Red - число Рейнольдса частицы.

Анализ системы (12-14) выявил основные закономерности движения и область возможного применения критериального подхода. При условиях B =1, WB = const и малых значениях числа Стокса система допускает структурное упрощение. В этом случае решения исходной системы приближенно описываются решением вырожденной системы квазистационарного приближения (правые части уравнений (12-14)). Такое движение частиц будет обладать свойством подобия траекторий и позволяет применить критериальный подход.

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

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

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

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

В процедуре идентификации производится пошаговое построение цепочки моделей Zn, аппроксимирующих C( ) - кривую отклика (здесь n - номер этапа, - нормированное по среднему время пребывания материала в аппарате). Выбор моделей проводился на основе F - критерия и с помощью анализа приближенных функций интенсивности L( ). Функции L( ) рассчитывались путем кубической сплайн-интерполяции точек кривых отклика C( ). По интерполированной кривой C( ) определялась интегральная функция F( ) и, далее, приближенная функция интенсивности L( ) = -d[ln(1- F( ))]/ d.

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

Аналитические решения для первых пяти моделей легко компонуются из известных решений для соответствующих потоков с учетом соотношения между передаточной функцией системы W(s) и передаточными функциями отдельных потоков Wi(s).

Решение для смешанной циркуляционной модели получено в виде ~ k K Ak, Rkj-1 exp(- Tk, j((1 + Rk )/(1 - Bk ))Nk)Ф(Tk, j) C( )= j k , (15) (1 - Bk )Nk j Г(Nk j) k =1 j= (Nk (Nk где Ak, j =(1 + Rk -1) j) Nk Nk j Tk j-1), Tk, j = ( - ( j - 1) Bk / Rk ), ~ K - число контуров в системе, Nk – число секций контура k в ядре потока, k – доля общего потока в контуре, Rk = Qr,k Qk – относительный поток рециркуляции, Bk – объемная доля блока идеального вытеснения, Mk – объемная доля секции в ядре потока,k - число членов ряда обратного преобразования Лапласа, Г - гамма-функция, Ф - функция Хевисайда.

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

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

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

Средний размер частиц d*, попавших в угловую точку цилиндрической границы зоны однородности потоков, характеризует эффективность разделения равную 50 %. Эффективность улавливания, время пребывания и среднее число Рейнольдса различаются для частиц, подаваемых в аппарат по разным входам.

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

Выполнена проверка возможности применения критериального подхода на основе критериев квазистационарного приближения: B 1, St <<1. Оценка критериев проводилась осреднением по траекториям для критического диаметра частиц d*. Показано, что критериальное обобщение наиболее обосновано можно применять в диапазоне характерных размеров частиц 0 d* 30 мкм.

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

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

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

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

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

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

слабого и сильного захвата. Определены верхние границы размеров частиц для зоны сильного захвата, составляющие при V0 = 4 м/с и K=0,39 соответственно 405 и 300 мкм для первого или второго каналов. Установлено, что в режиме разделения потоков, крупные частицы могут иметь низкочастотные колебания.

В стохастической модели использована так же система (12–14), однако r r r скорость газовой фазы моделировалась в виде V = U + ui, где U – осредненная, а ui – пульсационная составляющие скорости газа. Пульсационная составляющая ui рассчитывалась с помощью метода стохастического моделирования:

___ ___ 2 ui = i (ui'), (ui') = Cu [ReT µi2 fu (r)] L2(r). (16) где i – случайные числа нормированного распределения Гаусса, i = (z,r,) – индекс координат, L(r) и f (r) – демпфирующие функции, аналогичные корреляциям для масштаба вихрей в моделях турбулентности, Cu – параметр задачи.

Взаимодействие гетерогенной частицы с турбулентными пульсациями моделировалась на основе схемы сферического вихря А.Д. Госмена и Е. Ионидиса. Принималось, что газовый вихрь, сохраняет пульсационную скорость в течение времени своего существования. Время контакта частицы с вихрем определялось как наименьшее из времени жизни вихря T1 и времени пребывания частицы в вихре T2. Для T1 и T2 предложены уравнения ___ r r ' T1 = L(r) (u1), T2 = L(r) V -W. (17) Введены также пристеночные функции пограничного слоя для частиц, моделирующие эффекты взаимодействия с шероховатой стенкой и вязким подслоем. Взаимодействие со стенкой, помимо удара частиц, моделировалось в виде численного аналога односторонней пелены с радиусом Rpl = (R - k ), pl где – шероховатость стенки.

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

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

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

На рисунках 3, б и 3, в показаны траектории частиц диаметром 2 мкм, введенных по первому каналу, но попавших в различные структуры потоков.

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

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

а) б) в) Рис.3 Стохастические траектории движения частиц в аппарате в плоскости. H =2,15; =0,6; =5,2.

= S1 Sp а) К=0,62, 2 канал, dp=6 мкм; б) К=0,62, 1 канал, dp=2 мкм; в) К=0,39, 1 канал, dp=2 мкм.

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

Для оценки влияния структур потоков в ядре течения на фракционную эффективность аппарата предложено формировать и исследовать функции фракционной эффективности осаждения (ФЭО) на стенку аппарата. В ряде случаев эффекты ЛПС будут незначительно влиять на фракционную эффективность улавливания (ФЭУ), и тогда, функции ФЭО и ФЭУ совпадают.

Таким образом в основу расчетов ФЭО положена турбулентная модель несущей фазы (1-9) и стохастическая модель движения частиц (16-17). Анализ движения групп частиц осуществлялся по осреднению траекторий пробных частиц, представляющих группу. Все исходы движения кроме уноса в стохастической модели рассматривались как осаждение на стенку. Для одного вида пробной частицы осреднение проводилось по 700 траекториям. На рис. 4, а и 4, б показаны результаты теоретических расчетов функций фракционного осаждения ФЭО и экспериментальные кривые ФЭУ для идентичных условий. Как видно из рисунка 4, а при цилиндрической структуре (рис. 1, б) второй канал данного аппарата работает более эффективно, чем первый, что подтверждают и эксперименты. В целом, если иметь в виду, что в данном случае расчеты прогнозирующие, то можно считать качественное и количественное совпадение с экспериментом достаточно хорошим.

а) б) Рис.4 Результаты расчетов ФЭО в различных структурах потоков.

а) К=0,62 : 1 канал: – расчет, – эксперимент; 2 канал: – расчет, – эксперимент;, – соответственно расчетные ФЭУ по базовой 2P модели турбулентности [*] и экспериментальные данные[**]. б) К=0,39 : обозначения – см. а).

Таким образом, при цилиндрической структуре потоков методика позволяет рассчитать основную характеристику центробежного разделения – фракционную эффективность улавливания. На рис. 4, а приведены также расчеты близкого по характеристикам циклона, выполненные с помощью базовой CFD – модели [*] и соответствующие экспериментальные данные [**]. Как видно из рис. 4, а эта модель недостаточно хорошо описывает дисперсию ФЭУ.

На рис. 4, б даны теоретические функции ФЭО и экспериментальные фракционные характеристики для структуры разделения потоков. В этом случае характеристики каналов принципиально меняются: расчеты и эксперименты показывают, что при этой структуре более эффективен первый канал. Однако совпадение с экспериментом дает лишь функция ФЭО второго канала. ФЭО первого канала в данном случае прогнозирует очень высокие эффективности осаждения, в том числе для тонкодисперсных частиц в диапазоне 1–5 мкм.

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

Отличие расчета от эксперимента вызвано, в основном, не учитываемым эффектом скопления и взаимодействия частиц в ЛПС в районе точки отрыва [*] –Griffiths W.D., Boysan F. //J. Aerosol Science. 1996. V. 27. №2. P. 281– 304.

[**] – Kim C.H., Lee J. W. //J. Aerosol Science. 2001. V. 32. №2. P. 251– 269.

потока 2 (рис. 3, в).

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

В качестве обобщенных характеристик структуры потоков предложены параметры K*, H, K. Основной обобщенной характеристикой вихревого апp парата ВЗП предложено соотношение расходов потоков, K* при котором происходит слияние точек присоединения и отрыва первого потока от стенки (совпадение точек 1 и 2, обозначенных на рисунке 1, в). Величина H обозначает p высоту зоны разделения, то есть расстояние между точками 1 и 2. Через K обозначено такое соотношение потоков, что при K > K доля K второго потока всасывается восходящий поток в малой окрестности поверхности z = 0. В рабочем диапазоне соотношения расходов {Kmin, Kmax} выделены три структуры течения: разделения потоков Kmin < K < K*; структура равномерного перехода потоков K* < K < K (цилиндрическая); структура последовательного прохождения потоков (режим частичного байпаса) – K < K < Kmax.

Обобщенную характеристику, связывающую конструкцию аппарата и структуру потоков, предложено искать в виде зависимости K* = f (S1,S2, H ).

Определение этой характеристики, а также высоты зоны разделения проводилось по результатам экспериментов на 15 модификациях конструкции с различными значениями высоты аппарата и параметрами круток завихрителей. В результате обобщения опытных данных для K и H получены уравнения p 0,33 0, K = 0,77 S1 H 1+ 0,132 S2 + 0,77 S1 H, (18) H = 4,72 (K - K)0,67, (19) P Формулы (18) и (19) справедливы в диапазоне аргументов: S1 = 0,6 - 3,0 ;

S2 = 3,0 -12,0; H = 2 - 4. Среднеквадратичная погрешность аппроксимации для уравнения (18) составляет 2,7%, а для уравнения (19) – 6%.

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

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

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

Выполнено исследование структур потоков для следующих аппаратов:

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

а) б) Рис.5 Кризисные явления в вихревых аппаратах.

а) распад вихря в аппарате ВЗП: расчетные линии тока нанесены на экспериментальную визуализацию структуры вращающихся частиц; б) Кризис структуры потоков в циклоне: – экспериментальная L( )– функция интенсивности в точках замеров и между ними; – двухпоточная секционная модель; – ячеечная модель;

– относительная ошибка L экспериментальной функции интенсивности.

Расчет структуры двухфазного течения в циклоне показал, что для всех режимов и конструкций (кроме длинного аппарата) адекватной является ячеечная модель с числом ячеек N =35–50. Для длинного аппарата ( H =12 ) адекватV0* на двухпоточная модель, возникающая при критической скорости.

На рисунке 5, б представлены экспериментальные и расчетные кривые L – функций, характеризующие структуру потоков в длинном циклоне при H = 1.5 м, V0 =1,28 м/с. Как следует из этих результатов, в данном случае происходит распад потока на две ветви, причем вторая ветвь имеет большее время пребывания и дисперсию. Такими характеристиками должны обладать частицы, вынесенные из пограничного слоя, но затем вернувшиеся в основной поток и V0* уловленные. Получено уравнение для критической скорости циклона V0* =1,6 103 µ (Hц)-0,3 D0,2 S {2g (S -1)}, (20) p где , – плотность материала и газа; µ – вязкость газа; Hц = hц / D – p относительная высота цилиндрической части аппарата; S – параметр крутки.

Установлено, что в комбинированном аппарате также может происходить разделение течения на две ветви с разными характеристиками, однако в данном случае причина в высоте переливного порога. При увеличении высоты порога до величины h = 0,075 исчезает разделение потоков и максимально возрастает среднее время пребывания.

Показано, что в дисковой вихревой сушилке, структуры потоков зависит от расхода газа через камеру Q, причем вид модели изменяется при некотором критическом значении Q*. При Q = Q Q* 1 структура близка к идеальному смешению. В диапазоне 1 Q 1.5 происходит интенсивная перестройка структуры потоков, число ячеек увеличивается до N = 3.2. Далее, при расходах 1.5 Q 2 число ячеек N 4 и в системе возникают и увеличиваются с ростом Q циркуляционные колебания.

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

Представлены расчетные и экспериментальные данные, характеризующие циркуляционную структуру потоков в аппарате фонтанирующего слоя с дополнительными вводами теплоносителя. Установлено, что все стадии циркуляционных процессов хорошо описывает двухконтурная циркуляционная модель. При табличном значении критерия FT = 4,7 эта модель имеет F = 3,89 и является адекватной. Соотношение времени циркуляций в контурах составляет 2 1 = 2.17, что говорит о значительном снижении неравномерности в пристенной зоне при применении дополнительных вводов с закруткой потоков.

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

Экспериментальные исследования проводились на 20 модификациях пылеуловителей ВЗП и ВП; комбинированной системе пылеочистки; модификациях многофункционального аппарата для сушки тонкодисперсных материалов. Для определения характеристик пылеуловителей использовались молотые кварцевые пыли, а также специальные и промышленные пыли. Дисперсность пылей классифицировалась по удельной поверхности, составлявшей для пыли P4 – 10500 см2/г, и, соответственно для пылей P - P3 – 6500 см2/г – 3500 см2/г.

Для подготовки и анализа проб использовался лабораторный комплекс приборов “Fritsch”, фотоседиментограф “Analysette”, лазерный измеритель грансостава “Malvern”. Лабораторные аппараты для исследования технологических характеристик имели широкий диапазон варьирования параметров и были согласованы с конструкцией аэродинамического стенда (глава 3).

На рис. 6, а даны результаты экспериментального исследования общей эффективности улавливания при различных структурах потоков в лабораторном аппарате с характеристиками: D = 0,107 м; S1 = 0,6 ; S2 = 5,2 ; K* = 0,45.

а) б) Рис.6 Выходные характеристики вихревых аппаратов при различных структурах потоков, каналах подачи материала и конструкции.

а) Эффективность улавливания и в аппаратах ВЗП: ВЗП-107: 1 канал: – пыль P1;

– P2; – P4; 2 канал: – пыль P1; – P2; – P3; – расчет эффективности для пыли P2, канал 1.

Дисперсии ФЭУ 1 и 2 каналов: – 1; – 2; – дисперсия ФЭУ канала №при S1 =1,8. б) Расчетные и экспериментальные значения потерь давления в аппарате ВЗП107 при изменении крутки периферийного завихрителя.

Как видно из представленных данных существует диапазон соотношения расходов K1 K K2, (включающий значение K*), в котором общая эффективность улавливания первого канала 01 превосходит эффективность второго канала. Существование такого диапазона экспериментально подтверждает эффект туннельного турбофореза. Действительно, при значениях K2 K уменьшается зона распада вихря (рис. 1, б), ослабляется движущая сила туннельного турбофореза для первого потока, но при этом увеличивается движущая сила центробежного разделения для обоих каналов за счет большей крутки второго потока. В итоге в этой области эффективность второго канала становится выше, чем эффективность первого.

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

Установлено, что фракционная эффективность улавливания каналов описывается интегральной функцией F логарифмически нормального распределения. Фракционная степень очистки в аппарате в целом описывается уравнением (d) =100 (1- K) Flg(d / d1) lg(1) + K Flg(d / d2) lg(2), (21) где dk, k, k=1,2 – параметры функций фракционной эффективности улавливания по каналам аппарата.

Исследование улавливания в конструктивной модификации аппарата с увеличенным параметром крутки нижнего ввода S1 =1,8 дает увеличение эффективности при K = 0,6 и уменьшение при K = 0,4 по сравнению с данными на рис. 6, а. В этой конструкции параметр K* = 0,55, поэтому при K = 0,4 аппарат уже работает в режиме вихревого слоя, что вызывает уменьшение 0. Вместе с тем, повышение S1, как видно из рис. 6, а улучшает остроту разделения.

Таким образом, для разработки эффективных режимов разделения и для многофункциональных процессов надо согласованно варьировать предложенные обобщенные параметры K*,S1,S2,H, K.

p Исследованы зависимости потерь давления от различных параметров. Установлено, что они также зависят от вида структуры потоков и степени закрутки, то есть от обобщенных переменных K*,S1,S2,H. Значение потерь давления p от K имеет экстремальный характер, однако точка минимума энергозатрат в общем случае не совпадает с кризисным параметром K*.

Разработана модель потерь давления в вихревом аппарате, согласованная с положениями вариационного принципа (глава 2). Получено, что потери энергии пропорциональны действующему в данном сечении Z параметру крутки SZ, и обратно пропорциональны характерному радиусу R*t, связанному с радиусом выхлопной трубы. Для расчета коэффициентов сопротивления каналов получены уравнения о 2 Тk = R*tФТk R[(1- hk ФТk R*t)-1]+ fk0(ФТk ), k=1,2 (22) ФТk = Фвх (1+ hк Фвх / R*t), k = 1,2 (23) (к) (к) Фвх = 3S1R*t (2R1), Фвх = S2R*tR2 (R2 - RT)Rвх, (24) (1) (2) где S1 и S2 – конструктивные параметры крутки завихрителей, fk0 – функция, корректируемая по экспериментальным данным, RT – радиус выхлопной трубы, – коэффициент трения, h1 = H, h2 = H + hT. Как видно из рисунка 6, б наблюдается хорошее согласование расчетов по модели с экспериментом.

Для характеристики эффективности разделения предложен критерий F =100 P 0, где P –потери давления.

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

Как установлено в результате расчетов, режим распада вихря дает возможность захвата и удержания как частиц, размером от 1 до 20 мкм, так и частиц в диапазоне 20–500 мкм. Как показали эксперименты в этом режиме возрастает среднее время пребывания. Неоднородная структура обуславливает также интенсификацию процесса за счет нестационарного взаимодействия фаз вдоль траектории движения материала. При этом относительное число Рейнольдса, Red, увеличивается, особенно в зонах циркуляции частиц.

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

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

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

Описание процесса в зонах основано на уравнениях (12-14) и (21). Основные расчетные уравнения получены в виде С1 µ(1- k )Q2(2 fS -1)ln(LK ) d1 =, (25) R0M H (Q14 fS -2 -1) С2R02 µQ2 R1+ d2 = f2 , (26) 2 c M H2 4R02(R2 - R02) ' 2R0 Q11/(1-kK ) QQ1 = 1+, Lk =, f2 = 0,75 - 0,5(R02 / R)2 - 0,25(R02 / R)4, (27) ' 2 Q1 R1 + R' ' где Q1 = Q1 + k Q2, Q2 = (1- k )Q2 ; величины 0 и с зависят от интегральных параметров крутки по уравнению (23), C1, C2 – параметры задачи.

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

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

а) б) Рис.7 Вихревой аппарат типа ВП.

а) Принципиальная схема пылеуловителя типа ВП: 1 – дополнительный завихритель 2-го канала; 2 – вставка пережима; 3 – завихритель 1-го канала; б) Фракционные характеристики аппарата:, – параметры ФЭУ приведенные к характеристикам d=типового пылеуловителя ВЗП-Т:, – соответственно и первого каd=нала;, – и второго канала; K –соотношение расходов потоков.

d=Комплексные исследования, выполненные на ряде модификаций аппарата ВП диаметром D = 0,2 по стандартной методике испытаний пылеуловителей, позволили определить рациональные конструктивные параметры и исследовать характеристики полученной конструкции. Характеристики аппарата ВП позволяют, при одинаковых параметрах, снизить выбросы в атмосферу для типичных промышленных пылей в 1,4 – 1.8 раза. Аппарат имеет в 2 раза больший диапазон регулирования соотношения расходов, K = (0,35 - 0,65) чем типовая конструкция ВЗП-Т и, в этом плане занимает промежуточное положение между пылеуловителями ВПУ и ВЗП. При соотношении K = 0,5 обеспечивается минимум энергозатрат. Высота аппарата и металлоемкость примерно в 1,25 раза меньше чем у ВЗП-Т. Аппарат устойчив к колебаниям как общей нагрузки, так и соотношения расходов газа.

По результатам исследования разработаны промышленные типоразмеры аппаратов ВП диаметром от 200 до 1200 мм и установки ВП-700х2 и ВП-800х3, обеспечивающие очистку газов в диапазоне производительности от 500 до 22000 м3/час. Разработаны и внедрены аппараты ВП-300 в производстве трехсернистой сурьмы, пылеуловители ВП-700, установка ВП-700х2, комбинированная установка ВП-800х3 в литейном цехе машиностроительного предприятия.

В ряде производств внедрены пылеуловители типа ВЗП. В производстве красок и жировых смесей внедрена установка ВЗП-600х3, обеспечивающая очистку 22000 м3/час газа. В производстве синтетических моющих средств внедрены три пылеочистные установки, включающие семь аппаратов ВЗП-700. В производстве сложных минеральных внедрены пылеуловители ВЗП-2000, обеспечивающие очистку газов от пыли при расходе воздуха 70000-1200м3/час. Полученные в работе рациональные значения парметров реализованы в пылеуловителях ВЗП-2000 и ВЗП-2000У. Результаты работы были использованы при внедрении пылеуловителей с разделением потоков ИВРП и ВЗП-ЦРП 1200/700, при проектировании микрогранулятора ВЗП-300-МГ и сушилки СВЗП-600. Даны рекомендации по рациональным значениям режимно – конструктивных параметров пылеулавливающих и тепломассообменных аппаратов ВП и ВЗП.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ 1. Разработаны научные основы прикладной гидродинамики вихревых аппаратов с закрученными гетерогенными потоками, включающие теорию турбулентных вихревых течений (физические и математические модели), ЭйлероЛагранжеву детерминированную и стохастическую модели движения дисперсной фазы. Разработанные методики, основаны на учете аппаратных особенностей при применении фундаментальных методов для инженерных расчетов.

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

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

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

4. Предложены обобщенные параметры, связывающие неоднородную гидродинамику и конструкцию вихревых аппаратов: модельные параметры крутки Фм, критическое соотношение расходов потоков K*, высота зоны разделения H, критический диаметр частиц d* и другие. В результате обобщения опытp ных данных получены уравнения для параметров K и H.

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

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

7. Предложена методика расчета функций фракционной эффективности осаждения – ФЭО, характеризующих влияние структур потоков на эффективность сепарации частиц; разработаны эффективные методы определения ФЭУ.

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

9. Разработан ряд экспериментальных стендов, на которых выполнены испытания 20 модификаций вихревых аппаратов в широком диапазоне параметров:

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

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

11. Разработан и исследован новый вид вихревых пылеуловителей типа ВП, с рабочим соотношение входных расходов 1:1, широким диапазоном по соотношению расходов, малочувствительных к колебаниям нагрузок; типоразмеры аппаратов ВП обеспечивают очистку газов в диапазоне от 500 до 220м3 /час. Разработаны рекомендации по конструктивным характеристикам и эффективным режимным параметрам для пылеулавливающих и многофункциональных аппаратов ВЗП, комбинированных аппаратов.

12. По результатам работы внедрено в различных производствах более 15 вихревых аппаратов.

Условные обозначения, нерасшифрованные в тексте:

DT, ( RT ) – диаметр (радиус) выхлопной трубы; V – скорость газа по отношению к плановой скорости V0.

Основные научные результаты диссертации опубликованы в 1печатных работах, в том числе:

1. Белоусов А.С., Сажин Б.С., Отрубянников Е.В. Структура потоков в аппаратах со взвешенным слоем // Химическая технология. 2008. т. 9. №7. С. 332– 336.

2. Белоусов А.С., Сажин Б.С. Гидродинамические неоднородности вихревого слоя газовзвеси // Изв. вузов. Технология текстильной промышленности.

2007. № 6С. С. 140–144.

3. Белоусов А.С., Сажин Б.С. Закономерности структур течений в аппаратах для обработки волокнообразующих полимеров при активных гидродинамических режимах // Химические волокна. 2007. № 6. С. 40–43.

4. Сажин Б.С., Кочетов Л.М., Белоусов А.С. Удерживающая способность и структура потоков в вихревых аппаратах // Теоретические основы хим. технологии. 2008. т. 42. №2. С. 125–135.

5. Белоусов А.С., Сажин Б.С. Циркуляционные течения в организованном взвешенном слое //Изв. вузов. Технология текстильной промышленности. 2007.

№ 3С. С. 115–118.

6. Белоусов А.С., Сажин Б.С., Лопаков А.В. Кризисные явления в винтовых потоках газовзвесей //Изв. вузов. Технология текстильной промышленности.

2007. № 2. С. 110–114.

7. Белоусов А.С., Сажин Б.С., Сажин В.Б. и др. Распад потока газовзвеси в циклоне // Успехи в химии и химической технологии. Т. XXI. М.: РХТУ им.

Д.И. Менделеева. 2007. №9 (77). С. 111–113.

8. Сажин Б.С., Сажина М.Б., Белоусов А.С. и др. Основные характеристики и классификация дисперсных волокнообразующих полимеров как объектов сушки // Изв. вузов.Технология текстильной промышленности. 2007. № 3. С.

111–113.

9. Белоусов А.С., Сажин Б.С., Лопаков А.В. Обработка волокнообразующих полимеров в винтовом потоке // Химические волокна. 2007. № 6. С. 40–43.

10. Белоусов А.С., Сажин Б.С., Сажин В.Б. и др. Исследование структуры закрученных потоков с дисперсной фазой // Сборник тез. докл. междун. конф.

по химич. технологии ХТ’07. Т.2/ Под ред. А.А.Вошкина, Н.Н.Кулова, А.И.Холькина, Е.В.Юртова. М.: ЛЕНАНД. 2007. С. 134-136.

11. Сажин Б.С., Белоусов А.С. Гидродинамические характеристики двухфазных течений в вихревых устройствах.//Вестник МГТУ. М.: МГТУ им.

А.Н.Косыгина. 2007. С. 93-97.

12. Белоусов А.С., Сажин Б.С., Лопаков А.В. и др. Влияние параметров процесса на время пребывания дисперсной фазы в вихревом аппарате // Успехи в химии и химической технологии. Т. XX. М.: РХТУ им. Д.И. Менделеева. 2006.

№10 (68). С. 92–94.

13. Сажин Б.С., Сажина М.Б., Белоусов А.С. и др. Комплексный анализ дисперсных волокнообразующих полимеров как объектов сушки во взвешенном слое // Химические волокна. 2007. № 3. С. 24–26.

14. Белоусов А.С., Сажин Б.С., Лопаков А.В. и др. Гидродинамическое перемешивание дисперсной фазы в вихревом потоке // Изв. вузов. Технология текстильной промышленности. 2006. № 6. С. 104–109.

15. Белоусов А.С., Сажин Б.С., Отрубянников Е.В. Турбулентное смешение потоков в вихревом аппарате с двумя входами // Успехи в химии и химической технологии. Т. XX. М.: РХТУ им. Д.И. Менделеева. 2006. №10 (68). С. 94–96.

16. Сажин Б.С., Сажина М.Б., Белоусов А.С. и др. Классификация волокнообразующих полимеров как объектов сушки // Химические волокна. 2007. № 3.

С. 21–23.

17. Белоусов А.С., Сажин Б.С. Структура потоков в вихревых устройствах // Изв. вузов. Технология текстильной промышленности. 2006. № 5. С. 98–103.

18. Белоусов А.С., Сажин Б.С., Лопаков А.В. Дисперсия фракционной эффективности разделения в центробежных пылеуловителях // Успехи в химии и химической технологии. Т. XX. М.: РХТУ им. Д.И. Менделеева. 2006. №10 (68).

С. 96–98.

19. Белоусов А.С., Сажин Б.С., Лопаков А.В. и др. Структура потоков в циклоне // Успехи в химии и химической технологии. Т. XX. М.: РХТУ им. Д.И.

Менделеева. 2006. №9 (67). С. 121–123.

20. Белоусов А.С., Сажин Б.С.. Радиальный сток в центробежных пылеуловителях // Изв. вузов. Технология текстильной промышленности. 2006. № 4. С.

96–100.

21. Белоусов А.С., Сажин Б.С., Лопаков А.В. Влияние режима работы и конструкции на характеристики вихревого аппарата //Успехи в химии и химической технологии. Т. XX. М.: РХТУ им. Д.И. Менделеева. 2006. №2 (60). С. 94– 97.

22. Белоусов А.С., Сажин Б.С. Поля скоростей в вихревых аппаратах // Изв.

вузов. Технология текстильной промышленности. 2006. № 2. С. 100–105.

23. Белоусов А.С., Сажин Б.С., Лопаков А.В. Гидродинамическая модель пылеуловителя с разделяющимися потоками //Успехи в химии и химической технологии. Т. XIX. М.: РХТУ им. Д.И. Менделеева.2005. №10 (58). С. 117–119.

24. Белоусов А.С., Сажин Б.С. Закрутка потоков в вихревых аппаратах // Изв. вузов. Технология текстильной промышленности. 2005. № 3. С. 96–100.

25. Белоусов А.С., Сажин Б.С., Лопаков А.В. К расчету потерь давления в вихревом аппарате // Успехи в химии и химической технологии. Т. XIX. М.:

РХТУ им. Д.И. Менделеева. 2005. №10 (58). С. 62–64.

26. Белоусов А.С., Сажин Б.С. Диффузионная модель перемешивания в технологических аппаратах при малых числах Пекле // Изв. вузов. Технология текстильной промышленности. 2005. № 2. С. 96–100.

27. Белоусов А.С., Сажин Б.С., Кочетов Л.М. и др. Структура потоков двухфазной полидисперсной смеси в вихревой сушилке //Успехи в химии и химической технологии. Т. XIX. М.: РХТУ им. Д.И. Менделеева. 2005. №10 (58). С.

77–79.

28. Белоусов А.С. Адекватность моделей перемешивания структуре потоков в химико-технологических аппаратах // Современные наукоемкие технологии и перспективные материалы текстильной и легкой промышленности: сб. матер.

междунар. конф. Ч. 1. Иваново: ИГТА. 2004. С. 103–104.

29. Белоусов А.С., Сажин Б.С., Кочетов Л.М. и др. Аэродинамика турбулентного потока в вихревой сушилке // Успехи в химии и химической технологии.

Т. XVIII. М.: РХТУ им. Д.И. Менделеева. 2004. №6 (46). С. 126–128.

30. Белоусов А.С., Сажин Б.С., Миронов А.В. Интегральные параметры крутки потоков в вихревых пылеуловителях // Успехи в химии и химической технологии. Т. XVIII. М.: РХТУ им. Д.И. Менделеева. 2004. №7 (47). С. 86–89.

31. Belousov A., Sazin B. Application of Guided Vortex Breakdown for Drying and Separation of the Powder in Vortex Cyclone // Proceeding of The Second Nordic Drying Conference (NDC–03). CD. Copenhagen. Denmark. 2003. P. 1-5.

32. Белоусов А.С., Кочетов Л.М., Сажин Б.С. Гидродинамическая структура потоков в вихревой сушильной камере // Успехи в химии и химической технологии. Т. XVII. М.: РХТУ им. Д.И. Менделеева. 2003. №13 (38). С. 94–97.

33. Белоусов А.С., Яковлев А.Ю. Экспериментальное исследование процесса разделения газовзвесей в вихревом аппарате // Разработка и внедрение аппаратов с активными гидродинамическими режимами в текстильной промышленности и производстве химических волокон: межвуз. сб. науч. тр. М.: МТИ им. А.

Н. Косыгина. 1991. С. 43–47.

34. Сажин Б.С., Белоусов А.С., Коротченко С.И. и др. Использование текстильных фильтровальных материалов для предварительной концентрации очищаемых запыленных газов.// Актуальные проблемы техники и технологии в текстильной и легкой промышленности: межвуз. сб. науч. тр. М.: МТИ им. А.

Н. Косыгина. 1991. С. 99–102.

35. Сажин Б.С., Белоусов А.С., Лапшин А.Б. Математическая модель и метод расчета течения газа через фильтровальную ткань // Изв. вузов. Технология текстильной промышленности. 1990. № 2. С. 83–87.

36. Белоусов А.С. Расчет и моделирование вихревых аппаратов и установок для пылеочистки // Очистка газовых выбросов промышленных предприятий:

сб. тр. всесоюзн. конф. Тольятти: 1990. С. 116-117.

37. Белоусов А.С., Яковлев А.Ю. Исследование и разработка методов расчета характеристик классификации и разделения газовзвесей в вихревых аппаратах.// Разработка и совершенствование технологических процессов и худож.

оформления изделий в текстильной и легкой промышленности: межвуз. сб. науч. тр. М.: МТИ им. А. Н. Косыгина. 1988. С. 47–50.

38. Сажин Б.С., Белоусов А.С., Коротченко С.И. Моделирование вихревых течений газа в технологических аппаратах и устройствах.// Изв. вузов. Технология текстильной промышленности. 1990. № 1. С. 81–85.

39. Белоусов А.С. Метод расчета процесса разделения газовзвесей в вихревых аппаратах с переменной структурой потоков // Повышение эффективности технологических процессов и оборудования в текстильной промышленности и производстве химич. волокон и разработка систем пылеулавливания: межвуз.

сб. науч. тр. М.: МТИ им. А. Н. Косыгина. 1988. С. 35–37.

40. Сажин Б.С., Белоусов А.С. Степень разделения газовзвесей в вихревых аппаратах // Современные машины и аппараты химических производств: тез.

докл. всесоюзн. конф. Чимкент: 1988. часть 3. С. 115-116.

41. Белоусов А.С., Яковлев А.Ю. Определение фракционных характеристик классификаторов и пылеуловителей со встречными закрученными потоками // Повышение эффективности технологич. процессов и оборудования в текстильной пром-ти и производстве химич. волокон и разработка систем пылеулавливания: межвуз. сб. науч. тр. М.: МТИ им. А. Н. Косыгина. 1988. С. 3–5.

42. Белоусов А.С. Прогнозирование структуры течений в аппаратах со встречными закрученными потоками.// Методы кибернетики в химии и химич.

технологии: тез. докл.2-го всесоюзн. совещания. Грозный: 1984. С. 48.

43. Белоусов А.С., Фокин И.Ф. Режимная оптимизация промышленного аппарата со встречными закрученными потоками.// Аппараты с активными гидро динамич. режимами для текстильной промышленности и производства химич.

волокон: межвуз. сб. науч. тр. М.: МТИ им. А. Н. Косыгина. 1983. С. 33–35.

44. Сажин Б.С., Белоусов А.С. Турбулентные и вихревые течения в аппаратах со встречными закрученными потоками // Аппараты с активными гидродинамич. режимами для текстильной промышленности и производства химич. волокон: межвуз. сб. науч. тр. М.: МТИ им. А. Н. Косыгина. 1983. С. 17–21.

45. Сажин Б.С., Лукачевский Б.П., Белоусов А.С. Математическое моделирование аппаратов со встречными закрученными потоками // Всесоюзн. совещ.

по математич. моделированию и управлению высокотемпературными процессами в циклонных и вихревых аппаратах: тез. докл. Одесса. 1980. С. 13–14.

46. Белоусов А.С. Расчет коэффициентов переноса многокомпонентных газовых смесей // Гидродинамика и тепломассоперенос в технологических процессах: межвуз. сб. науч. тр. М.: МТИ им. А. Н. Косыгина. 1981. С. 95–98.

47. Патент РФ на изобретение № 2312286. Многосекционная вибрационная сушилка для дисперсных и адгезионных материалов / Б.С. Сажин, О.С. Кочетов, А.С. Белоусов и др.; опуб. Б.И. № 34 от 10.12.07.

48. А.с. 1585018 СССР. Сепаратор / Б.С. Сажин, А.С. Белоусов, А.Ю Яковлев и др.; опубл. Б.И. № 30 от 15.08.90.

49. А.с. 1595570 СССР. Вихревой пылеуловитель / Б.С. Сажин, А.С. Белоусов, С.И. Коротченко и др.; опубл. Б.И. № 36 от 30.09.90.

50. А.с. 1623731 СССР. Способ очистки газа от пыли / Б.С. Сажин, А.С. Белоусов, С.И. Коротченко; опубл. Б.И. № 4 от 30.01.91.






© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.