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

• Волков, Валерий Иванович
• доктор технических наук
• 1998

Оглавление диссертации доктор технических наук Волков, Валерий Иванович

Оглавление

Введение

Некоторые исследования по гидродинамики и процессам

переноса

1 Аппаратура для исследования двухфазных потоков

1.1 Установки для гидродинамических исследований

1.2 Измерительные схемы

1.2.1 Контактные методы измерения скорости и концентрации

1.2.2 Оптические методы исследования

1.3 Физико-химический анализ иммерсионных жидкостей

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

2.1 Распределение давления в упаковках

2.2 Лазерная анемометрия структуры течения

2.2.1 Измерения в кубической упаковке

2.2.2 Поле скоростей в хаотичной упаковке

2.3 Неустойчивость течения жидкости в засыпках

2.3.1 Нестационарность пристеночного давления

2.3.2 Нестационарность при измерениях скорости жидкости

2.3.3 Свободноконвективная неустойчивость в засыпках

2.4 Гидродинамическая обстановка в стационарных упаковках

3 Математические модели фильтрации и переноса

3.1 Некоторые модели течения

3.2 Расчет модели переноса

3.3 . Сравнение теоретических и экспериментальных иссле-

дований

4 Исследование разрушающихся упаковок

4.1 Исследование распределения давления и скорости

4.1.1 Измерение давления

4.1.2 Измерения скорости за упаковкой

4.2 Тепловые исследования разрушаемых упаковок

5 Газогенераторы

5.1 Газогенераторы с магнитной подачей

5.2 Газогенерирование больших объемов

6 Рациональное использование природного топлива

6.1 Сжигание угольной пыли

6.2 Исследование топок с кипящим слоем

6.3 Исследование процессов шлакования и методы воздействия на пористую среду

Заключение

Примечание

Литература

Введение диссертации (часть автореферата) на тему «Исследование аппаратов двухфазных потоков в ресурсосберегающих технологиях»

Введение

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

1. Методика и техника измерения параметров потока.

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

3. Решение проблем масштабного перехода от модели к натуре.

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

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

Автор защищает:

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

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

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

4.Создание нового направления в газогенераторах на гибридном то-

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

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

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

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

8.Получение новых экологически - безопасных иммерсионных жидкостей. Их синтезирование и физико-химический анализ.

Некоторые исследования по гидродинамики и процессам

переноса

Отличительной чертой движения жидкости через зернистые слои является неравномерность распределения скоростей по сечению слоя. Подобная неоднородность потока приводит не только к снижению эффективности работы аппарата с зернистым слоем, но и часто к локальному перегреву и запеканию зерен слоя в горячем газовом потоке или к замораживанию отдельных участков рабочего элемента в теплообменниках, к усилению капельного уноса в фильтрующих аппаратах, что иногда вызывает полное нарушение рабочего цикла и выход аппарата из строя. Исследованию гидродинамики и переносу в зернистых слоях посвящено большое количество работ [1] - [34]. В работах Стру-минского с сотрудниками показано наличие крупномасштабных неод-нородностей в распределении скорости по сечению опытного участка и предложены некоторые методы устранения этих неоднородностей. С помощью визуализации Струминским и Павлихиной была зафиксирована неоднородность профиля скорости в зернистых слоях [4]. Следует отметить, что профиль скорости за упаковкой может отличаться от профиля скорости внутри ее. Так в работах Кириллова, Кузьмина и др. [6], Леру и Фромента [7] показано, что профиль скорости за пористой средой существенно зависит от ориентации и упаковки двух-трех последних слоев засыпки или от прогиба сетки, ограничивающей

засыпку. Работ, посвященных измерению скорости внутри пористой средкц сравнительно мало, что объясняется сложностью проведения подобных экспериментов. Так известны две работы (Миклей, Смит, Корчак [9], Ван дер Мерве и Гаувин [8]) по измерению термоанемометром профиля скорости и интенсивности турбулентности потоков в правильных упаковках. Миклей и др. измеряли профиль локальной скорости, интенсивность турбулентности и спектр пульсаций скорости в ромбической упаковке для чисел Re, равных 4780 и 7110 при D/d=8, где D - ширина квадратного канала, d - диаметр шарика. В этой работе получено, что профили интенсивности турбулентности и локальной скорости становятся более пологими при увеличении числа Re, так отношения максимальных скоростей в ячейках к расходной скорости при числе Re=4780 и 7010 составляли 6,1 и 4,5 соответственно. Кроме того, скорость вблизи стенки упаковки оказалась на 10 % выше скорости в центре, хотя были приняты меры по устранению неравномерности пористости у стенки, в которой крепились половинки сфер. Профиль средней скорости, полученный ими, имел два пика, в соответствии с количеством пропускных каналов для элементарной ячейки ромбической упаковки. Уменьшенное значение скорости между пиками соответствовало теневой области за шариком, в которой возможно образование возвратных течений, необнаруженных авторами. С увеличением числа Re максимальные значения скоростей в пиках сближаются, т.е. профиль становится более равномерным в отдельных ячейках упаковки. Результаты Ван дер Мерве и Гаувина по измерению интенсивности турбулентности на правильной кубической упаковке при числах Re = 2500, 5000, 10000, 27000 и D/d = 4 качественно совпадают с результатами Миклея и др., но в спектре пульсаций скорости, полученным Ван дер Мерве и Гаувином, имелись пики, тогда как в спектре, полученном Миклеем и др. пики отсутствовали. Кроме того, Ван дер Мерве и Гаувин отмечают увеличение крутизны профи-

ля продольной скорости с увеличением числа 11е. Эти результаты расходятся с выводами работ Евсеева, Накорякова, Романова [10, 11], которые измеряли спектр пульсаций и профиль скорости внутри кубической и октаэдрической упаковок с помощью лазерного доплеров-ского анемометра. Известны еще две работы по измерению профиля скорости внутри хаотической упаковки электродиффузионным методом (Акехато и Сато [12], Кириллов, Кузьмин и др. [6]). Результаты этих двух работ в основном согласуются между собой, но отличаются от результатов Миклея и др. Принципы электродиффузионной диагностики подробно описаны в работе [3]. Акехато и Сато использовали электродиффузионные датчики в виде 5 концентрических колец, выполненных из платиновой проволоки диаметром 0,1 мм. Эти кольца располагались в одном сечении стеклянной трубы, заполненной шариками, диаметр которых варьировался от 1 до 7 мм. Кириллов, Кузьмин и др. использовали рабочий участок, диаметр которого был равен диаметру стеклянной трубы в экспериментах Акехато, Сато (50 мм). В рабочий участок загружались трех миллиметровые стеклянные шарики. В качестве электродиффузионных датчиков Кириллов, Кузьмин и др. так же как и Акехато, Сато использовали пять концентрических колец из никелевой проволоки диаметром 0,28 мм, расположенных в одном сечении рабочего участка на радиусах б, 11, 15, 20 и 24 мм. Кроме кольцевых датчиков Кириллов, Кузьмин и др. использовали никелевые шарики диаметром 3.8 мм и проволочные щупы диаметром 0.28 мм и длиной 0,5 - 1,5 мм. Профиль скорости, замеренный Акехато и Сато для чисел Яе= 0.4 - 140 и Б/с1= 16 - 54 и Кирилловым и Кузьминым для Яе= 20 - 600 и Б/с1=16, найден плоским. Но для отношения Б/с1=7.2 и 9.4 в экспериментах Акехато и Сато профиль скорости оказался выпукло-вогнутым у стенки ¿¿а расстоянии, близком к диаметру зерна. Как отмечают авторы область с минимальной средней скоростью у стенки соответствует области с минимальной по-

