Development and research of induction heating system in technological complexes of oil transportation

Cover Page

Cite item

Abstract

The problem of development and research of an induction system for heating viscous non-conductive liquids during transportation through pipelines is considered. A specific feature of the research object is the combination of the processes of heating a laminar fluid flow in a cylindrical inductor and mixing in a three-phase inductor - mixer. Research of interrelated electromagnetic, thermal and electromechanical processes in heating and stirring systems has been carried out. The results of a numerical calculation of electromagnetic and thermal fields and the electromagnetic moment of rotation of the rotor in a heated liquid are presented. A relationship has been established between the thermal power released in the elements of the system and the power used to rotate the rotor. The electromechanical (control) characteristics of the induction mixer are given. It is shown that with an appropriate choice of design and operating parameters of the induction system, it is possible to find the required ratio between the thermal power released in the elements of the system and the power used to stir the liquid. The results of the research are intended to solve the problems of designing the design of induction systems for heating viscous liquids and effective control of the heating process.

Full Text

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

Эффективность трубопроводного транспорта зависит от ряда факторов. К ним относится вязкость перекачиваемого продукта, от которой зависит производительность трубопроводной системы и энергозатраты на перекачку [1–4]. Одним из эффективных способов снижения энергозатрат на компримирование и увеличения пропускной способности трубопровода является снижение величины вязкости жидкости путем периодического подогрева с помощью путевых подогревателей, устанавливаемых с определенным интервалом вдоль магистрали. Подогрев осуществляется различными способами. Наиболее надежным и экологически чистым способом подогрева высоковязких и высокозастывающих нефтей в настоящее время является электроподогрев [6, 7].

В работах [5–9] рассматриваются вопросы расчета, проектирования и практического применения систем электрообогрева трубопроводов, резервуаров и технологического оборудования в нефтегазовой промышленности. Приведены описание и характеристики установок электрообогрева, работающих в различных отраслях промышленности. Однако вопросы моделирования и расчета охватывают ограниченный класс нагревателей с источниками энергии в виде греющих кабелей или индукционных систем с однослойной нагрузкой в виде металлической стенки резервуара или трубопровода.

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

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

Как следует из результатов исследований, такие свойства транспортируемой жидкости, как высокая вязкость, исключающая возможность создания турбулентности при технологически возможных скоростях потока, и низкая теплопроводность обусловливают необходимость применения многосекционных индукционных нагревателей, длина которых составляет 12÷16 м при требуемой температуре на выходе 50÷60 град. [11, 13]. Задача усложняется наличием дополнительного ограничения на максимальную температуру пограничных, наиболее нагретых слоев жидкости.

Решение задачи снижения массогабаритных показателей и повышения эффективности индукционных нагревательных устройств возможно за счет интенсификации процесса теплообмена между источниками тепла и нагреваемой жидкостью. Интенсифицировать процесс теплообмена и таким образом повысить эффективность нагревательной системы можно либо путем создания турбулентности потока за счет высоких скоростей, либо путем перемешивания жидкости в нагреваемом потоке [14–17]. При практически реализуемых скоростях потока нефти в трубопроводных системах имеет место ламинарное течение жидкости [17]. Средняя скорость потока составляет 2,8–3,3 м/с, причем вследствие высокой вязкости скорость по сечению потока распределена существенно неравномерно. Высокая скорость в средней части потока быстро падает к стенкам трубы. Пристенные слои потока прилипают к стенкам, а значительная величина вязкости способствует высоким касательным напряжениям между слоями и их торможению. При таких условиях создать турбулентность не представляется возможным.

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

 

Рис. 1. Схема нагревательного модуля: 1 – трубопровод; 2 – нагреваемая жидкость; 3 – статор трехфазного индуктора; 4 – ротор; 5 – однофазные индукторы; 6, 7, 8 – регулируемые источники питания; 9 – система управления нагревом

 

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

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

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

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

 

Рис. 2. Трехфазный индуктор с магнитопроводом: 1 – магнитопровод; 2 – обмотка статора; 3 – ротор; 4 – обмотка ротора

 

Для перемешивания слоев используется трехлопастная полая крыльчатка диагональной конструкции [19–21]. Центральная часть потока жидкости в полой крыльчатке лопастного колеса отклоняется к стенке трубы. Таким образом, объем жидкости, вошедший в лопастное колесо, попадает в пристенную область, вытесняя разогретый слой. Температура пристенного слоя снижается. Таким образом обеспечивается более интенсивная теплоотдача в жидкость. При этом вся энергия, подведенная к индуктору смесителя, также в конечном итоге преобразуется в тепло, обеспечивая дополнительный нагрев жидкости.

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

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

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

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

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

