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

Петровичева, Елена Александровна

Турбулентное течение смешивающихся жидкостей в малогабаритных трубчатых аппаратах химических производств тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

Петровичева, Елена Александровна АВТОР

кандидата технических наук УЧЕНАЯ СТЕПЕНЬ

Казань МЕСТО ЗАЩИТЫ

2006 ГОД ЗАЩИТЫ

01.02.05 КОД ВАК РФ

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

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

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

ПЕТРОВИЧЕВА ЕЛЕНА АЛЕКСАНДРОВНА

ТУРБУЛЕНТНОЕ ТЕЧЕНИЕ СМЕШИВАЮЩИХСЯ ЖИДКОСТЕЙ В МАЛОГАБАРИТНЫХ ТРУБЧАТЫХ АППАРАТАХ ХИМИЧЕСКИХ ПРОИЗВОДСТВ (ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ)

01.02.05 - Механика жидкости,-газа и плазмы

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

Казань - 2006

Работа выполнена на кафедре « Прикладная химия » в Альметьевском государственном нефтяном институте

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

Данилов Юрий Михайлович

Официальные оппоненты - доктор технических наук, профессор

Фафурин Андрей Викторович - доктор физико-математических наук, профессор Зарипов Шамиль Хузеевич

Ведущая организация - Исследовательский центр проблем энергетики

Казанского научного центра РАН

Защита состоится часов на заседании

звета Д 212.1

диссертационного совета Д 212.079.02 в Казанском государственном техническом университете им. А.Н. Туполева по адресу: 420111, г. Казань, ул. К.Маркса 10.

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

Автореферат разослан /£С^/Ти^^Ь1? 2006г.

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

к.т.н., доцент .А.Г.Каримова

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

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

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

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

1. Грант Президента РФ № 96-15-97179 по теме «Моделирование процессов полимеризации при производстве синтетических каучуков».

2. Программа 05 ГКНТ 12 «Разработка методов моделирования и расчета принципиально новых малогабаритных реакторов для осуществления быстрых химических реакций, эффективной теплопередачи и массообмена в турбулентных потоках с проведением опытных • и промышленных испытаний».

3. Программа Республики Татарстан по развитию приоритетных • направлений науки по теме № 19-12/99 (Ф) «Научные основы технологических процессов производства синтетических каучуков. на предприятиях нефтехимического комплекса Республики Татарстан». ■ ' .

4. Программа Республики Татарстан по развитию приоритетных. направлений науки по теме №. 07-7.5 - 27/2001 (Ф) «Ресурсосберегающие и экологически безопасные трубчатые аппараты и технологические процессы для нефтехимической промышленности».

Целями и задачами работы являются:

1) численное моделирование турбулентного течения смешивающихся жидкостей или газов;

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

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

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

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

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

3. Установлено, что за счет изменения формы турбулизатора в трубчатом аппарате можно значительно улучшить качество смешения.

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

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

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

• выбрать наиболее эффективную форму образующей в области горловин секций промышленных аппаратов диффузор — конфузорного типа (ТРДКТ).

• выбрать наилучшую схему подвода компонентов в проточную часть аппарата.

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

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

Апробация работы. Материалы исследования и его результаты неоднократно обсуждались на заседаниях кафедры прикладной химии Альметьевского государственного нефтяного института и докладывались на научных конференциях: на научно-технической конференции « АлНИ-2000», г. Альметьевск, 2001 г; на научной сессии КГТУ им. С.М. Кирова, г.Казань, 2001г.; на III Международной научно-практической конференции ПГУ, г. Пенза, 2001г.; на научно-технической конференции « АлНИ- 2002», г. Альметьевск, 2003г.; на Всероссийской научно-технической конференции «Большая нефть; реалии, проблемы, перспективы», г. Альметьевск, 2001г.; в Альметьевском государственном нефтяном институте, 2003,2004,2005 гг.; в Казанском государственном техническом университете им. Туполева, 2003г.; на

4-й Школе-семинаре молодых ученых и специалистов под руководством акад. В.Е. Алемасова, г.Казань, 2004 г.; на научно-технических семинарах в Казанском государственном технологическом университете, 2003, 2004, 2005 гг.; на Международной научной конференции « Математические методы в технике и технологиях — ММТТ-18», г. Казань, 2005 г.

Публикации. По материалам диссертационной работы опубликовано 10 работ.

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

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

2. Результаты математического моделирования турбулентных течений смесей несжимаемых жидкостей в малогабаритных трубчатых аппаратах диффузор - конфузорного типа (ТРДКТ).

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

Краткое содержание работы

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

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

а б

Г«сЛ.Труб-т*тпцти _ Г) ТТ~TT" 1 TT" Y/"4»

цилиндрического (а) -- Ш Ш tXTlJlIlli

н диффузо)>-«о*»ф>тор1юго (5) типа л ^

К настоящему времени такого типа реактор получил свою конфигурацию благодаря работам академика A.A. Берлина, а также кафедры процессов и аппаратов химической технологии (ПАХТ) Казанского государственного технологического университета. Наиболее детально гидродинамические исследования проведены в работах Р.Г. Тахавутдинова и А.Г. Мухаметзяновой. При исследовании процессов авторы использовали коммерческий пакет PHOENICS.

Предпочтительны следующие параметры трубчатого турбулентного реактора для проведения реакции полимеризации в условиях отсутствия радиального температурного и концентрационного градиентов: скорость течения основного потока У=0,3 +• 1 м/с, 11=0,025 + 0,06 м; Ь=1+1,5 м. При этом лучшие результаты обеспечивает подача потока катализатора перпендикулярно к потоку мономера, а также использование второго поколения трубчатых турбулентных реакторов конфузор - диффузорного типа (рис.1.а). Разработка и внедрение в промышленное производство последних явились базой для создания энерго- и ресурсосберегающих технологий нового поколения реакторов повышенной экологической чистоты не только для весьма быстрых химических, но и ряда массообменных физических процессов. На этой основе разработан принципиально новый способ получения полимеров бутиленов с молекулярной массой от 200 до 120000, позволивший резко упростить технологию и аппаратурное оформление производственного процесса получения олиго- и полиизобутиленов.

На рис. 2 приведены возможные схемы подачи компонентов в рабочую зону МТГА, которые и рассмотрены далее при выполнении данной диссертационной работы.

............

Г1 РЗ — Ч) .»Л 1 1

-----------

Рис. 2

В работах Ю.М.Данилова и И.М.Ильиной диапазон исследованных форм образующей в окрестности горловины был расширен. Было установлено что наибольшей эффективностью обладают МТТА с турбулизаторами в виде плоских круглых диафрагм.

Подводя итог обзору информации о состоянии исследований турбулентных течений в малогабаритных трубчатых аппаратах (МТТА), можно сделать следующие выводы:

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

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

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

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

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

Приведенные выводы позволили сформулировать тему диссертации, ее научную и практическую направленность.

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

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

Наиболее удобное выражение для коэффициента вязкости газовой смеси было дано Уилки

Мт = 1>

1+—У с. ,..

л ¿-I кр ¡к

Здесь суммирование проводится по всем компонентам смеси и ^ коэффициент вязкости ьго компонента. С точностью до величин первого порядка

[1 + Ц/т.Г

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

Пусть в объеме V содержится (ш1+гп2) кг смеси компонентов 1 и 2, имеющих плотности р! и рг соответственно. Тогда У=т1/р1+т2/р2. Средняя плотность смеси тогда рср=(ш 1 +т2)/У=(т 1+т2)р ] р2/(т 1 р2+т2 р 1). Тогда