ристостью, находящейся примерно на расстоянии радиуса зерна от стенки. При экспериментах с разными Б/с1 Акехато и Сато не меняли положение датчиков в трубе, поэтому, когда размер зерна становился меньше, чем расстояние до второго датчика от стенки, профиль скорости становился плоским, т.е. датчики не чувствовали изменение потока на первом зерне от стенки. В этом случае их результаты естественным образом совпадали с результатами работы Кириллова, Кузьмина и др., при использовании ими датчика величиной больше, чем диаметр зерна. В измерениях Кириллова, Кузьмина и др. проволочным щупом, длина которого меньше, чем диаметр зерна, скорость жидкости в центре упаковки оказалась несколько меньше скорости жидкости у стенки. Из этих двух работ следует, что если неравномерность профиля скорости у стенки и существует, то она невелика и, при усреднении по размеру зерна получается среднее значение скорости. Аналогичный результат, говорящий о том, что профиль в хаотичной упаковке плоский, получил Кубота [13]. Кубота измерял поле скоростей по теплообмену от единичной сферы диаметром 9,6 мм, в точности равном диаметру зерна засыпки. Он измерил поле скоростей на хаотичной упаковке в диапазоне чисел Ке= 7 - 173 в прямоугольном канале, поперечное сечение которого, отнесенное к квадрату диаметра зерна, составляло 9 х 26. Кубота отметил, что хотя локальные скорости жидкости в центре упаковки и у стенки равны, тем не менее, расход жидкости у стенки существенно выше, за счет большого количества пустот у стенки по сравнению с центром хаотичной упаковки. Рассмотрим несколько работ по бесконтактным методам определения скорости в зернистом слое. Так во времяпролетном методе Колесанова скорость газа определялась с помощью радиоактивной метки (радона), вводимой взрывным образом в некоторых участках доменной печи [15]. Хотя в работе не приводятся полные профили скорости по сечению зернистого слоя, тем не менее, представляет ин-

терес исследование автором извилистости и проходного сечения для газа в зернистом слое из кускового материала. Колесанов отметил, что при низких числах Бе газ как бы "блуждает" по слою, "разыскивая" каналы с наибольшим проходным сечением. Чем меньше число 11е, тем по меньшему количеству каналов движется газ в зернистом слое. Это приводит к увеличению извилистости и уменьшению проходного сечения для газа в зернистом слое при понижении числа Яе. С увеличением расхода газа извилистость падает, а проходное сечение растет до определенного предела. Этот предельный расход газа может быть найден из условия перехода во взвешенный слой и соответствует прохождению газа практически по всем каналам в слое материала.

В работах Аэрова, Умника [25], Табунщикова [16] использовалась пропорциональность скоростей газа и фронта сорбции какой-либо примеси из газа. Визуальные наблюдения за движением фронта сорбции позволяли определять скорость газа для не слишком высоких чисел Яе, при которых еще фронт сорбции имеет отчетливые границы. Этим методом до чисел Яе < 40 Аэров, Умник получил, что скорость газа у стенки на 30-70% выше скорости газа в средних слоях. Табунщи-кову удалось применить этот метод для более высоких чисел Яе. Он установил, что 30% превышение скорости у стенки практически не изменяется при изменении числа Яе от 3 до 160. В работах Абаева, Попова, Смирновой и др. [17, 18] скорость газа определялась из однозначной связи между временем регенерации катализатора и количеством прошедшего через зернистый слой воздуха. Авторы получили рост крупномасштабных неоднородностей в профиле скорости с понижением числа Яе. Так при числах Яе= 150 и 3000 скорость газа у стенки превышала скорость газа в центре опытного участка в 6 и 2.5 раза соответственно. Известна также работа Каринса и Праузнит-ца [14], которые измеряли профиль скорости в хаотической упаковке для чисел Яе= 9 - 1700 и Б/с1=15.6. В этой работе использовалась

времяпролетная методика. В начале упаковки впрыскивалась проводящая метка. Чувствительные электроды были вмонтированы в скрепленные между собой шарики, имеющие диаметр, близкий к диаметру зерна засыпки (3,2 мм). Конструкция чувствительного элемента не позволила Каринсу и Праузнитцу подойти к стенке ближе, чем на два диаметра зерна. В этих измерениях авторы получили плоский профиль скорости с незначительным увеличением скорости в центре упаковки, но значение скорости отличалось от расходной. Это отличие автора отнесли к влиянию пристенной части упаковки, в которой они не могли измерить скорость. Следовательно, в пристенной части упаковки скорость у них оказалась выше, чем в центре. Таким образом, во всех бесконтактных методах определения скорости в зернистом слое получено повышенное значение скорости у стенки по сравнению с центром опытного участка. Но для каждого метода определения скорости величина превышения скорости у стенки и зависимость этого превышения от числа Ле оказалась различной. Тем не менее, большинство исследователей утверждало, что с увеличением числа Яе профиль средней скорости становился более плоским по сечению упаковки. Так Карине, Праузнитц получили для чисел Ие < 1700 почти 20% превышение скорости у стенки, независящее от числа 11е. Аэров, Умник Табунщиков получили для чисел 11е < 160 30% превышение скорости у стенки, практически независящее от числа Ее, но для чисел Яе > 600 Аэров, Умник предполагали уменьшение неоднородности профиля скорости до 10 - 20%. Абаев, Попов, Смирнова и др. также получили, что с ростом числа Т1е профиль скорости становится более однородным. Хотя и при числах Ие—3000 скорость газа вблизи стенки у них оказалась на 150% выше скорости в центре опытного участка. Следует отметить, что хотя результаты Каринса, Праузнитца противоречат результатам Куботы, Акехато, Сато, Кириллова, Кузьмина и др., в работах которых для хаотичных упако-

вок получены плоские или почти плоские профили скорости, тем не менее, следует учесть, что времяпролетная методика имеет немаловажное отличие от методик определения скорости, основанных на тепломассообмене. Так времяпролетная методика Каринса, Праузнитца требует измерения проводимости между двумя электродами, расположенными на расстоянии 1 мм. Проводимость измеряется после прекращения подачи проводящей метки и зависит не только от локальной скорости жидкости, но и от скорости, усредненной по сечению потока между электродами. В диагностиках, основанных на измерениях диффузии тепла или массы через пограничный слой теплового или электродиффузионного датчиков, фактически определяется скорость его обтекания. Поэтому следует ожидать, что времяпролетные методики будут более чувствительны к изменению пористости или расходной скорости, чем методики, основанные на тепломассообмене. Далее принимая во внимание замечание Куботы об увеличенном значении расходной скорости у стенки упаковки по сравнению с центром, при равенстве локальных скоростей в центре и у стенки, можно частично объяснить противоречие отмеченных работ. Из приведенных работ можно заметить, что локальная структура потока в зернистых слоях недостаточна ясна и результаты некоторых работ противоречивы. Это связано с известными недостатками контактной методики такими, как искажение пористой среды и потока, зависимость показания датчиков от расстояния до твердой поверхности шариков или стенки канала, сложность в определении вектора локальной скорости потока. С другой стороны, возможно и влияние вида упаковки, геометрии рабочего участка, режима течения жидкости. Так, например, Миклей и др. изучали поле скоростей на правильной упаковке при достаточно высоких числах Яе > 4500, т.е. в турбулентной области, где они нашли, что локальная скорость вблизи стенки на 10% выше скорости в центре. Кириллов, Кузьмин и др. изучали поле скоростей в хаотич-