μaHr,x,t,THr,x,tt1rrr1σTHr,x,trxr1σTHr,x,tr=0;

H0,x,τr=0;

cTγTr,x,θ,ττ=1rrrλTTr,x,θ,ττ+1r2θλTTr,x,θ,τθ+xλTTr,x,θ,τθ+1cTγWr,x,θ,τ.

Здесь HR,L,t – вектор напряженности магнитного поля; μa – магнитная проницаемость; σT – электропроводимость; Tr,x,θ,τ – температурное распределение в системе сопряженных тел; τ – время; λT – удельная теплопроводность; Wr,x,θ,τ – удельная мощность тепловыделения; cT – удельная теплоемкость; γ – плотность, θ – угловая координата.

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

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

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

x1μyAx+y1μxAyjωgA=jext; AL=0.

Здесь: A – векторный магнитный потенциал; μx, μy – относительная магнитная проницаемость материала по осям x, y; j – мнимая единица; ω – круговая частота тока; g – удельная электрическая проводимость; L – граница расчетной области (линия удаленной поверхности); jext – плотность стороннего тока. На базе приведенных систем уравнений строится конечно-элементная формулировка для плоской двумерной области [23].

Граница раздела магнитных сред описывается системой соотношений:

B1n=B2n;

μ1H1n=μ2H2n;

H1τH2τ=dIdl.

Последнее выражение учитывает скачкообразное изменение вектора напряженности H¯ на границе раздела сред. Кроме условий сопряжения для получения однозначного решения уравнений Максвелла в форме напряженности электрического поля E¯ и напряженности магнитного поля H¯ в области QR3 с границей S необходимо задать:

 – уравнения поверхностей, отделяющих друг от друга среды i и j, fijx,y,z=0;

 – начальные величины E0x,y,z, H0x,y,z в момент времени t0 в произвольной точке исследуемого объема QR3 с границей S;

 – касательные составляющие вектора E¯ или H¯ в произвольной точке поверхности в произвольном временном интервале от t0 до t, или распределения полей E¯ и H¯ вне исследуемого объема Q;

 – функциональные зависимости магнитной проницаемости μ и удельной проводимости γ от пространственных координат или от напряженности поля.

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

B¯=rotA¯;

E¯=A¯tdivV.

Чтобы функция A¯ была определена, нужно определить значение ее дивергенции. Для этого добавляется условие, которое называется калибровкой Кулона:

divA¯=0.

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

rot1μrotA¯+σA¯t=J¯;

rotA¯=B¯;

divA¯=0.

Используя соотношение

rotrotA¯=graddivA¯2A¯,

при μ=const получим уравнение

2A¯jωσA¯J¯=0.

Уравнение Пуассона дополняется граничными условиями Дирихле и Неймана на различных участках границы:

A¯=0 на S1; A¯n=0 на S2.

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

Численные эксперименты выполнены для устройства с параметрами, приведенными в табл. 1.

 

Таблица 1. Исходные данные

Наименование

Величина

Диаметр трубы, мм

169

Материал трубы, мм

Нержавеющая сталь

Толщина стенки трубы, мм

4,5

Диаметр цилиндрического ротора внешний, мм

82

Толщина стенки цилиндра, мм

12

Диаметр цилиндрического ротора внутренний, мм

58

Материал цилиндрического ротора

Ферромагнитная сталь

Размер паза (высота×ширина), мм

4×7

Размер зуба (высота×ширина), мм

4×7

Обмотка ротора

Шинка алюминиевая, 4×6

Количество пазов, шт.

12

Диаметр вала, мм

20

Диаметр статора внутренний, мм

86

Диаметр статора внешний, мм

158

Материал статора

Электротехническая сталь

Длина статора, мм

400

Размер паза (высота×ширина), мм

16×10

Размер зуба (высота×ширина), мм

16×10

Частота тока индуктора, Гц

50

Кол-во пазов, шт.

12

 

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

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

 

Рис. 3. Распределение плотности тока в элементах конструкции: 1 – в катушке индуктора; 2 – в короткозамкнутой обмотке ротора; 3 – в стальном теле ротора

 

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

 

Таблица 2. Характеристики трехфазного индуктора

Частота источника питания

50

Ток, А

6300

Напряжение, В

4,9

Момент, Н/м, кз

106,9

Частота вращения, об/мин, идеального холостого хода

314

Полная мощность, ВА

39 633

Мощность реактивная, ВАр

31 300

Полная активная мощность, Вт

26 372

Мощность на нагрев в роторе, Вт

15652

Мощность на вращение, Вт