Рст=т2 р! Р2/(С2(т! р2+т2 ро, что можно далее переписать в виде

Рсш= р1р2/(С1р1+С2р2). Здесь с|=1-с2 - массовая доля (концентрация) первого компонента. Окончательно имеем рст= р1р2/(ср!+(1-с)р2).

Проводя аналогичные рассуждения, получим формулу для вычисления средней вязкости смеси 1дст= Ц1Ц2/(сц1+(1-с) ц 2).

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

Будем далее называть мономер, подлежащий полимеризации, первым компонентом, растворитель с катализатором - вторым. Тогда с = ш2/(ш1+гп2) концентрация второго компонента в единице объема (0 < с £1), а в соответствии с вышеизложенным

р(х,г) =--,м (Х,г)=--

Теперь среднее значение коэффициента молекулярной (ламинарной) кинематической вязкости

уи(х,г)=р(х,гУр(х,г)

Важным для работы является предположение о том, что коэффициент турбулентной диффузии может быть найден из определения числа Шмидта Бс = \iZpDt = у/ТН , представляющего собой отношение импульса, перенесенного за счет вязкости, к импульсу, перенесенному за счет диффузионного массопереноса. Для большинства жидких и газообразных веществ (за исключением нефти и жидких металлов) Бс близко к единице, что позволяет находить коэффицент турбулентной диффузии в виде

= ц, / р Бс = ^ / Бс Здесь и v, коэффициенты вязкости для турбулентного течения.

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

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

1. Движение стационарное, осесимметричное.

2. Компоненты являются несжимаемыми жидкостями или газами.

3. Компоненты - ньютоновские жидкости.

4. Течение происходит без тепло— и массообмена и без химических реакций.

5. Локальные значения скоростей компонентов равны скорости смеси.

6. Компоненты не имеют поверхностей раздела (смешение идеальное)

7. Массовые силы не учитываются.

Для описания движения сплошной среды при сделанных допущениях можно использовать дифференциальные уравнения, выражающие законы сохранения массы и переноса количества движения либо в переменных (\|/, со), либо в простейших физических переменных (u, v, р). При численном моделировании турбулентного смешения часто применение переменных (u,v, р) оказывается более предпочтительным. Это отчасти связано с неудобством использования функции тока \\г при постановке граничных условий в аппаратах с изменяемыми условиями ввода растворителя, содержащего катализатор. Недостаток (vjr, со) систем - невысокая точность локального выполнения уравнения неразрывности при точном выполнении условия заданного расхода через канал. Это связано с тем, что скорости u, v вычисляются как производные от а при конечном шаге интегрирования ошибка в вычислении производных неизбежна, несмотря на то, что при бесконечно малых шагах использование \jr обеспечивает точное выполнение уравнения неразрывности.

В удобном для использования виде система уравнений движения бинарной смеси запишется следующим образом:

(/?ur)x+(pvr)r=0 (puu)x+(pvu)r+ Рх=(( Цц+ й ) ux)x + (1/г) (г (цц + ц,) ur)r— puv/r (puv)x+(pw)r+ Рг=(( Цц + ц,) Vx)x + (1/г) (г (цц + ц,) vr )r— pvv/r -( + ц, )/r2 (puc)x+( pvc)r=(Dcx)x+(Dcr)/r

Система незамкнута из-за наличия коэффициента вязкости v=/i/p и коэффициента диффузии D . Так как для однофазных течений жидкостей при близких значениях вязкости число Шмидта Sc, как это уже упоминалось, близко к единице, то

Здесь V = Уц + Ут - эффективный коэффициент кинематической вязкости, молекулярная (ламинарная), а Ут - турбулентная его составляющие, зависящие от плотности смеси р(х,г), уд - определяется физико-химическими свойствами компонентов.

Для нахождения Ут используются модели турбулентности, которых разработано достаточно большое количество. Одна из наиболее известных работ, связанных с сравнением эффективности современных моделей турбулентной вихревой вязкости — работа [Меп1ег Р.Я.]. В ней выполнена оценка четырех современных моделей. Показано, что ни одна из рассмотренных моделей не дает удовлетворительных результатов для присоединяющихся течений. Там же содержится замечание о том, что для развитых сдвиговых течений стандартная к-е модель дает вполне приемлемые результаты.

В связи с этим, здесь используется стандартная к-е модель турбулентности, уравнения которой для несжимаемой жидкости записаны в следующей форме (р = const)

(urK)x +(vr • K)r =(v-rKx)x +(vrКД +rSK

(urf)x +(vrf)r = - K2

re.

r-er

W =■

F„ =2

1 f\\2>

£ £ = cei ~' ^т^К ~ Ce2

Sk=vtFK-£

(TS=L3, cd = 0.09, cei = 1.44, cf2=1.92 Для течений с переменной плотностью: (purK)x + (¿>vrK)r = (//crr • Kx )x + (//cfrKr)r + pr ■ SK

(pure)x +(/?vre)r=\ цц We, + //„ rrr

+pr Se

crr J ) Ц (7, Pt=Vt P

В уравнениях используются относительные величины:

/> = —, г = -

U V Р D

u=—, v =—, P = — , P =

U„

Pouo '

^"Т*. г=-

и010 "о "о 10 10

Верхний индекс «~» - признак действительных (не относительных) значений, р0, и0 - плотность и скорость на входе, 10 = с1 (диаметр входного сечения аппарата).

Граничные условия

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

-о---

^--т . -е-|

° и.

/-V-

I i 'к

-Т '?L

I I

Рис.3. Типичные формы каналов МТТА. а - локальный радиальный или спутный подвод компонентов; б - зонная подача второго компонента

На входе в канал (сечение О - А, рис.2.2): и(0,г) = и(и(г); у(0,г) = 0; с(0,г) = 0 (0 < г < <1/2),

если используется схема подвода А, то и(0,г) = иог(г); v(0,r) = 0; с(0,г) = 1 (&\И < г < (У2)

Если не задано распределение К(0,г) и е(0,г), то они полагаются равными нулю.

На стенках канала - условия прилипания: и(х,гст) = 0; у(х,гст) = 0. Для давления и концентрации используется условие равенства нулю

нормальной производной: =0.Для К и е: на стенке К(х, гст)

= 0, £г(х,гй6) = V-—rJ на пристеночном слое разностной сетки - их значения в on

соответствии с «методом пристеночных функций».

На оси симметрии; v(x,0) = 0; производные ur, рг, Кг, е,., сг равны нулю.

В выходном сечении (С - Д, рис.2.2) задаются мягкие граничные условия вида fxx = 0.

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

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

Здесь не приводится подробное описание алгоритма SIMPLE. Отметим лишь его особенность, состоящую в использовании «шахматных» сеток. Введение

Э2К

этих сеток связано с необходимостью не допустить появления нефизичного волнообразного поля параметров u(x,r),v(x,r),p(x,r) при решении уравнений движения. Это объясняется тем, что в SIMPLE приходится решать уравнение Пуассона для поправки к давлению с использованием центрально-разностной схемы при аппроксимации вторых производных на сетке, показанной на рис.4. Видно, что значения давления в центральном узле сеточного шаблона, обозначенном квадратиком, определяется через узлы, обозначенные кружочками, и наоборот.

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

□

о-

ii+l О

A,J

□

t 7

О--O^--□

П П О

оооооооо П О————О ГЗ

П.1 и 1+1Л

Рис. 4. Возникновение взаимно-проникающих сеток при численном решении уравнения Пуассона для поправки к давлению

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

На рис. 5,6 показано сравнение распределения составляющей скорости в направлении оси ОХ и кинетической энергии турбулентности к (х,г) в поперечном сечении канала с параллельными стенками с внезапным расширением. Показаны результаты, полученные в ЦИАМ и полученные в наших расчетах (Re = 2,17-10'). Заметим, что на гладкой стенке канала (верхняя

Ък

на рисунке) в программе CANAL использовано условие симметрии: — = 0,

Эг

— = 0. Эти условия являются естественными в такой ситуации. В связи с этим Эг

на стенке канала получены конечные значения для к иг. Сплошной линией показаны результаты Майоровой А.И., пунктирной - наши.

Рис.5. Сравнение распределения составляющей скорости в направлении оси ОХ в поперечном сечении канала с параллельными стенками с внезапным расширением

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

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

На рис.7 приведены сравнительные данные по результатам расчета. Кружочки — экспериментальные данные Chatarvedi М.С. , пунктирная линия — расчет Turgeon Е., Pelletier D., Jgnat L., сплошная линия - расчет Р. Тахавутдинова, А. Мухаметзяновой, крестики — наши вычисления.

Рис. 7. Зависимость и 2 от ПРИ угле образующей диффузора-

конфузора а=90° ; сечение расположено на расстоянии одного калибра от входа

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

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

Число секций от 1 до 6: 1 5,1 ^^^3.0,10° <£ог£90°,0,01 ^ Re^lO7,

(re-m).

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

50 <i Ni £ 500

Размеры разностной сетки:

Точность вычислений соответствует:

- 0,01% по сохранению расхода (ошибка аппроксимации уравнения неразрывности);

- п0 уравнениям переноса (использована разностная схема I12

порядка аппроксимации «против потока»). Это соответствует 2+6 %.

Итоги тестирования позволили сделать вывод о возможности использования комплекса CANAL для решения поставленных в работе задач. Удовлетворительные результаты получались в диапазоне изменения отношений массовых расходов компонентов m2/mi от 0 до 4, отношения плотностей компонентов р2/ pi от 0 до 5, коэффициентов динамической вязкости цг !цх от

0.1 до 10, чисел Шмидта от 0.1 до 10. Следует отметить,что время решения задачи о течении в «эталонном» аппарате (без смешения) Ni *Nj=8000: FLUENT (Кульментьева Е.И.,КГТУ) -180 с, CANAL - 75 с.

Четвертая глава содержит результаты решения следующих задач:

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

б) оценить влияние геометрической формы элементов МТТА на качество смешения;

в) исследовать влияние на качество смешения физических свойств компонентов и режимов течения;

г) оценить потери полного давления в различных условиях работы МТТА.

В качестве критерия эффективности принято распределение коэффициента смешения уСм (х) по длине рабочей зоны аппарата

yCm(x)=l-(l/S(x)) J| 0(x,r)-C(x,r)|dS

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

х 2

У=озЬ> J Г an (ж) А

Хг начало х2 - конец рабочей зоны аппарата. Падение полного давления

Я XV 2

Apo^mlrci J ( p+pw2/2)oA г dr +(m2/mi) J(p+pw2/2)vdx-

0 Jtvl

- J ( p+pw2/2)cD r dr ) 0

mi- расход основного (первого) компонента (условно — воды)

В качестве «эталонного» для сравнения выбран аппарат (рис. 7а) с следующими соотношениями геометрических размеров (по рекомендациям работ Тахавутдинова Р.Г., Мухаметзяновой А.Г.) L/D= 1,7; D/d-1,6; ai=45°; a2=45°

Первый (основной) компонент

Рис. 7а

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

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

с(х,г) у=0.901

12 3 4

Подвод в с в* ими №

Рис. 8. Распределение концентрации и линии тока в эталонном МТТА и сравнение эффективности при подводе компонента в разных его секциях

Таблица 1

формы образующей МТТА в. окрестности горлоЕин.

_ Тип ♦орла Прнэич »и в» -

А > Ч/ м 0-25Ь<0.5

В ь \м / д 1 £ 0-2 <Ь <0.5 0<1<0.5Ь

Е 0<1<Ь

Вычисления показали, что МТТА с профилем типа «Е» значительно более эффективны по сравнению с другими. По этой причине далее исследования проводились именно на таких аппаратах. С двумя секциями предварительной турбулизации (\ZKvx )ио = 0,25 (Рис. 8а).

С(х.г)

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

1. Радиальный подвод второго компонента эффективнее спутного.

2. Аппарат «эталонной» формы не является наиболее эффективным.

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

4. Увеличение числа секций для предварительной турбулизации более двух слабо сказывается на эффективности МТТА.

5. Наибольшей эффективностью обладают МТТА с профилем образующей турбулизатора типа «Е»

Влияние физических свойств компонентов

1. Отношение плотностей компонентов

Уже говорилось о возможности значительного различия плотностей и вязкости компонентов. Нами была проведена серия расчетов смешения в МТТА при изменении отношения плотностей компонентов от р2/р1 от 0.25 до 2.5. Число Рейнольдса, рассчитанное по значению коэффициента вязкости на входе в сечении ОА, принималось равным 1 ООО. Коэффициенты вязкости компонентов принимались одинаковыми. Здесь и далее при проведении вычислений предполагалось отсутствие предварительной турбулизации потока перед входом в МТТА.

у=0.843 л?о^-0.058

р(х.г )

р2/р 1=0.25

Рис.9. «Эталонный» МТТА без предварительной турбулизации. Ду0 - приращение коэффициента смешения по сравнению с «Эталонным» при рг/р1=1-

Рис.10. Смешение в МТТА с предварительной турбулизацией при р2/рг = 0.25

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

2. Отношение коэффициентов вязкости (р2/р1=1,Ш2/гп)=1,8с=1) В этой серии численных экспериментов плотности обоих компонентов принимались равными. Коэффициент кинематической вязкости первого (основного) компонента принимался постоянным, а отношение коэффициента кинематической вязкости второго к первому варьировалось в пределах от 0.05 до 10. Результаты показаны на рис.11. Видно, что качество смешения существенно ухудшается с ростом вязкости второго компонента. Однако, как правило, второй компонент (растворитель) имеет меньшую, по сравнению с основным, вязкость. Поэтому в рассмотренной схеме аппарата влияние вязкости несущественно, хотя в других схемах, где более вязкий компонент вводится в периферийную часть потока, соприкасающуюся со стенкой, это влияние может быть значительным.

09 0.8 0.7 Ё

Ау=—0.013

е--------=------:------ |у=0.912_—

С(х,г)

ц2/ц 1=0.5

С(х,г)

Ау=—0.063

ц2/ц 1=2.0

Рис.11. Смешение в МТТА с предварительной турбулизацией при различных отношениях коэффициентов динамической вязкости

Результаты исследований (при ш2/ш1=1, Ук =сош1 => Бух -переменна) приводятся в виде диаграммы. Как это видно из (рис. 11а), диаграмма имеет вид седловой поверхности с явно выраженной точкой ШШМАХ-а с координатами (2.0; 1.0).

□ 0,930-0,940

■ 0,920-0,930

■ 0,910-0,920

□ 0,900-0.910

■ 0,890-0,900 и 0,830-0,890

■ 0,870-0.880

□ 0,860-0,870

□ 0.850-0.860

■ 0,840-0,850

■ 0,830-0,840

Рис. 11а

3. Число Шмидта ( Бс = V, / Б, )

Б, - коэффициент турбулентной диффузии V, - коэффициент турбулентной вязкости

Ранее было отмечено, что при проведении численных экспериментов предполагалось, что Бс- число Шмидта для рассматриваемых компонентов близко к единице, и поэтому все приведенные выше результаты получены при 8с=1. Как это известно, число Шмидта определяется соотношением между интенсивностью процессов переноса за счет вязкости и за счет диффузии. При выполнении работы были проведены вычисления с измененными его значениями от Бс^ОЛ до Бс^Ю.О с целью определить чувствительность распределения параметров течения по отношению к этим изменениям. Как это видно на рис.12 влияние числа Шмидта на характер поля концентраций весьма существенно.

С(х,г)

вс=1и.и / *

Вынос второго компонента вверх по потоку

8с=0.1

'.О 4» . !. О Г>в «.О ел ?,0 Г.» го М ».О О 4 10 0

Рис. 12. Поля концентраций и зависимость у — Г (Бс) для аппарата с предварительной турбулизацией

Выводы:

1. В «эталонном» аппарате уменьшение плотности второго компонента сильно ухудшает качество смешения (у).

2. В аппаратах с предварительной турбулизацией уменьшение плотности

второго компонента незначительно изменяет у

3. Уменьшение плотности второго компонента вызывает вынос его вверх по потоку от места ввода.

4. С увеличением вязкости качество смешения снижается.

5. При повышенной вязкости второго компонента происходит вынос его на

стенку за местом его ввода.

6. Изменение числа Шмидта сильно влияет на качество смешения.

4. Влияние режима течения Влияние числа Ие

Решалась задач выяснения влияния режима течения (ламинарный, переходный, турбулентный при малых, средних и и при больших числах Рейнольдса (Яе= 103 -106).Число И.е вычислялось по условиям на входе в канал МТТА (сечение ОА). Вычисления проведены при т2=Ш1, р2= р1, у2= VI Яе=и0 с!/V

1Ло - скорость на входе в аппарат, с! - диаметр входа, V - вязкость первого компонента (сечение О-А)

Рис. 12а

Отношение массовых расходов компонентов (ш2/ш1-8ух=сопз1, Ук=переменна)

переменно,

у=0.939

С(х,г)

1712/1111=4

у==0.836

Ау=+0.014

Лу=—0.089

Рис. 13. Качество смешения в зависимости от соотношения расходов компонентов(при постоянной скорости подвода второго компонента Ук=сопБ1)

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

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

Коэффициент смешения

■ 0,950-1,000

□ 0,900-0,950

□ 0,850-0,900

■ 0,800-0,850

■ 0.750-0,800

Ук/ио

Рис. 14

Падение полного давления

й XVI

Ар0 =2 71 ( | ( р+р\у2/2)ОА г с1г +(т2/т,) Я |(р+р\у2/2)у ёх -

0 XV)

- \ ( р+р\у2/2)со г ёг ) Ш1/ (гп1+ш2)

Шр расход основного (первого) компонента (условно — воды)

т2)'ш1

Рис. 15

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

1. Проанализировано состояние . исследований гидродинамических процессов в технических устройствах, близких по характеру протекающих в них процессов к малогабаритному трубчатому реактору диффузор -конфузорного типа (МТТА). Сделан вывод о том, что систематических исследований влияния формы образующей' канала рабочей части. на его характеристики нет;

2. Разработан алгоритм решения внутренней задачи о турбулентном течении вязкой несжимаемой жидкости (или газа) и смешении компонентов в малогабаритных трубчатых аппаратах (МТТА) на основе системы уравнений гидродинамики, уравнения диффузии и уравнений стандартной к -е модели.

3. Разработан специализированный программный комплекс программ CANAL для ЭВМ, позволяющий моделировать движение, потока двухкомпонентной смеси в каналах МТТА.

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

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

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

Основное содержание диссертации опубликовано в работах

1. Петровичева Е.А. Использование реакционно-массообменных процессов в химической промышленности. /Макарова Т.П., Ильина И.М.// Материалы научно-технической конференции «АлНИ-2002».-Альметьевск, 2003. С. 149.

2. Петровичева Е.А. Численное моделирование турбулентных течений в ТРДКТ. /Ильина И.М.// Научная сессия. — г. Альметьевск, 2004 г. С. 44.

3., Петровичева Е.А. Турбулентное смешение компонентов в трубчатых аппаратах. /Данилов Ю.М., ■ Мухаметзянова А.Г.// Химическая промышленность. 2004. Т. 81, №9. С. 451-457.

4. Петровичева Е.А. Математическое моделирование гидродинамики смешения компонентов в трубчатых турбулентных аппаратах. /Данилов Ю.М.// Труды Международной научной конференции «Математические методы в технике и технологиях», Казань: 31 мая-2 июня 2005, С. 312-313. Казанский государственный технологический университет.

5.Петровичева Е.А. Расчет осесимметричных турбулентных течений вязкой несжимаемой жидкости на основе программы CANAL 4S (5S). /Ильина -И.М.// Ученые записки. Сборник научных трудов. Т.З, Альметьевск: АГНИ, 2005. С. 236-258.

6. Петровичева Е.А. Проблемы смешения однофазных компонентов в трубчатых турбулентных аппаратах. /Макарова Т.П., Ильина И.М.// Материалы научной сессии по итогам 2004 г., Альметьевск: АГНИ, 2005.

7. Петровичева Е.А. Методы численного решения задач для турбулентных течений смешивающихся жидкостей в МТТА. /Корепанова Л.Ф., Ильина И.М.// Ученые записки. Сборник трудов, Альметьевск: АГНИ, .

2005, С. 119-132.

8. Петровичева Е.А. Двухуровневая математическая модель процесса сополимеризации в газожидкостном каталитическом реакторе. /Тахавутдинов Р.Г., Мухаметзянова А.Г., Кульметьева Е.И.// Труды Международной научной конференции «Математические методы в технике и технологиях», Казань, 2005 г. Казанский государственный технологический университет.

9. Петровичева Е.А. Численное моделирование турбулентных течений с помощью пакета PHOENICS. /Макарова Т.П., Ильина И.М.// Материалы научной сессии по итогам 2005 г., Альметьевск: АГНИ, 2006.

10. Петровичева Е.А. Турбулизация потока в трубчатых аппаратах как способ повышения эффективности смешения компонентов. /Ильина И.М., Макарова Т.П.// Ученые записки. Сборник трудов. Том 4. Альметьевск: АГНИ,

2006, С. 225-232.

Соискатель

Е.А. Петровичева

Подписано в печать 27.09.2006 г.

Формат 60x84/16 Печать RISO 1,75 уч.-изд.л..' 1,8 ус. печ .л. Тираж 100 экз. Заказ №77

ТИПОГРАФИЯ АЛЬМЕТЬЕВСКОГО ГОСУДАРСТВЕННОГО

НЕФТЯНОГО ИНСТИТУТА 423452, Татарстан, г. Альметьевск, ул. Ленина, 2

Содержание диссертации автор исследовательской работы: кандидата технических наук, Петровичева, Елена Александровна

Введение.

Глава 1. Турбулизация потока в трубчатых аппаратах как способ повышения эффективности смешения компонентов

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

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

1.3. Схемы турбулизаторов и их эффективность.

1.4. Турбулизация потока в каналах типа диффузор-конфузор.

1.5. Роль отрывной зоны в турбулизации потока в трубах.

1.6. Сравнение экспериментального, теоретического и численного подходов для решения задач гидродинамики.'.

1.7. Численное моделирование турбулентных течений с помощью коммерческих пакетов программ.

1.8. Методы численного решения задач математической физики.

1.9. Численное исследование течений в трубчатых турбулентных реакторах диффузор-конфузорного типа.

Глава 2. Система уравнений движения бинарной смеси.

2.1. Основные упрощающие допущения.

2.2. Система уравнений.

2.3. Метод пристеночных функций (закон стенки).

Глава 3. Численное решение задач о турбулентном течении вязкой несжимаемой жидкости в осесимметричных каналах.

3.1. Выбор метода решения краевых задач.

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

3.3. Описание алгоритма и программы CANAL.

3.4. Тестирование и настройка программного комплекса CANAL.

3.5. Тестирование программы.'.;.

Глава 4. Исследование влияния способов подвода компонентов, геометрической формы элементов МТТА, физических свойств компонентов и режимов движения на качество смешения.

4.1. Влияние способа организации смешения и формы турбулизаторов.

4.2. Многозонный подвод второго компонента.(подвод второго компонента через центрально-расположенную трубку).

4.3. Влияние профиля образующей в окрестности горловины.

4.4. Влияние физических свойств компонентов.

4.4.1. Отношение плотностей компонентов.

4.4.2. Отношение коэффициентов вязкости.

4.4.3. Число Шмидта.,.

4.4.4. Влияние режима течения.

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

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

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

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

3. Программа Республики Татарстан по развитию приоритетных направлений науки по теме № 19-12/99 (Ф) «Научные основы технологических процессов производства синтетических каучуков на предприятиях нефтехимического комплекса Республики Татарстан». 4. Программа Республики Татарстан по развитию приоритетных направлений науки по теме № 07-7.5 - 27/2001 (Ф) «Ресурсосберегающие и экологически безопасные трубчатые аппараты и технологические процессы для нефтехимической промышленности». Целями работы являются:

1. Разработать методику численного решения уравнений гидродинамики для расчета осесимметричных турбулентных течений бинарных смесей несжимаемых жидкостей (или газов) в МТТА на основе наиболее надежных, известных к настоящему времени, численных методов.

2. Произвести математическое моделирование турбулентных течений смесей несжимаемых жидкостей в малогабаритных трубчатых аппаратах диффузор -конфузорного типа (МТТА) для выяснения их основных особенностей.

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

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

1. Разработан эффективный численный алгоритм решения уравнений гидродинамики, бинарных смесей жидкостей или газов в трубчатых аппаратах, (в частности, в длинных трубах L/D » 1).

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

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

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

Результаты, полученные при выполнении диссертационной работы, позволяют:

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

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

• Оценить потери энергии на организацию смешения.

• Разработанный программный комплекс CANAL для ПЭВМ, ориентированный на расчет турбулентных течений в МТТА, может быть использован для определения характеристик течения и, в силу незначительных затрат времени на решение и простоты обращения с ним может рассматриваться как инструмент для инженерных расчетов наряду с другими известными инженерными методами.

Апробация работы. Материалы исследования и его результаты неоднократно обсуждались на заседаниях кафедры прикладной химии Альметьевского. государственного нефтяного института и докладывались на научных конференциях: на научно-технической конференции « АлНИ-2000», г. Альметьевск, 2001г; на научной сессии КГТУ им. С.М. Кирова, г.Казань, 2001г.; на III Международной научно-практической конференции ПГУ, г. Пенза, 2001г.; на научно-технической'конференции « АлНИ- 2002», г. Альметьевск, 2003г.; на Всероссийской научно-технической конференции «Большая нефть: реалии, проблемы, перспективы», г. Альметьевск, 2001г.; в Альметьевском государственном нефтяном институте, 2003, 2004, 2005 гг.; в Казанском государственном техническом университете им. Туполева, 2003г.; на 4-й Школе-семинаре молодых ученых и специалистов под руководством акад: В.Е. Алемасова, г.Казань, 2004 г.; на научно-технических семинарах в Казанском государственном технологическом университете, 2003, 2004, 2005 гг.; на Международной научной конференции « Математические методы в технике и технологиях-ММТТ-18», г. Казань, 2005 г.

По материалам диссертационной работы опубликовано 10 работ.

1. Петровичева Е.А. Использование реакционно-массообменных процессов в химической' промышленности. /Макарова Т.П., Ильина И.М.// Материалы научно-технической конференции «АлНИ-2002».-Альметьевск, 2003. С. 149.

2. Петровичева Е.А. Численное моделирование турбулентных течений в ТРДКТ. /Ильина И.М.// Научная сессия. - г. Альметьевск, 2004 г. С. 44.

3. Петровичева Е.А. Турбулентное смешение компонентов в трубчатых аппаратах. /Данилов Ю.М., Мухаметзянова А.Г.// Химическая промышленность. 2005. Т. 81, №3. С. 152-157.

4. Петровичева Е.А. Математическое моделирование гидродинамики смешения компонентов в трубчатых турбулентных аппаратах. /Данилов Ю.М.// Труды Международной научной конференции «Математические методы. в технике и технологиях», Казань: 31 мая-2 июня 2005, С. 312-313. Казанский государственный технологический университет.

5. Петровичева Е.А. Расчет осесимметричных турбулентных течений вязкой несжимаемой жидкости на основе программы CANAL 4S (58)./Ильина И.М.// Ученые записки. Сборник научных трудов. Т.З, Альметьевск: АГНИ, 2005. С. 236-258.

8. Петровичева Е.А. Двухуровневая математическая модель процесса сополимеризации в газожидкостном каталитическом реакторе. /Тахавутдинов Р.Г., Мухаметзянова А.Г., Кульметьева Е.И.// Труды Международной научной конференции «Математические методы в технике и технологиях», Казань,2005 г. Казанский государственный технологический университет.

10. Петровичева Е.А. Турбулизация потока в трубчатых аппаратах как способ повышения эффективности смешения компонентов. /Ильина И.М.,

Макарова Т.П.// Ученые записки. Сборник трудов. Том 4. Альметьевск: АГНИ, 2006, С. 225-232. .

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

Заключение диссертации по теме "Механика жидкости, газа и плазмы"

Заключение

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

1. Проанализировано состояние . исследований гидродинамических. процессов в технических устройствах, близких по характеру протекающих в них процессов к малогабаритному трубчатому реактору диффузор -конфузорного типа (МТТА). Сделан вывод о том, что систематических исследований влияния формы образующей канала рабочей части на его характеристики нет.

3. Разработан специализированный программный комплекс программ CANAL для ЭВМ, позволяющий моделировать движение потока двухкомпонентной смеси в каналах МТТА.

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

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

Список источников диссертации и автореферата по механике, кандидата технических наук, Петровичева, Елена Александровна, Казань

1. Берлин А,А., Минскер К.С., Дюмаев К.Н. Новые унифицированные энерго - и ресурсосберегающие высокопроизводительные технологии повышенной экологической чистоты на основе трубчатых турбулентных реакторов. М.: ОАО «НИИТЭХИМ», 1996.-188 с.

2. Минскер К.С., Берлин А.А., Захаров В.П., Дьяконов Г.С., Мухаметзянова А.Г., Заиков Г.Е. Быстрые процессы при синтезе полимеров 7/ Журнал прикладной химии. 2003.Т.76. №2. С.272-278.

3. Захаров В.П., Мухаметзянова А.Г., Тахавутдинов Р.Г., Дьяконов Г.С., Минскер К.С. Создание однородных эмульсий в трубчатых турбулентных аппаратах диффузор конфузорной конструкции // Журн.прикл.химии. 2002.Т.75.№9. С. 1462-1465.

4. Тахавутдинов Р.Г., Мухаметзянова А.Г., Дьяконов Г.С., Минскер К.С., Берлин А.А.Трубчатые турбулентные предреакторы для проведения процессов инициирования при каталитическом синтезе полимеров// Высокомолекулярные соединения. 2002. Т.44. №7. С. 1094 1100.

5. Минскер К.С. , Дьяконов Г.С., Тахавутдинов Р.Г., Мухаметзянова А.Г., Захаров В.П., Берлин А.А. Многофазные течения в трубчатых аппаратах диффузор конфузорной конструкции // Доклады РАН. 2002. Т.382. №4. С. 509-512.пб

6. Тахавутдинов Р.Г. Теоретические основы и методы повышения эффективности промышленных полимеризаторов в производстве синтетических каучуков. Дисс. на соиск. уч. степ, доктора технических наук. Казань. 2000, 380 с.

7. Мухаметзянова А.Г. Малогабаритные трубчатые аппараты в производстве синтетического каучука СКЭП (Т). Дисс. на соиск. уч. степ, кандидата технических наук. Казань,2002. 149 с.

8. Тахавутдинов Р.Г., Дьяконов Г.С., Мухаметзянова А.Г. Интенсификация диспергирования в трубчатых турбулентных аппаратах при производстве синтетических каучуков//Химическая промышленность, №1, 2000, с. 1-6.

9. Лойцянский J1.T. Механика жидкости и газа. М.: «Наука», 1987, 840 с.

10. Мигай В.К. Повышение эффективности современных теплообменников. Л.: Энергия. Ленинград, отделение, 1980, 144 с.

11. Калинин Э.К., Дрейцер Г.А., Ярхо С.А. Интенсификация теплообмена в каналах. Издание 3-е, перераб. и дополн. М.: Машиностроение, 1990, 208 с.

12. Ван-Дайк М. Альбом течений жидкости и газа. М. «Мир». 1986, 182 с. (Wan Dyke М. Au , Album of Fluid Motion. The Parabolic Press, Stanford, California, 1982.)

13. Runchal A.K. 1969 International Seminar "Heat and Mass Transfer in Flow with Separated Regions and Measurement Techniques", Herceg-Novi, Yugoslavia, September 1-13, 1969.

14. Майорова А.И., Ягодкин В.И. Методика и результаты расчетов течений в каналах с внезапным расширением. М. Труды ЦИАМ, №883, 1979, 31с.

15. Chatarvedi М.С. Flow characteristics of axusimmetruc exspansions // Yonrnal of the Hydraulics Division. Proceedings of the American Sosiety of Cuvil Engineers, 1963, v.89, p.61.

16. Turgeon E., Pelletier D., Jgnat L. Effects of adaptivity on varions finite element schemes for turbulent heat transfer and flow predictions// Papers of of the 36

17. Самарский А.А. Теория разностных схем. М.: Наука. ГРФМЛ, 1983, 616 с.

18. Флетчер К. Численные методы на основе метода Галёркина. М.: «Мир», 1988,352 с.

19. Mitchell A.R. Varitional principles and the finite elements method //J.instit.Math. and Its Appl.1972, vol 9, p. 378-389.

20. Бреббия К.И. и др. Методы граничных элементов М.: Мир, 1987, 520 с.

21. Григорьев М.М., Фафурин А.В. Метод граничных элементов для течений жидкости при высоких числах Рейнольдса. ИВУЗ «Авиац.техника», 1996, №4, с.33- 42.

22. Данилов Ю.М. Инвариантные операторы для численного решения краевых задач гидродинамики. Казань: КГТУ, 1999, 140 с.

23. Данилов Ю.М. Численное решение стационарных уравнений гидродинамики в дозвуковй области течения.// ИВУЗ. Авиац.техника. 1980, №3, с.42 45.

24. Orlandi P.Vorticity — velocity formulation for High re Flows// Computers & Fluids, vol.15. №2, 1987, pp. 137-149.

25. Кочин H.E., Кибель И.А., Розе H.B. Теоретическая гидромеханика. ч.1, И. М.: ГИФМЛ, 1963.

26. Рычков А.Д. Математическое моделирование газодинамических процессов в каналах и соплах. Н ск.: «Наука» Сибирское отделение АН СССР, 1988, 221 с.

27. Белов И.А., Исаев С.А. Моделирование турбулентных течений. Учебное пособие. СПб. гос. техн. ун-т. 2001,108 с.

28. Launder В.Е., Spalding D.B., The numerical computation of turbulent flow // сотр. Math. Appl. Mech. Eng. 1974. vol.3, №2, p.269-289.

29. Smith L.M., Reynolds W. С. On the Yakhot Orszag Renormalization group method for driving turbulence statistic and models. Phys. Fluids. A, 1992, v.4, m, p.364.

30. Yahot V., Orszag S.A. Renormalization group analysis of turbulence. J.Sci.1. Comput. 1986, v.l,p.3.

31. Роуч П. Вычислительная гидродинамика. М.: Мир, 1980, 616 с.

32. Данилов Ю.М., Дегтярёва О.М., Хасанов Р.Х. Расчет газодинамики и теплообмена в осесимметричных каналах сложной геометрической формы. // Межвузовский сборник «Оптимальные задачи авиационной техники» Казань, КАИ, 1990. С. 105.

33. Томпсон. Методы расчета сеток в вычислительной гидродинамике // Аэрокосмическая техника. 1985, №8, с. 141 171.

34. Теория пограничного слоя / Под ред. Л.Г. Лойцянского М.: Наука, 1974, 711с.

35. Launder В.Е., Spalding D.B. Mathematical Models of Turbulence. London: Acad. Press, 1972.

36. Фафурин B.A. Расчет вращающегося турбулентного потока на основе модифицированной к-е модели. Тепломассообменные процессы и аппараты химической технологии. Межвузовск. тематическ. сб. научн.трудов. Казань: КГТУ, 2000, с.27 32.

37. Данилов Ю.М., Ильина И.М., Ситдикова И.П. Численное решение трехмерных задач гидродинамики. Материалы научно-технической конференции « АлНИ- 2002». -Альметьевск, 2003.-147с.

38. Макарова Т.П., Ильина И.М., Петровичева Е.А. Использование реакционно-массообменных процессов в химической промышленности. Материалы научно-технической конференции « АлНИ- 2002».-г.Альметьевск, 2003.-149с.

39. Данилов Ю.М., Ильина И.М., Ситдикова И.П., Бергман А.Н.Решение трехмерных задач о течении вязких несжимаемых жидкостей в каналах прямоугольной формы. «Естественные и технические науки»,Ж.№3, 2003г.

40. Данилов Ю.М., Дьяконов Г.С., Мухаметзянова А.Г., Бергман А.Н., Ильина И.М. Оптимизация проточной части трубчатых турбулентных реакторов.ВестникКГТУ., г.Казань, 2003г.

41. Шамсутдинов A.M., Хапугин И.Н., Махоткин А.Ф., Ильина И.М., Г.В. Каштанова. Исследование массоотдачи в газовой фазе в вихревых контактных устройствах.

42. Ильина И.М., Шамсутдинов A.M. Тепловая устойчивость закрученных потоков и разработка теплообменного аппарата. Материалы научно-технической конференции « АлНИ-2000».-Альметьевск, 2001.-65с.

43. Хапугин И.Н., Шамсутдинов A.M., Махоткин А.Ф., Ильина И.М. Исследование моссоотдачи в газовой фазе в вихревых контактных устройствах. Проблемы энергетики.//Известия высших учебных заведений №11-12.-Казань:Изд-во КГЭУ,2001.-С.11-20.

44. Ильина И.М., Петровичева Е.А. Численное моделирование турбулентных течений в ТРДКТ. Научная сессия. г. Альметьевск, 2004г.

45. Данилов Ю.М., Мухаметзянова А.Г., Дьяконов Г.С., Кульментьева Е.И.

46. Математическое моделирование течений в малогабаритных трубчатых турбулентных аппаратах // Химическая промышленность. 2004. Т.81, №9. С.451-457.

47. Дьяконов Г.С, Данилов Ю.М., Мухаметзянова А.Г., Бергман А.Н., Ильина

48. И.М. Численное моделирование течений в трубчатых аппаратах. Вестник КГТУ, № 1-2, г.Казань, 2002г.-С.267-272.

49. Шамсутдинов A.M., Махоткин А.Ф., Ильина И.М., Шамсутдинов М.А.,

50. Новширванов А.Г., Михайлов Н.М. Разработка абсорбционной вихревой технологии. Научная сессия 2001г. Аннотация сообщений 2000г.

51. Васильцов Э.А., Ушаков В.Г. Аппараты для перемешивания жидких сред. Справочное пособие. JI.: Машиностроение, 1979г.

52. Дрейцер Г.А. Модель процесса солеотложения при обтекании охлаждающей водой труб с кольцевыми турбулизаторами // Современные проблемы гидродинамики и теплообмена в элементах энергетических установок и криогенной технике.- М.:1988.- С. 69 77.

53. Бодров В.И., Дворецкий С.И., Дворецкий Д.С. Оптимальное проектирование энерго- и ресурсосберегающих процессов и аппаратовхимической технологии // Теоретические основы химической технологии. 1997. Т.31. №5. С. 542-548.

54. Берлин А.А.,Прокофьев К.В., Минскер К.С. и др. // Химия и технология топлив и масел. 1988. №7, с.8.

55. Крехова М.Г., Минскер С.К., Прочухан Ю.А. Влияние турбулентности на эффективность смешения потоков разной плотности // Теоретические основы химической технологии, 1994. Т.28. №3,с.271-273.

56. Минскер К.С., Берлин А.А,, Свинухов А.Г.и др. .// Докл. АН СССР. 1986. •Т.286,№5, с. 1171.

57. Берлин А.А., Минскер К.С. и др. // Докл. АН СССР. 1986. Т.287.№1, с. 1.45.

58. Берлин А.А., Минскер К.С., Прокофьев К.В. и др. // Нефтепереработка и нефтехимия. 1988.№2. с.25.

59. Минскер К.С., Прокофьев К.В., Прочухан Ю.Л. и др. // Там же. 1989.№1. с.17.

60. Минскер К.С., Прочухан Ю.А., Колесов С.В. и др. // Там же,1989.№1 l.c.40.

61. Reynolds О. On the dynamical theory of incompressible viscous fluids and the determination of the Criterion // Phi. Trans. Rou. Soc. London. Ser.A. 1884.V.186.P.123-16L

62. Taylor G.I. Production and dissipation of vorticity in a turbulent fluid //Proc.

63. Roy. Soc. Ser/А/ 1938.V.164.P. 15-23.

64. Каминский B.A., Рабинович А.Б., Федоров А.Я., Фрост В.А. Моделирование турбулентного микросмешения. // Теоретические основы химических технологий. 1997.Т.31. №31. С.243.

65. Каминский В.А., Федоров А.И., Фрост В.А. Методы расчетов турбулентных потоков с быстрыми химическими реакциями.// Теоретические основы химической технологии. 1994.Т.28. № 6. С. 591.68. . Каминский В.А., Рабинович А.Б., Федоров А.Я., Фрост В.А.

66. Физикохимия микросмешения в турбулентных потоках с химическими реакциями. //Журнал физической химии. 1995. Т.69. №8. С. 1456.

67. Котов С.В., Берлин А.А., Прокофьев К.В. и др. // Химия и технология топлив и масел. 1990.№ 6. С. 10.

68. Котов С.В., Атманджаков В.Ё., Мйнскер К.С. и др.'// Нефтепереработка и нефтехимия. 1992.№1.С. 38.

69. Minsker K.S., Berlin A.A., Fast Polymerization Process: Gordon and Breach Publ.Inc., 1996.146р.

70. Берлин A.A., Дюмаев K.M., Мйнскер К.С. и др. // Химическая промышленность. 1995.№9.С.550.

71. Берлин А.А., Мйнскер К.С., Дюмаев К.М. и др. Проблемы протекания быстрых химических реакций синтеза низкомолекулярных продуктов в потоке. Новая технология. // Химическая промышленность. 1997.№5.С. 329-332.

72. Колмогоров А.Н. Уравнения турбулентного движения несжимаемой жидкости: //Изв. АН СССР. 1942.Т.6. №1/2.С. 56-58.

73. Колган В.П. Применение принципа минимальных производных к построению конечно-разностных схем для расчета разрывных решений газовой динамики. //Уч.зап.ЦАГИТ.З.№6.1997.С.68-72.

74. Prandtl L., Wlighardt К. Cher eir neues Formelsystem fur die ausgebildete Turbulenz // Nachr. Ges. Wiss. Math.-Phys. K.l. Gottingen. 1945. Bd. 11 A. S.6-19.

75. Лисейкин В.Д. Обзор методов построения структурных адаптивных сеток. // Ж. Вычисл. матем. и матем. физики. Т.36. №1. 1996. С. 3-41.

76. Копылева Б.Б., Павлушенко И.С. Теория и практика перемешивания в жидких средах. М.: НИИ ТЭХми. 1973.

77. Баренблатт Г.О. О движении взвешенных частиц в турбулентном потоке. //UMM. 1953.Т.17.С.261-274.

78. Крамере X, Вестертерн К. Химические реакторы. М.: Химия. 1967.

79. Villermaux J., Falk L. A generalized mixing model for initial contacting of reactive fluids // Chem. Eng. Sci. 1994.V.49.P.5127.

80. Зайчик Л.И. Модели турбулентного переноса импульса и тепла в дисперсной фазе, основанные на уравнении для вторых и третьих моментов пульсаций скорости и температуры частиц. // ИФЖ, 1992.Т.63. С.404-413.

81. Balduga J., Bourne J.R., Pllarn S.J. Interaction between chemical reactions and mixing on varions scales // Chem. Eng. Sci. 1997. V.52. P.457.

82. Balduga J., Bourne J.R. Simplification of micromixing calculations/ Chem. EngngJ. 1989.V.42., p. 83-101.

83. Малкин А.Я., Эпенштейн Г.А., Березовский A.B. и др. Течение полимеризующейся жидкости в трубчатом реакторе. // Теоретические основы химической технологии. Т.20.№43. 1986.

84. Tsai К., Fox R.O. PDF Modeling of turbulentmixing effects on initiator efficiency in a tubular 1ЛЭРЁ reactor. AI Ch E Journal, 1996. Vol. 42.№ 10. P. 2926-2940.

85. Шрайбер A.A., Гавин Л.Б., Наумов B.A., Яценко В.П. Турбулентные течения газовзвеси. Киев, 1987. 240с.

86. Мухаметзянова А.Г., Захаров В.П., Тахавутдинов Р.Г., Дьяконов Г.С., Минскер К.С. Движение многофазных потоков в трубчатых каналах диффузор-конфузорной конструкции // Вестник Башкирского университета. 2002.№1.С.28-31.

87. Мухаметзянова Д.Г., Тахавутдинов Р.Г., Дьяконов Г.С. Численный расчет гидродинамики двухфазного турбулентного потока в диффузор-конфузорном аппарате // Научная сессия КГТУ. Аннотация сообщений.Казань. 2001.С.91.

88. Колмогоров А.Н. ДАН СССР. 1949.Т.66.С. 825.

89. Баранаев М.К., Теверский Е.Н., Трегубова Э.Л.- ДАН СССР. 1949.Т.66. С.272.

90. Патанкар С. Численные методы решения задач теплообмена и динамики жидкости. Пер.с англ.-М.:Энергоатомиздат, 1984. С.152.

91. Патанкар С., Сполдинг Д. Тепло и массообмен в пограничных слоях М.: Энергия. 1971. С. 127.

92. Нигматуллин Р.И. Динамика многофазных сред. М.: «Наука». 1987.

93. Кафаров В.В. Основы массопередачи. М.: Высш.школа. 1979.

94. Calderbank Р.Н. // Trans.Inst. Chem. Engrs. 1959.V.37.№2.P.131.

95. Павлушенко И.С., Брагинский Л.Н., Брыкав В.Н. II Журн.прикл.химии. 1972.Т.6.№5.С.821.

96. Тахавутдинов Р.Г., Дьяконов Г.С., Дебердеев Р.Л., Мухаметзянова А.Г.,

97. Минскер К.С. Турбулентное смешение в малогабаритных трубчатых аппаратах // Тепломассобменные процессы и аппараты химической технологии. Межвуз. тем.сб.науч.трудов. Казань, КГТУ.2000. С. 38-50.

98. Литманс Б.А., Кукуреченко И.С., Бойко И.Д. и др. //Теор. основы хим. технол. 1975.Т.6. № 5. С.821.

99. Размолдин Л.П., Коротков А. Л., Кузмичев Ю.Б. Математическая модельмассообмеиа пузырька пара с жидкостью при наличии градиентов поверхностного натяжения.// Теор. основы хим. технологии. Т.24.№4. 1990.С.570-572.

100. Whitman W.G. // Chem. Met.Engng.l923.V.29.P.146.

101. Nernst N.// Z.Phus.Chem.l904.V.47.P.52.

102. Прочухан Ю.А., Минскер К.С. и др. Влияние способов смешения на характер протекания сверхбыстрых полимеризационных процессов //• Высокмолек.соед. 1988. Т.30.№6.С.1250.

103. Берлин А.А., Коноплев А.А. Минскер К.С. и др. Влияние геометрии течения и способа ввода реагентов на характеристики смешения в проточных реакторах. //Докл. АН. 1989.Т.305.№5. С. 1143.

104. Ландау Л.Д., Лифшиц Е.Н. Теоретическая физика.Т.6. Гидродинамика.М.: Наука. 1986.

105. Брагинский Л.Н., Бегачев В.И., Барабаш В.Н. Перемешивание в жидких средах. Изд.:Химия. 1983. С.255.

106. Касаткин А.Г. Основные процессы и аппараты химической технологии. М.: Химия. 1971.С.784.

107. Кафаров В.В. Основы массопередачи. М.: Высш.шк.1979.

108. Литманс Б.А., Кукуреченко И.С., Туманов Ю.В. Исследование массоотдачипри больших затратах удельной мощности на перемешивание // Теория и практика перемешивания в жидких средах. М.: НИИТЭХИМ, 1973.С. 137.

109. Акберов P.P., Понявин В.И., Фафурин В.А. Численное моделирование течений в осесимметричных каналах методом конечных элементов // Тепломассобменные процессы и аппараты химической технологии. Межвуз.тем.сб.научн.трудов. Казань. КГТУ.1998.С.160-167.

110. Дьяконов Г.С., Тахавутдинов Р.Г., Курочкин И.М. Мухаметзянова А.Г. Влияние диффузионных сопротивлений на перенос мономеров умежфазной поверхности. // Теоретические основы химических технологий. 2000.Т.34.№1 .С.82-86.

111. Дьяконов Г.С., Тахавутдинов Р.Г., Курочкин И;М., Мухаметзянова А.Г.

112. Тепломассоперенос при химических превращениях в полимеризаторе синтеза этилен-пропиленового каучука. //2-ая Всероссийская национальная конференция по теплообмену. Москва. 1998.Т.З.С.198-201.

113. Хинце И.О. Турбулентность, ее механизм и теория. -М.: Физматгиз, 1963.

114. Шорин С.Н. Теплопередача. -М.: Высшая школа, 1964.

115. Шлихтинг Г: Теория-пограничного слоя. М.: Наука, 1969.

116. Брэдшоу П. Введение в турбулентность и ее измерение. М.: Мир, 1974. -280с.

117. Варфоломеев И.М., Глебов Г.А., Гортышов Ю.Ф. и др. Структура турбулентного отрывного течения в прямоугольной полости // Пристенные струйные потоки. Новосибирск, 1984. С. 86-92.

118. Гарсиа, Сперроу. Турбулентный теплообмен за участком резкого сужения канала типа обращенной вперед ступеньки // Теплопередача. 1988.- №2. -С.60.

119. Минскер С.К., Голубева Т.В., Коноплев А.А. и др. Формирование плоскогофронта реакции при проведении весьма быстрых химических процессов в турбулентных потоках // Доклады АН. 1990.Т.ЗГ4,№6.С. 1450-1454.

120. Минскер С.К., Коноплев А.А., Минскер К.С., Прочухан Ю.А.Домпаниец В.З., Берлин А.А. Организация фронта реакции в турбулентном потоке // Теорет. Основы хим.технолог. 1992.Т.1992, № 5.С. 686-691.

121. Патанкар С. Численные методы решения задач теплообмена и динамики жидкости: Пер. с англ. М.: Энергоатомиздат, 1984. 152 с.

122. Ахвердиев Р.Ф. Турбулентное движение продуктов горения в пиролизных установках. Кандидатская диссертация. Казань, КГТУ, 2004, 147 с.

123. Ильина И.М. Повышение эффективности трубчатых аппаратов на основечисленного моделирования турбулентных течений в их пробочной части, кандидатская диссертация. Казань, КГТУ, 2004 г.

124. F.R. Menter. Сравнение некоторых моделей турбулентной вихревой вязкости A Comparison of somt Recent Eddy-Viscosity Turbylence Models. Trans. ASMEJ. Fluids End. 1996-118. №3-c. 514-519.-Англ.

125. Петровичева Е.А. Турбулентное смешение компонентов в трубчатых аппаратах. /Данилов Ю.М., Мухаметзянова А.Г.// Химическая промышленность. 2005. Т. 81, №3. С. 152-157.

126. Петровичева Е.А. Расчет осесимметричных турбулентных течений вязкой несжимаемой жидкости на основе программы CANAL 4S (5S). /Ильина И.М.// Ученые записки. Сборник научных трудов. Т.З, Альметьевск: АГНИ, 2005. С. 236-258.

127. Петровичева Е.А. Проблемы смешения однофазных компонентов в трубчатых турбулентных аппаратах. /Макарова Т.П., Ильина И.М.// Материалы научной сессии по итогам 2004 г., Альметьевск: АГНИ, 2005.

128. Петровичева Е.А. Методы численного решения задач для турбулентных течений смешивающихся жидкостей в МТТА. /Корепанова Л.Ф., Ильина И.М.// Ученые записки. Сборник трудов, Альметьевск: АГНИ, 2005, С. 119-132.

129. Петровичева Е.А, Численное моделирование турбулентных течений с помощью пакета. PHOENICS. /Макарова Т.П., Ильина И.М.// Материалы научной сессии по итогам 2005 г., Альметьевск: АГНИ, 2006.

130. Петровичева Е.А. Турбулизация потока в трубчатых аппаратах как способ повышения эффективности смешения компонентов. /Ильина И.М., Макарова Т.П.// Ученые записки. Сборник трудов. Том 4. Альметьевск: АГНИ, 2006, С. 225-232.

131. В.Д.Шашин исемендэге «ТАТНЕФТЬ» АЧЫК акционерлык ж;эмгыятеоткрытое акционерное ОБЩЕСТВО "ТАТНЕФТЬ"имени В.Д.Шашина

132. НЕФТЬ ьем ГАЗ ЧЫГАРУ ИДАРЭСЕ1. ЖЭЛИЛНЕФТЬ»

133. НЕФТЕГАЗОДОБЫВАЮЩЕЕ УПРАВЛЕНИЕ «ДЖАЛИЛЬНЕФТЬ»