ных упаковках при числах 11е < 600, а Акехато, Сато и Кубота при числах Ие < 200, т.е. ближе к ламинарному течению, и получили плоский или почти плоский профиль локальной скорости. Эти результаты не согласуются с результатами Каринса, Праузнитца, которые практически для всего диапазона чисел Яе= 9 - 1700 получили скорость у стенки на 20% выше скорости в центре упаковки, и противоречат результатам Ван дер Мерве, Гаувина, у которых с повышением числа Яе профиль скорости становится более крутым.

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

где Юи = к Я2 и\/481)1 - наличие пристенного эффекта аналогичного влиянию "тейлоровской диффузии"; /г-численный коэффициент; £>о = Ве/Т - коэффициент диффузии для зернистого слоя; е - пористость; Т- извилистость; ИК = к\щй - конвективная составляющая для дисперсии; д, - диаметр зерна; - релаксационная составляющая дисперсии, для больших чисел 11е 1)р ~ щ<1, для малых Ир ~ и^/В; -Оф ~ Щ<1 ~ флуктуационная составляющая дисперсии.

Объединяя все составляющие в одну формулу, в первом приближении можно записать:

где и некоторые константы.

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

Вьф - Д, + £>к + £>р + 4- £>ф,

(0.1)

Дэф/1)1 = к2 + кзИеБс,

(0.2)

много исследований [25] - [34]. В основном все они базируются на двух различных подходах, обусловленных разной гидродинамикой внутреннего и внешнего обтекания. Так капиллярная модель связана с течением внутри капилляров и дисперсионные эффекты в ней объясняются различными пропускными способностями и ориентацией капилляров [27]. Как отмечено в работе [2], капиллярная модель применима больше к сцементированным грунтам, чем к зернистым средам. Насыпные слои, по всей видимости, ближе к ячеистой модели пористой среды, в которую входят такие признаки внешнего обтекания, как наличие неравнодоступных объемов или застойных зон, особенно вблизи точек контакта элементов среды [28] - [33]. Как видно из соотношений (0.1, 0.2) в эффективные коэффициенты переноса тепла и массы в зернистых средах могут одновременно входить несколько физических механизмов переноса, что затрудняет как разработку теоретических моделей, так и сравнение теории с экспериментом. Еще более усложняется картина в нестационарных зернистых слоях, разрушаемых в процессе течения жидкости или газа, например, в охладителях или газогенераторах при наличие в потоке частиц твердой фазы.

В заключение обзора по гидродинамике и переносу в зернистых слоях приведем известные данные по обтеканию одиночного шарика и укладки из шариков при различных числах Ие [19], [24], [25]. Для хаотичной укладки из шариков при числе Ие < 20 режим течения жидкости ламинарный, с плавным обтеканием жидкостью отдельных элементов. Механизмы переноса тепла и массы в основном молекулярные. Коэффициент продольной и поперечной теплопроводности при Ые < 10 равны друг другу. При числе Ие 10 - 20 за шариками появляются вихри, струи между шариками медленно изменяют свое направление. При числе Ые больше 20, но меньше 50 вихри то зарождаются, то исчезают.

При обтекании одиночного шарика стационарный кольцевой вихрь

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

При числах Яе > 20 и примерно до числа Яе ~ 100 коэффициенты переноса тепла как вдоль, так и поперек потока в зернистом слое не зависят от числа Яе. В диапазоне чисел Яе 50-140 существует устойчивый ламинарно-вихревой режим, с четко очерченной транзитной струей, ограниченной вихревой структурой. При числе Яе = 115 - 140 происходит переход к турбулентному режиму течения.

За одиночным шариком вихревое кольцо увеличивается в размерах и начинает осциллировать при числе Яе=130. За кольцевым вихрем образуется ряд искаженных вихревых петель, несимметричных относительно оси.

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

Для одиночного шарика при числе Яе=500 осцилляции становятся периодическими, от кормовой области с определенной частотой, растущей с числом Яе, отрываются вихревые кольца и уходят вниз по потоку в виде вихревой дорожки Кармана. При числе Яе ~ 3 105 турбулизуется пограничный слой и уменьшается коэффициент сопротивления шарика.

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

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

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

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

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

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

Заключение диссертации по теме «Охрана окружающей среды и рациональное использование природных ресурсов», Волков, Валерий Иванович

Выводы.

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

Выявлены температурные интервалы наиболее активного шлакования для Кузнецкого, Березовского, Экибастузкого, Карагандинского углей.

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

1.15 Т/Т1

1.10

1.05

-Г 1 - * 2- • 1

- ^У \

1.00

0.0 0.5 1.0 1.5 2.0

Ь, кг/кг

Рис.6.1

Относительная темепература воспламенения в зависимости от весового соотношения между топливом и воздухом. 1 и 2 - кузнецкий и березовский угли, фракции 63-100 мкм, Т1 - 1000 К и пыль фракции 90 мкм меньше 15%, Т1-810 К.

1200 Т,К

1100

1000

900

800

-

-

1 - ж - 2- А 3- •

10 f, Гц

Рис.6.2

Темепература воспламенения в зависимости от частоты пульсаций расхода топлива. 1, 2 и 3 - кузнецкий, экибастузкий и березовский угли, фракции 63-100 мкм, избыток воздуха 1.2, фракции 100-160 мкмк, избыток воздуха 2.4 и пыль фракция 90 мкм меньше 15%, с недостатком воздуха 0.5.

Рис. 6.3

Профили скорости и концентрации в запыленном потоке.

1 - вертикальная скорость, К 1=2, К2=5 , КЗ=8 м/с;

2 - расходная концентрация, К1=0.2, К2= 0.5, К3=0.8 кг/м (-3 - счетная концентрация, N/N111, где Ы-разность между количеством частиц, двигающихся вверх и вниз, №п - ьаксимальное значение концентрации, К1=0.2, К2=0.5, К3=0.8.

5,

ММ 1,2

0,8

0,4

0 0,2 0,4 0,6 х, м

Рис. 6.4

Распределение среднего размера частиц пыли на разных высотах форсированного кипящего слоя. 1, 2, 3 - расстояния Ъ от диффузора равны 0.6, 1,

1.8 м.

230

Рис. 6. 5

Распределение скорости запыленного потока на разных высотах форсированного кипящего слоя. 1, 2, 3 - расходные массовые концентрации частиц равны О,

0.5, 2 кг/кг.

Х5м

Рис.6.6

Профили скорости воздушного потока, а и б - без пыли и с пылью. 1 - 03=0.345 м /с, 0^0.034 м/с, 0.099 м1/с; 2- 0=0.45 м'/с, (2г=0.13 м'/с, 02=0;

3 и 4 - 02=0, О0"=0.45 м'/си 0=0.58 м/с; 5 - 0 (Г 0.45 м/с, 0=0, 0=0.13 м'/с.

Рис.(6.7

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

Рис. 6. В

Шлакование Карагандинского угля при температуре 1270 К. 1, 2 и 3, 4 - без дутья и с дутьем воздуха. 1, 3 и 2, 4 - расстояния от сопла 0.5 и 15 см.

Рис. 6.9

Шлакование углей без дутья воздуха. 1, 2 - Кузнецкий уголь, Т=1070 К, расстояния от сопла 2.5 и 15 см. 3,4, и 5 -Карагандинский, Экибастузкий и Березовские угли при температурах 1070, 1020 и 820 К при расстояниях от сопла 1.5, 10 и 10 см.

Заключение

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

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

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

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

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

Впервые получены бесконтактными измерениями в хаотичных засыпках профили скорости внутри различных упаковок из шариков и распределение пристенного перепада давления по длине засыпок из шариков и таблеток. Полученные результаты были использованы на моделях перспективных экологически безопасных и малоэнергоемких ядерных реакторов на шаровых твэлах в виде правильной кубической упаковки и включили в себя новые данные по полям скоростей, полученные с помощью лазерной анемометрии в диапазоне чисел Яе < 1000 и распределение пристеночного давления в хаотичных упаковках до чисел 11е=20000. Экспериментально получено, что при числе Рейнольдса меньше 100 влияние входного эффекта на распределение пристеночного давления в хаотичной упаковке простиралось на расстояние, примерно равное 30 диаметрам элемента упаковки. При этом обнаружилось влияние низкочастотных акустических пульсаций (частотой меньше 100 Гц) на распределение пристеночного давления до не слишком высоких чисел 11е (меньше 1000), при которых амплитуды акустических пульсаций давления меньше амплитуды турбулентных пульсаций. Дискретная составляющая спектра внутри правильной кубической упаковки, соответствовала числу Струхаля 3.5, и была существенно выше, чем для обтекания одиночного шарика, что, по всей видимости, обусловлено меньшим размером вихрей из-за загромождения потока.

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

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

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

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

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

5.Выявлены оптимальные частоты пульсации первичного воздуха с пылью Кузнецкого, Карагандинского, Березовского, Бикинского и Итатского углей, связанные с температурой воспламенения топливной смеси и временем ее прогрева. Получено уменьшение температуры воспламенения угольной пыли при ее пульсирующем расходе на 50-100 градусов. Разработаны конкретные устройства пульсаций первичного и вторичного воздуха, защищенные 7 патентами России, основанные на управлении, как средиерасходными характеристиками потока и температуры, так и на изменении пульсаций и профиля скорости первичного и вторичного воздушных потоков.

6.Разработаны методики измерения скорости и счетной концентрации запыленного потока, с помощью которых удалось зафиксировать в форсированном кипящем слое сильное расслоение частиц по фракционному составу: крупные частицы мигрировали к стенкам, а мелкие концентрировались по оси потока и подвергались значительному выносу из модели топки. Это приводило к жестким требованиям к вну-тритопочным уловителям циркулирующих частиц. Установлено, что увеличение массовой расходной концентрации циркулирующих частиц до 0,5 кг/кг, также как и боковой вдув, устраняли неоднородность профиля скорости воздушного потока, обусловленную эффектом Ко-анда. Дополнительный вдув через боковые решетки (до 20% от полного расхода) существенно увеличивал придонную циркуляцию частиц и практически не сказывался на равномерности профиля скорости в надслоевом пространстве модели. Полученные данные необходимы при реконструкции слоевых котлов (ДКВР) производительностью до

10 т/ч.

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

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

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

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

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

Список условных обозначений

Б - ширина квадратного канала или его диаметр; коэффициент диффузии.

Я - радиус канала; радиус шарика; параметр. с! - диаметр шарика или таблетки, плотность жидкости.

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

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

Во - коэффициент диффузии для зернистого слоя. £ - пористость.

Т - извилистость; параметр; температура газа. ИК - конвективная составляющая для дисперсии. £)р - релаксационная составляющая дисперсии. Дф - флуктуационная составляющая дисперсии. 5с - число Шмидта.

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

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

С - емкость измерительного прибора и измерительной линии; коэффициент проницаемости слоя.

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

М - масса частицы, ударившейся о поверхность датчика; усредненная масса частицы; параметр. V - скорость частицы.

- средняя плотность частиц; расход воздуха через установку с кипящим слоем.

N - число частиц.

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

Яе = щс1/р - число Рейнольдса. Аг - критерий Архимеда. д , 1>„ р8 - ускорение свободного падения, кинематическая вязкость, плотность частиц. р, с1 - плотности жидкости; частицы. и, и* текущее и максимальное значение скоростей. щ -пристеночное значение скорости с учетом поправки. ьу - постоянная скорость подъема пузырька. их - продольная скорость жидкости. у - поперечная (продольная) координаты; параметр. х поперечная (продольная) координаты; параметр. 2 - поперечная (продольная) координаты; параметр. т - параметр; массы измерительных рамок датчиков; относительная расходная концентрация циркулирующих частиц; (мг/см3) массовая концентрация.

Ь - максимальное смещение пузырька по продольной координате или длина опытного участка (диаметр), параметр. и* - максимальное значение скорости жидкости. п - численное значение показателя степени для пуазейлевского профиля скорости, количество шариков, укладывающихся на одной из сторон опытного участка, показатель преломления, г - рефракция. w - радиус лазерного луча в измерительном объеме, гго = (n —nl)/n - относительное отличие показателей преломления стеклянных шариков и жидкости.

Re3 = 4и/иа - эквивалентное число Рейнольдса; а = ао(1 — s)+3/D - удельная поверхность; ао = 6[2/3 + d/(SH)]/d - удельная поверхность таблетки.

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

Uuz - расходная скорость, приведенная к узкому сечению кубической упаковки.

U* - максимальное значение скорости в центральной ячейке кубической упаковки.

Ра - акустическое давление. Р - давление.

К - относительная мощность звука. и' -амплитуда пульсаций скорости. о - предельное значение акустической интенсивности. ой - циклическая частота пульсаций. t - время.

Ra - критерий Релея. - коэффициент теплового расширения жидкости.

AT - перепад температур.

Рг - критерий Прандтля.

Ra\ - критерий Релея для зернистого слоя.

Л - коэффициент теплопроводности жидкости; газа.

Ао - эффективный коэффициент теплопроводности жидкости.

I - динамическая вязкость жидкости. do - эквивалентный диаметр зернистого слоя. ф - минимальное относительное проходное сечение засыпки. щ, щ - скорости потока с двух сторон щели. и[, и'2 - скорости вихрей.

Z, h - ширина щели и толщина пластины.

Wi(l).W&(b) - силы трения на единицу длины I, Ь. pi, р2 - параметры. gi, q2 - параметры.

Gi(l).Gs{b) - диффузионный поток на единицу длины /, Ь.

Ci, С2 - концентрация в потоке с двух сторон щели. и[, и'2 - концентрация в вихрях.

3(1) - средний коэффициент массопереноса.

Re\ = u\l/v] Re2 = U2I/V.

Jb hi h ~ потоки массы через соответствующие слои смешения. 0"ъ ^г? - толщины турбулентных слоев смешения. Di, D2, D^ - эффективные коэффициенты турбулентной диффузии. Dr - эффективный коэффициент поперечной диффузии. Vot - относительная скорость, равная отношению измеренной скорости к расходной в незаполненном опытном участке. J - поток энтальпии.

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

Ti - температура окончания выхода летучих. АТ - разность температур на входе и выходе потока из датчика. T'd - производная по времени от температуры измерительной рамки датчика.

6 - начальный размер частицы угольной пыли. N и - тепловой критерий Нуссельта.

Список условных обозначений к рисункам рис.1.10, ХА и ВР - хромель-алюмель и вольфрам-рений, рис.1.12, (¿2 и (¿о -расход воздуха через боковые решетки и через центр аппарата. рис.1.15, и мВ - напряжение на выходе пьезоэлемента, А - амплитуда сигнала, Т - длительность сигнала. рис. 1.18, N - количество ударов или отсчетов импульсов с пьезоэлемента в тысячах; и - соответствующее напряжение на выходе в вольтах. рис.1.21, С} (см3/с) - расход жидкости, / • 10~4 (Гц) - частота на выходе измерительного контура. рис.1.22, 0 (см3/с) - расход жидкости, N - показания на шкале ротаметра РС-5. рис.1.29, и/и* - отношение локальной скорости к максимальной в цилиндрической стеклянной трубе, у - поперечная координата, г -радиус трубы, п - показатель преломления жидкости, Ь - длина измерительного объема лазерного доплеровского анемометра, А - центр измерительного объема. рис. 1.30, Ь/с1 - локальный линейный просвет, отнесенный к диаметру шарика; у - поперечная координата в хаотичной упаковке, диаметр которой в 7.5 раза больше диаметра шарика. рис. 1.34, V, V* - локальная и максимальные скорости в канале, х и Ь - продольное и максимальное смещения пузырьков по потоку, у -поперечная координата, II - радиус канала. рис.1.36, рис. 1.37, т/отн кинематическая вязкость раствора, отнесенная к вязкости воды. рис.1.37, х - смещение лазерного луча, при попадании на поверхность шарика, у - соответствующая поперечная координата, II - радиус шарика. рис.2.4г, Ат - амплитуда изменения перепада давления на расстоянии 2 см, отнесенная к среднему перепаду давления на этом расстоянии. рис.2.9-рис.2.18, - продольная скорость, отнесенная к расходной, х,у - координаты в плоскости поперечной внешнему градиенту давления. рис.2.19-рис.2.22, Усо! - поперечная скорость, отнесенная к расходной, х,у - координаты в плоскости поперечной внешнему градиенту давления. рис.2.23-рис.2.32, Уо1 - продольная скорость, отнесенная к расходной в хаотичной упаковке, диаметр которой в 7.5 раза больше диаметра шарика. рис.2.33, Усо! - поперечная скорость, отнесенная к расходной в кубической упаковке. рис.2.34, Усо! - поперечная скорость, отнесенная к расходной в отдельных точках кубической упаковки, t - время измерения в минутах. рис.2.35, относительные спектры продольных пульсаций скорости в правильных упаковках в зависимости от числа Струхаля, равному произведению частоты на диаметр шарика, поделенному на скорость потока в точке измерения. рис.3.1, щ, и2 - скорости потока с двух сторон щели, и'1; и2 - скорости вихрей в щели, с\, с2 - концентрации с двух сторон щели, с'ь с'2

- концентрации внутри щели; а, Ь - размеры вихрей, 1, Ь - толщина пластины и длина щели. рис.3.2, Б - коэффициент диффузии в жидкости, 11е и Эс - числа Рейнольдса и Шмидта. рис.4.5-4.7, Уо1 - продольная скорость, отнесенная к расходной, Б

- диаметр засыпки с таблетками. рис.5.3, Н - магнитное поле на оси соленоида при разной коммутадии катушек, II - постоянное напряжение, II - сопротивление катушки. рис.6.1, Ь, кг/кг -расход топлива, приходящийся на единицу расхода воздуха. рис.6.5, уг - вертикальная составляющая скорости запыленного потока. рис.6.6, - обозначения приведены на рис.1.12.

Список литературы диссертационного исследования доктор технических наук Волков, Валерий Иванович, 1998 год

Литература

[1] Аэров M.Э.,Тодес ОМ. Гидравлические и тепловые основы работы аппаратов со стационарным и кипящим зернистым слоем. -Л.:Химия, 1968. -512 с.

[2] Аэров М.Э.,Тодес О.М., H арийский Д. А. Аппараты со стационарным и кипящим зернистым слоем. -Л.:Химия, 1979. -176 с.

[3] Накоряков В.Е. Исследование турбулентных течений двухфазных сред. - Новосибирск: ИТФ СО АН СССР, 1973. -152 с.

[4] Струминский В.В., Павлихина М.Я. Экспериментальное исследование поля скоростей за стационарным зернистым слоем катализатора. -В сборнике " Аэродинамика в технологических процессах". М.: Наука, 1981, с.63-74.

[5] Филиппов В.М. Измерение поля потока за неподвижным зернистым слоем. -В сборнике "Аэродинамика в технологических процессах". М.: Наука, 1981, с.13-29.

[6] Кириллов В.А., Кузьмин В.А., Пьяное В.И., Канаев В.М. О профиле скорости в неподвижном зернистом слое. Доклады АН СССР, 1979, том 245, N 1, с.159-162.

[7] Lerou J.J. and Froment G.F. Velocity, temperature and conversion profiles in fixed bed catalytic reactors. - Chem. Engng. Sei., 1977, v.32, N 8, p.853-861.

[8] Van der Merve D.F. and Gauvin W.H. Velocity and turbulence measurements of air flow througth a packed bed. -AIChE journal, 1971, v.17, N 3, p.519-528.

[9] Mickley H.S., Smith K.A. and Kochak E.I. Fluid flow in packed beds. - Chem. Engng. Sci., 1965, v.20. N 3, p.237-246.

[10] Евсеев A.P., Накоряков B.E., Романов H.H. Локальная структура фильтрационного потока в кубической упаковке шаров при больших числах Рейнольдса. - Изв. СО АН СССР. Сер. техн. наук. - 1989. - Вып.1, с.51-56.

[11] Yevseyev A.R., Nakorykov V.E. and Romanov N.N. Experimental investigation of a turbulent filtrational flow. -Int. J. Multiphase Flow, 1991, v.17, N 1, p.103-118

[12] Akehato F., Sato K. Flow distribution in packed beds. Kagaku Kogaku (Chem. Eng., Japan.), 1958, v.22, N 7, p.430-436.

[13] Kubota H., Ikeda M., Nishimura V. Note on flow-profile in packed beds. Kagaku Kogaku (Chem. Eng., Japan.), 1966, v.4, N 1, p.58-61.

[14] Cairns E.I. and Prausnitz I.M. Velocity profiles in packed and fluidized beds. -Ind. and Engng. Chem., 1959, v.51, N 12, p.1441-1444.

[15] Колесанов Ф.Ф. Движение газа через слой кусковых материалов. -М.¡Издательство литературы по черной и цветной металлургии, 1956, с.43-65.

[16] Табунщиков Н. П. Поле скорости газа в шахтных известковых печах. ЖПХ, 1956, том 29, N1, с.32-40.

[17] Краснушкина Н.В., Дильмап В.В., Сергеев СЛ. Влияние профиля скорости потока в слое катализатора на производительность реактора синтеза аммиака. -В сборнике " Аэродинамика в технологических процессах". М.: Наука, 1981, с.74-79.

[18] Абаев Г.Н., Попов Е.К. и др. Результаты исследования аэродинамики зернистого слоя на стендах и промышленных реакторах синтеза мономеров для синтетического каучука. -В сборнике "Аэродинамика в технологических процессах". М.: Наука, 1981, с.79-91.

[19] Беляшевский H.H., Бугай Н.Г. Гидравлическая структура потока в отдельной поре при нелинейной фильтрации. -В сборнике "Фильтрация воды в пористых средах". Киев, 1978, с.15-24.

[20] Шейдеггер А.Э. Физика течения жидкости через пористые среды. -М.:Гостоптехиздат, 1960. с.38-42.

[21] Коллинз Р.Течение жидкости через пористые материалы. -М.;Мир, 1964. с.277.

[22] Christiunsen С. Untersuchungen ü her die optischen Eigenschaften von fein vertheilten К Ö rpern. Ann. Phisik. Cherri., 1884, v. 23, p.298-310.

[23] Johnston W., Dybbs A., Edvards R. Measurement of fluid velocity inside porous media with a laser anemometer, he physics of fluids. 1975, v. 18, N 7, p.913-914.

[24] Бэтчелор Дж. Введение в динамику жидкости. ~М.:Мир, 1973. с.329.

[25] Аэров М.Э., Умник H.H. Коэффициенты теплопроводности в зернистом слое. ЖТФ, 1951, том 21, N И, с.1351-1352.

[26] Голубев B.C. Гидродинамическая дисперсия и динамика сорбции в пористой среде с застойными зонами. Доклады АН СССР, 1978, том 23, N 25, с.1161-1162.

[27] Николаевский В.Н. Капиллярная модель диффузии в пористых средах. - Изв. АН СССР, 1959, N 4, с.146-147.

[28] Иоффе И.И., Письмен JI.M. Инженерная химия гетерогенного катализа. -Л.:Химия, 1972. с.39.

[29] Касаманян М.А., Кириллов В.А., Матрос Ю.М. Перенос вещества в свободном объеме неподвижного зернистого слоя, зернистом слое. ИФЖ, 1973, том 25, N 1, с.36-37.

[30] Левин В.Г., Маркин B.C., Чизмаджиев Ю.А. О гидродинамическом перемешивании в модели пористой среды с застойными зонами. Доклады АН СССР, 1966, том 166, N 6, с.1401-1402.

[31] Левин В.Г., Письмен Л.М., Кучанов С.Н. О гидродинамическом перемешивании в зернистом слое. Физическая модель застойных зон. Доклады АН СССР, 1966, том 168, N 2, с.392-393.

[32] Письмен Л.М., Кучанов С.Н., Левич В. Г. Поперечная диффузия в зернистом слое. Доклады АН СССР, 1967, том 174, N 3, с.650-651.

[33] Кучанов С.Н., Левич В.Г., Письмен Л.М. Поперечная диффузия и теплопроводность в зернистом слое. ПМТФ, 1967, N 3, с.45-46.

[34] Prausnitz I.M. and Wilhelm R.H. Turbulent concentration fluctuations in a packed bed. -Ind. and Engng. Chem., 1957, v.49, N 6, p.978-984.

[35] Волков В. И. Изобретательские задачи в процессах переноса. Барнаул. ¡Издательство Алтайского госуниверситета, 1997, -150 с.

[36] Волков В.И., Мухин В.А., Накоряков В.Е. Исследование структуры течения в пористой среде. ЖПХ, 1981, том 34, N4, с.838-842.

[37] Трембовля В. И. и др. Теплотехнические испытания котельных установок. -М.:Энергия, 1977, -296 с.

[38] Волков В.И., Колмогоров А.И., Кисляк С.М. Прибор для определения параметров двухфазного потока. Заводская лаборатория. Диагностика материалов. 1998, том 64, N 2, с.39-40.

[39] Малое 5.5.Пьезорезонансные датчики.-М.:Энергоатомиздат, 1989, -272 с.

[40] Кремлевский 77. Я. Расходомеры и счетчики количества. Справочник -JI.: Машиностроение, 1989, -701 с.

[41] Волков В.И., Волков М.Ю., Кисляк С.М. Измеритель скорости запыленного потока. Сборник трудов Всероссийской научно-технической конференции: Экспериментальные методы в физике структурно неоднородных сред. Издательство Алтайского госуниверситета, Барнаул, 1997, том 4, с.10-12.

[42] Волков В.И., Волков М.Ю., Максимова С.Н. Устройство для измерения концентрации твердой фазы двухфазного потока. Пат.2059223, 1996, бюлл.12.

[43] Волков В.И., Колмогоров А.И., Соболев М.А. Устройство для определения характеристик двухфазного потока. Пат.2059222, 1996, бюлл.12.

[44] Волков В.И., Колмогоров А. И. Калибратор скорости частиц.

Пат. 2082175, 1997, бюлл.17.

[45] Волков В.И., Колмогоров А.И., Кисляк С.М. Определение параметров твердой фазы. Сборник трудов Всероссийской научно-технической конференции: Экспериментальные методы в физике структурно неоднородных сред. Издательство Алтайского госуниверситета, Барнаул, 1998, том 1, с.66-69.

[46] Волков В.И., Волков М.Ю., Славская Е.Г. Ротаметр. А.С.1682793, 1991, бюлл.37.

[47] Волков В.И., Волков М.Ю., Губенко С.М., Утемесов М.А. Электромагнитный объемный насос. А.С.1536035, 1990, бюлл.2.

[48] Волков В.И., Волков М.Ю. Электромагнитные преобразователи линейных перемещений. Сборник тезисов докладов Региональной научно-технической конференции: Порошковые материалы и покрытия. Издательство Алтайского госуниверситета, Барнаул, 1990, с.104-106.

[49] Волков В.И. Расходомер. Пат. 2104496, 1998, бюлл.4.

[50] Волков В. И. Изотермическое течение жидкости в упаковке их сфер. ИФЖ, 1985, т.ХЫХ, N5 с.827-833.

[51] Волков В.И., Мухин В.А., Титков В.И. Измерение пористости хаотичной упаковки с помощью ЛДИС. Сборник трудов Всероссийской научно-технической конференции: Экспериментальные методы в физике структурно неоднородных сред. Издательство Алтайского госуниверситета, Барнаул, 1998, том 1, с.69-71.

[52] Химические реактивы и высокочистые химические вещества Каталог. -М.: Издательство "Химия", 1990, с.214.

[53] Татарский В. Б. Кристаллооптика и иммерсионный метод. -М.:Наука, 1965, -125 с.

[54] Иоффе Б.В. Рефрактометрические методы химии. -Л.: Химия, 1974, -400 с.

[55] Волков В.И:, Мухин В.А. Жидкость для оптических исследований. А.С.948994, 1982, бюлл.29.

[56] Справочник химика. /Ред. Б.Л.Никольский: Т.2. -Л.: Химия, 1971.-851 с.

[57] Анисимов К.Г., Анисимова Е.А., Волков В.И. Экспериментальное исследование показателя и характеристики переноса тяжелых жидкостей. Известия АГУ, N 1, 1997, с.58-59.

[58] Волков В.И., Анисимов К.Г., Перов Э.И. Исследование гидролиза высокопреломляющих сред. Сборник трудов Всероссийской научно-технической конференции: Экспериментальные методы в физике структурно-неоднородных сред. Издательство Алтайского госуниверситета, Барнаул, 1998, том 1, с.71-74.

[59] Волков В.И., Волков М.Ю., и др. Иммерсионная жидкость. Пат. 2051940, 1996, бюлл.1.

[60] Россоти Ф., Россотти X. Определение констант устойчивости и других констант равновесия в растворах. -М.: Мир, 1965. -210 с.

[61] Гатчек Э. Вязкость жидкостей. -М.,Л.: ГТТИ, 1932, -215 с.

[62] Эрдеи-Груз Т. Явление переноса в водных растворах. -М.: Мир, 1976. -597 с.

[63] Анисимов К.Г., Волков В.И. Экспериментальное исследование пристеночного перепада давления при переходе к турбулент-

ности. Сб. трудов V Международного семинара: Устойчивость гомогенных и гетерогенных жидкостей, часть I, 1998, с. 11-15.

[64] Волков В.И., Волков М.Ю. Исследование перепада давления в зернистом слое. ЖПХ. Деп.К5654-1988, 1988 .

[65] Анисимов К.Г., Волков В.И. Экспериментальное исследование перепада давления в системе связанных вихрей в зернистом слое. Известия АГУ, N 1, 1998, с.12-13.

[66] Воришанский В.М., Виноградов A.C., Лузин И.П., ТигаревИ.П. Гидравлическое сопротивление засыпок из сферических частиц. Теплоэнергетика, 1980, N 1, с.61-63.

[67] Волков В. И. Исследование гидродинамики и процессов переноса в пористых средах. Канд. дисс., Институт теплофизики СО АН СССР, -Новосибирск, 1980, с.156

[68] Волков В. И. Исследование гидродинамики и процессов переноса в пористых средах: Автореф. дис. ... канд. техн. наук. - Новосибирск: ИТФ СО АН СССР, 1982. -21 с.

[69] Волков В.И., Мухин В.А., Накоряков В.Е. и др. Исследование гидродинамики пристенного слоя на модели кубической упаковки. ПМТФ, 1980, N6, с.58 - 64.

[70] Волков В.И., Накоряков В.Е., Мухин В.А. Исследование поля скоростей в кубической упаковке с помощью лазерного до-плеровского анемометра. -В сб.тезисов докладов 2-го Всесоюзного семинара "Теплофизика гидроохлаждаемых реакторов." МВТУ им. Баумана и ин-т Атомной энергии им. Курчатова. Москву, 1980, с.7.

[71] Волков В. И. Исследование гидродинамики в пористой среде с помощью ЛДИС. -В сб.: Физическая гидродинамика и тепловые процессы. Новосибирск, 1980, с.5-9.

[72] Волков В.И., Накоряков В.Е., Мухин В.А. Исследование гидродинамики на модели пористого пласта с помощью лазерной анемометрии. -В сб.тезисов докладов 2-го Всесоюзной научно-технической конференции "Проблемы горной теплофизики". Ленинградский Горный институт. Ленинград 1981, с. 11.

[73] Волков В.И., Накоряков В.^. Экспериментальное исследование поля скоростей в неподвижном зернистом слое. -В сб. трудов: Процессы переноса в аппаратах энергохимических производств. Институт Теплофизики, Новосибирск, 1985, с.15-23.

[74] Волков В.И., Мухин В.А., Титков В.И., Томсонс Я.Я. Исследование поля скоростей в пористой среде с помощью лазерного доплеровского измерителя скорости. Автометрия, 1982, N3, с.82 - 86.

[75] Волков В.И., Мухин В.А. Исследование структуры течения в пористой среде с помощью лазерной анемометрии. -В сб.: Аннотация докладов V Всесоюзного съезда по теоретической и прикладной механике. Алма-Ата, 1981, с.96.

[76] Голъдштик М.А. Теория концентрированных дисперсных систем. В сб.: Материалы международной школы по процессам переноса в неподвижных и псевдоожиженных зернистых слоях. Минск, 1977, с.49-84.

[77] Карпенко Ю.В., Коротченко Г.И., Ноздрин Г.Н., Фомичев М.С. Экспериментальное исследование колебаний давления жидкости в циркуляционном контуре при работе центробежного насоса. Теплоэнергетика, 1987, N 1, с.69-71

[78] Римский-Корсаков A.B., Баженов Д.В., Баженова Л.А. Физические основы образования звука в воздуходувных машинах. -М.: Наука, 1988. -173с.

[79] Крайнов В. П. Качественные методы в физической кинетике и

гидрогазодинамике. -М.: Высшая школа, 1989, с.202.

[80] Волков В.И., Мухин В.А., Типаков В.И. Неустойчивость течения жидкости в засыпках. Сб. трудов V Международного семинара: Устойчивость гомогенных и гетерогенных жидкостей, часть I, 1998, с.43-49.

[81] Волков В.И., Волков М.Ю. Способ управления тепломассооб-менным аппаратом. A.C. 1286265, 1987, бюлл.4.

[82] Волков В.И., Волков М.Ю. Способ управления тепломассооб-менным аппаратом. А.С.1567261, 1990, бюлл.20.

[83] Шец Дж. Турбулентное течение. Процессы вдува и перемешивания. -М.: Мир, 1984, с.46.

[84] Ишии Ч., Хонами С. Нестационарные структуры отрывного трехмерного турбулентного пограничного слоя. Сб. работ под ред. Фернхольца X. и Краузе Е. Трехмерные турбулентные пограничные слои.// -М.: Мир, 1985, с.117.

[85] Монин A.C., Яглом A.M. Статистическая гидромеханика. Часть 2. -М.: Наука, 1967, с.441.

[86] Кадер Б.А., Яглом A.M. Спектры анизотропных турбулентных пульсаций скорости и температуры в пристеночных турбулентных потоках.// Проблемы турбулентных течений. -М.: Наука, 1987, С.70.

[87] Гершуни Г.З., Жуховицкий Е.М. Конвективная устойчивость несжимаемой жидкости. - М.: Наука, 1972, с.74.

[88] Brinkman Я. С. A calculation of the viscous force exerted by a flowing fluid on a dense swarm of particles. -Appl. Sci. Res., 1947,

Vol.Al, pp. 27-34.

[89] Бадатов E.B., Слинько М.Г., Накоряков B.E. Математическое моделирование процессов переноса в отрывных течениях с ламинарной областью смешения. Теоретические основы химической технологии. 1970, том 4, N5, с.687 - 693.

[90] Бадатов Е.В., Слинько М.Г., Накоряков В.Е. Математическое моделирование процессов переноса в отрывных течениях с турбулентной областью смешения. Теоретические основы химической технологии. 1970, том 4, N6, с.868 - 864.

[91] Волков В.И., Сеначин П.К., Утемесов М.А., Потапов М.Г. и др. Исследование теплофизических параметров двухфазного потока методом тонкой стенки. Сб.трудов Всесоюзной конференции "Актуальные вопросы физики аэродисперсных систем", Изд-во ОГУ, Одесса, 1986, т.2., с.23.

[92] Volkov V.I. Pressure drop distribution and redistribution in packed columns. Workshop, Transport phenomena in twophase flou. 1997. c.37.

[93] Волков В.И., Хрусталев Ю.В. Исследование теплового потока в разрушающихся упаковках. Теплофизика и аэромеханика. 1998, том 5, N 1, с.125-128.

[94] Волков М.Ю., Сакович Г.Г., Волков В.И., Шандаков В.А. Генератор газа. А.С.1517997, 1989, бюлл.40.

[95] Волков В.И., Волков М.Ю., Губенко С.М., Утемесов М.А. Генератор газа. А.С.1498548, 1989, бюлл.29.

[96] Волков В.И., Волков М.Ю. Устройство для получения газа. А.С.1637865, 1991, бюлл.12.

[97] Волков В.И., Волков М.Ю., Утемесов М.А., Бобыръ П.Н. Устройство для получения газа. А.С.1538051, 1990, бюлл.З.

[98] Волков В.И., Волков М.Ю.. Утемесов М.А., Шандаков В.А., Сакович Г.Г. Газогенератор газа. А.С.1524912, 1989, бюлл.44.

[99] Волков В.И., Волков М.Ю. Газогенератор. А.С.1674975, 1991, бюлл.ЗЗ.

[100] Волков М.Ю., Волков В.И., Шечков Г. Т. Применение микрока-псулирование в химических генераторах газа. Сборник трудов Всероссийской научно-технической конференции: Экспериментальные методы в физике структурно неоднородных сред. Издательство Алтайского госуниверситета, Барнаул, 1997, том 4, с.12-14.

[101] Бурдуков А.П., Накоряков В.Е.. Горение угольной пылевзвеси в ультразвуковом поле. Сб.трудов "Кинетика и горение ископаемых топлив". Издательство института Теплофизики СО АН СССР, Новосибирск, 1963, с.5-7.

[102] Бурдуков А.П., Накоряков В.Е.. Теплообмен от цилиндра в звуковом поле при числах Грасгофа, стремящихся к нулю. ПМТФ, N 1, 1965, с.7-10.

[103] Бурдуков А.П., Накоряков В.Е.. О переносе массы в звуковом поле. ПМТФ, N 2, 1965, с.11-14.

[104] Бурдуков А.П., Накоряков В.Е.. Влияние колебаний на массо-отдачу от сферы при больших числах Прандтля. ПМТФ, N 3, 1967, с.23-27.

[105] Бурдуков А.П., Зауличный В.Г., Накоряков В.Е.. Теплооотдача от цилиндра в звуковом поле. Известия Сибирского отделения АН СССР N 6, вып. 2, 1965, с.18-20.

[106] Хзмалян Д.М., Каган Я.А. Теория горения и топочные устройства. -М.: Энергия, 1967, с.355.

[107] Сеначин П.К.,Волков В.И., Утемесов М.А. Влияние пульсаций воздуха на горение пыли натурального твердого топлива. Тезисы докладов II всесоюзного научно-технического семинара "Сжигание топлив с минимальными вредными выбросами". Академия наук Эстонской ССР, Таллин, 1978, с.42.

[108] Волков В.И., Сеначин П.К. Управление процессами зажигания и горения пылевидного твердого топлива в парогенераторах. Отчет по научно - исследовательской работе. Алтайский государственный университет, N 77076475, инв.Б680292, от 11 июля 1978, 76 с.

[109] Мальцев A.A., Бажов В.И., Павлов Я.В., Пилягин В.Ф., Jlumu-нецкий В.Я., Шингиль И.А., Волков В.И., Сеначин П.К., Утемесов М.А. Способ сжигания твердого топлива. A.C. 868256, 30.09. 1981, бюлл. 36.

[110] Волков В.И., Накоряков В.Е. Способ термического воздействия на нефтяной пласт. А.С.928858, 1982, бюлл.18.

[111] Мальцев A.A., Бажов В.И., Павлов Н.В., Пилягин В.Ф., Лити-нецкий В.Я., Шингиль И.А., Гордеев В.Г., Сеначин П.К., Волков В.И., Утемесов М.А. Устройство для сжигания топлива. А.С.857642, 1981, бюлл.31.

[112] Волков В.И., Волков М.Ю. Пульсатор потока тепломассообмен-ного аппарата. A.C. 1182249, 1985, бюлл. 36.

[ИЗ] Волков В.И., Волков М.Ю., Губенко С.М., Утемесов М.А. Теплообменник с регулируемым теплосъемом. A.C. 1483234, 1989, бюлл.20.

[114] Волков В.И., Вобырь П.Н., Кисляк С.М., Потапов М.Г., Лушев В.П. Хрусталев Ю.В. Коллектор теплообменника. A.C.1276901, 1986, бюлл.46.

[115] Алексеенко C.B., Срывков C.B., Ефименко A.M. и др. Изотермическое моделирование фонтанно - вихревой топки. -Электрические станции, 1992, N 11, с.20 - 25.

[116] Волков В.И., Волков М.Ю., Кисляк С.М. Исследование аэродинамики в кипящем слое. Сборник трудов Всероссийской научно-технической конференции: Экспериментальные методы в физике структурно-неоднородных сред. Издательство Алтайского госуниверситета, Барнаул, 1997, том 4, с.8-10.

[117] Волков В.И., Кисляк С.М. Моделирование аэродинамики топочных устройств с развитой внутритопочной циркуляцией частиц. Известия АГУ, 1997, N 1, с.44-46.

[118] Волков В.И., Кисляк С.М., Колмогоров А.И. Экспериментальное исследование двухфазного потока в циркулирующем кипящем слое. Сборник трудов Всероссийской научно-технической конференции: Экспериментальные методы в физике структурно-неоднородных сред. Издательство Алтайского госуниверситета, Барнаул, 1997, том 4, с.5-7.

[119] Волков В.И., Рунг В.Ю. Управление теплообменом в кипящем зернистом слое. Сборник трудов Всероссийской научно-технической конференции: Экспериментальные методы в физике структурно-неоднородных сред. Издательство Алтайского госуниверситета, Барнаул, 1997, том 4, с.58-60.

[120] Алексеенко C.B., Процайло М.Я., Юрлагин A.B. Исследование турбулентного течения жидкости в прямоугольном канале с боковыми входами. -Известия СО АН СССР, Сер. техн. наук, 1988, N 7, Вып.2

[121] Рыжаков A.B. Исследование гидродинамических и фазо-переходных процессов в установках с кипящим слоем. -Электрические станции, 1996, N 8, с.56 - 61.

[122] Прасолов P.C. Массо и теплоперенос в топочных устройствах. -М.: Энергия, 1964, с.18.

[123] Коварский Л. Г. Защита паровых котлов от шлакования и заноса золой. - М.: Энергия, 1964, с.37.

[124] Кузнецов Н.В, Лужнов Г.И., Кропп Л.И. Очистка поверхностей нагрева котельных агрегатов. - М.: Энергия, 1966, с.15.

[125] Вдовенко М.И. Минеральная часть энергетических углей. -Алма-Ата: Наука, 1973, с.8.

[126] Вухман C.B., Крылова Н.П. Вспучивание частиц натуральных углей в процессе термической обработки. Сб.трудов "Проблемы теплоэнергетики и прикладной теплофизики". Выпуск N 8, Алма-Ата, Наука, 1972, с.6-9.

[127] Вухман C.B., Волков В.И., Сеначин П.К. Исследование возможности шлакования поверхностей в камере сгорания при низких температурах. -Отчет по научно - исследовательской работе. Алтайский государственный университет, N 75034409, инв.Б616172, от 12 августа 1977, 45 с.

[128] Зимон А.Д. Адгезия пыли и порошков. - М.: Химия, 1976, 432 с.

[129] Волков В.И. Способ очистки поверхности нагрева. A.C.634080, 1978, бюлл.43.

[130] Волков В.И., Сеначин П.К. Устройство для сброса шлаковых отложений с поверхности нагрева. А.С.642599, 1979, бюлл.2.

[131] Волков В.И., Сеначин П.К. Способ удаления шлака с поверхности нагрева. А.С.652430, 1979, бюлл.10.

[132] Волков В. И. Устройство для вибрационной очистки труб котло-агрегатов. А.С.798473, 1981, бюлл.З.

[133] Волков В. И. Устройство для воздействия на нефтяной пласт. А.С.733333, 1980, бюлл.9.

[134] Волков В.И. Способ обработки нефтяного пласта. А.С.753188, 1980, бюлл.З.

[135] Волков В.И., Накоряков В.Е., Шрейбер И. Р., Мукук К.В. Устройство для воздействие на нефтяной пласт. А.С.972891, 1982, бюлл.20.

[136] Волков В.И., Накоряков В.Е. Способ и устройство для обработки нефтяных скважин. А.С.927983, 1982, бюлл.18.

[137] Валаклеевский В.И., Волков В.И., Накоряков В.Е., Шрейбер И.Р. Способ обработки нефтяной скважины. А.С.979261, 1983, бюлл.17.

[138] Волков В.И., Накоряков В.Е., Шрейбер И.Р. Способ очистки внутренней поверхности трубопровода. А.С.918773, 1982, бюлл.13.

[139] Волков В.И., Накоряков В.Е., Шрейбер И.Р., Мукук К.В., Кузнецов В. В. Устройство для вибрационной очистки груб. А.С.926509, 1982, бюлл.17.