10 744

Доля мощности на нагрев

0,44

cosφ

0,65

 

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

 

Рис. 4. Механические характеристики трехфазного идуктора

 

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

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

Еще одна проблема заключается в согласовании электромеханических характеристик трехфазного индуктора с характеристиками смесителя. Сложность состоит в отсутствии методики расчета момента вращения осевой крыльчатки в вязких средах. Необходимый для перемешивания нефти момент вращения рассчитывается по методике, предложенной в работе [23, 24]. Методика расчета основана на использовании коэффициентов приведения, устанавливающих связь между моментом вращения в воде и моментом вращения в вязкой среде. Алгоритм расчета представлен на рис. 5. На первом этапе выполняется расчет по методике, в которой в качестве рабочей жидкости рассматривается вода. Затем с использованием коэффициентов приведения производится пересчет мощности и момента для вязкой жидкости с учетом ее реальных характеристик.

 

Рис. 5. Алгоритм расчета мощности и момента вращения

 

Исходные данные, необходимые для расчета, приведены в табл. 3.

По исходным данным определяется напор, создаваемый насосом

Hос=υк22gkн2

и скорость

Vc=kν2gHос,

где Kн2=0,0244ns0,66, ns – коэффициент быстроходности.

 

Таблица 3. Исходные данные для расчета момента вращения

Наименование

Величина

Размерность

Обозначение

Скорость потока до установки

0,05

м/с

vн

Скорость потока после установки

0,1

м/с

vk

Давление потока до установки

10

кгс/см2

pн

Давление потока после установки

12

кгс/см2

pk

Коэффициент быстроходности

30

 

ns

Угловая скорость

37

рад/с

ω

Плотность среды (вода)

1000

кг/м3

ρв

Плотность среды (нефть)

920

кг/м3

ρн

Вязкость (нефть)

4,5

см/сек2

υ

 

Полезная мощность:

Pпол=QидpвgHос102,

где Qид – идеальная подача осевого насоса:

Qид=π4D2d2VZ,

где D – диаметр колеса; d – рекомендуемый диаметр ступицы; VZ – осевая скорость; d=0,4.....0,5D; VZ=kυ2gHос; kυ=0,055ns0,66.

Мощность установки с учетом потерь:

P=1,2Pпол.

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

 

Таблица 4. Данные для перерасчета на нефть

Наименование

Величина

Размерность

Обозначение

Толщина лопатки

0,03

м

b

Ширина лопатки

0,25

м

q

Число рядов лопаток

2

шт.

z

 

Определяется коэффициент стеснения:

K2=πDqzπD=t2qt2,

где t2 – шаг лопаток.

Эквивалентный диаметр рабочего колеса:

D2экв=4dbK2.

При максимальной идеальной подаче насоса находится число Рейнольдса для вязких жидкостей:

Re=QрD2эквυ.

Дальнейший расчет выполняется с помощью опытных поправочных коэффициентов Kн, Kq, Kη, которые при найденных числах Рейнольдса находятся из графика [23, 24] (рис. 6).

 

Рис. 6. Зависимость коэффициентов от числа Рейнольдса

 

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

Qв.ж.=KqQ; Hв.ж.=KHH; ηв.ж.=Kηη.

Полезная мощность в среде с вязкой жидкостью:

P=Qв.ж.Hв.ж.gρв.ж.75ηв.ж..

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

По результатам расчета для исследуемого устройства с приведенными выше параметрами момент трогания составляет 51,4 Нм. По результатам расчета электромагнитной задачи момент трогания (пусковой момент), развиваемый устройством, составляет 106,9 Нм, что обеспечивает работу смесителя в пусковом и установившемся режимах.

Вывод

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

×

About the authors

Aleksandr I. Danilushkin

Samara State Technical University

Email: aidanilushkin@mail.ru

Dr. Sci. (Techn.), Professor

Russian Federation, 244, Molodogvardeyskaya st., Samara, 443100

Vasiliy A. Danilushkin

Samara State Technical University

Author for correspondence.
Email: aidanilushkin@mail.ru

PhD (Techn.), Associate Professor

Russian Federation, 244, Molodogvardeyskaya st., Samara, 443100

References

  1. Gubin V.E., Gubin V.V. Pipeline transport of oil and oil products. M.: Nedra, 1982. 296 pp.
  2. Pipeline transport of oil: SB. nauch. trudov. Ufa: VNIIS-Ptneft, 1987. 136 pp.
  3. Tugunov P.I. Non-Stationary modes of oil and oil products pumping. M.: Nedra, 1984. 224 pp.
  4. Nadirov N.To., Chugunov P.I. Pipeline transportation of viscous NEF-Tay. Alma-ATA: Nauka, 1985. 146 pp.
  5. Konesev S.G., Khlyupin P.A., Sadikov M.R. Analysis of the effectiveness of heating systems when pumping viscous oils electric Drive, electrical technologies and electrical equipment of enterprises // SB. nauch. Tr. III All ROS. nauch. – tehn. Conf. (from the international school). Ufa: USNTU, 2011. P. 211–218.
  6. Konesev S.G., Khlyupin P.A. Environmental heating systems for oil transport and storage facilities // Life safety, 2012. № 7. P. 35–42.
  7. Trofimenko K.V., Shishkin N.D. Design development and evaluation of parameters of an induction heater for oil and petroleum products // AGTU Bulletin, 2012. № 1(53). P. 78–83.
  8. Shishkin N.D. Application of induction heating of oil during its transportation from fields in the Northern Caspian // Vestnik AGTU. 2011. № 3. P. 52–56.
  9. Strupinsky M.L., Khrenkov N.N., Kuvaldin A.B. Design and operation of electric heating systems in the oil and gas industry. M.: Infra-Engineering, 2015. 272 рp.
  10. Danilushkin A.I., Bazarov A.A., Zinnatulin D.A. Investigation of electromagnetic and thermal fields in the installation of technological heating of oil // Vestnik SamGTU. Technical Sciences, 2004. Vol. 24. P. 171–173.
  11. Danilushkin V.A. Optimization of design and operation modes of induction heaters of high-viscosity oils during transportation through pipelines // Vestnik SamGTU. Technical Sciences, 2004. № 20. P. 176–179.
  12. Bazarov A.A. System of induction heating of a moving liquid // Vestnik SamGTU. Technical Sciences. 2005. № 37. P. 12–17.
  13. Bazarov A.A., Danilushkin A.I., Danilushkin V.A., Vasiliev I.V. Modeling of electromagnetic processes in multilayer three-phase induction cylindrical system // Vestnik SamGTU. Technical Sciences. 2017. № 3(55). P. 50–60.
  14. Ciofalo M., Brucato A., Grisafi F., Torraca N. Turbulent flow in closed and free surface unbaffled tanks stirred by radial impellers // Chemical Engineering Science, 1996. Vol. 51. P. 3557–3573.
  15. Haque J.N., Mahmud T., Roberts K.J., Rhodes D. Modeling turbulent flows with free surface in unbaffled agitated vessels // Industrial and Engineering Chemistry Research, 2006. Vol. 45. P. 2881–2891.
  16. Gazizullin N.A. Mixing of a viscous liquid with a free surface in a device with a Rushton turbine agitator // Vestnik SamGTU. Technical Sciences, 2014. № 3(42). P. 146–155.
  17. Vasiltsov E.A., Ushakov V.G. Apparatuses for mixing liquid media. L.: Mashinostroenie, 1979. 272 pp.
  18. Danilushkin A.I., Danilushkin V.A., Maksimova M.A., Surkov D.V. Development and research of a three-phase induction device for heating and mixing a liquid // Vestnik SamGTU. Technical Sciences, 2019. № 3(63). P. 120–132.
  19. Cherkassky V.M. Pumps, fans, compressors. M.: Energoatomizdat, 1984. 424 pp.
  20. Bashurov B.P. Ship pumps and fans. I/O Marterella, 1983. 32 pp.
  21. Belousov A.N., Musatkin N.F., Radko V.M. Theory and calculation of aircraft shovel machines. Samara: Samara house press, 2003. 336 pp.
  22. ELCUT. Modeling of two-dimensional fields by finite element method. User manual. Version 5.9. SPb: TOR production cooperative, 2009.
  23. Karpushkin S.V., Krasnyansky M.N., Borisenko A.B. Calculations and selection of mechanical mixing devices for vertical capacitive apparatuses: textbook. Tambov: TSTU, 2009. 168 pp.
  24. Yuriev A.S., Pirogov S.Yu., Nizovtsev V.M., Grachev I.G., Presnov A.I., Savishchenko N.P., Sokolova A.N. Reference book on the calculations of hydraulic and ventilation systems. SPb: ANO NPO “Peace and Family”, 2001. 1154 pp.

Supplementary files

Supplementary Files
Action
1. Fig. 1

Download (161KB)
2. Fig. 2

Download (211KB)
3. Fig. 3

Download (99KB)
4. Fig. 4

Download (163KB)
5. Fig. 5

Download (154KB)
6. Fig. 6

Download (129KB)

Copyright (c) 2022 Samara State Technical University

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies