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

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

Миньков Леонид Леонидович

ИССЛЕДОВАНИЕ ПРОЦЕССОВ РАЗДЕЛЕНИЯ СУСПЕНЗИИ В ЦЕНТРОБЕЖНЫХ УСТРОЙСТВАХ

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

АВТОРЕФЕРАТ

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

Томск - 2012

Работа выполнена на кафедре математической физики ФГБОУ ВПО Национальный исследовательский Томский государственный университет

Научный консультант: доктор физико-математических наук, профессор Шрагер Эрнст Рафаилович

Официальные оппоненты: доктор физико-математических наук, академик РАН Липанов Алексей Матвеевич доктор физико-математических наук, профессор Воеводин Анатолий Федорович доктор физико-математических наук, профессор Бубенчиков Алексей Михайлович

Ведущая организация: ФГБОУ ВПО Национальный исследовательский Томский политехнический университет

Защита диссертации состоится л28 декабря 2012 г. в 10 ч 30 мин на заседании диссертационного совета Д 212.267.13 при Томском государственном университете по адресу: 634050, г. Томск, пр. Ленина, 36.

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

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

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

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

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

Исследование процессов гидроциклонирования активно проводилось в СССР и России Поваровым А.И., Акоповым М.Г., Найденко В.В., Кутеповым А.М., Терновским И.Г., Непомнящим Е.А., Адельшиным А.Б. и др. За рубежом данное направление развивали Шуберт Г., Неессе Т., Хайсканен К.

и др.

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

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

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

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

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

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

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

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

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

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

5. Выявление формы сепарационной кривой классификационного аппарата в области мелких частиц.

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

7. Анализ воздействия инжекционного потока на характеристики классификационного аппарата.

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

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

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

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

4. Объяснено аномальное поведение сепарационной кривой классификационного аппарата в области мелкоразмерных фракций.

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

6. Установлено влияние расхода инжекционного потока, способа инжекции на сепарационную функцию классификационного аппарата.

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

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

Результаты исследования были использованы при моделировании вихревого сепаратора для улавлдивания паров фторуглеродной смазки в технологическом процессе ЗРИ ОАО СХК г. Северск.

Основные результаты диссертации получены при проведении исследований по гранту Немецкой службы академических обменов (2002 г.), гранту Баварского исследовательского фонда (2003Ц2004 гг, PIZ 19/03), гранту Министерства образования Германии (2008 г, BMBF, Project № RUS 08/A01), гранту ФЦП Научные и научно-педагогические кадры иновационной России (2012 г., № 14.B37.21.0872), а также в соответствие с приоритетным направлением инновационной образовательной программой Индустрия наносистем и материалов (2007 г.) и частично по грантам РФФИ №02-01-01022-а, №05-08-01396-а, №08-08-12029-офи, №11-08-00370а.

Положения, выносимые на защиту.

1. Математическая модель увлечения мелких частиц крупными при их совместном оседании.

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

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

4. Результаты численного моделирования оседания частиц полидисперсной суспензии в тарельчатой центрифуге.

5. Объяснение fish-hook эффекта сепарационной кривой классификационного аппарата в области мелких частиц.

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

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

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

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

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

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

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

а) Международных:

4th International Conference on Multiphase flow (May 27 - June 1, New Orlean, LA, USA, 2001);лScience & Technology of Filtration and Separations for the 21st Century (Tampa, Florida, USA,1Ц4 May, 2001);Байкальские чтения-II по моделированию процессов в синергетических системах (Улан-УдэЦТомск, Россия, 18-23 июля 2002); л5-th International Conference on Multiphase Flow (Yokohama, Japan, May 30-June 4, 2004); л5-th International conference on Transport Phenomena in Multiphase Systems (HEAT 2008, Bialystok, Poland, June 30 - July 3, 2008); International Сonference on Physical SeparationТ09 (Falmouth, UK, June 16-17, 2009);л7-th World Conference on Experimental Heat Transfer, Fluid Mechanics and Thermodynamics (ExHFT-7, Krakow, Poland, June - 03 July 2009); л19-th International Congress of Chemical and Process Engineering (CHISA 2010, Prague, Czech Republic, 28 August - 1 September 2010);л21-st International Symposium on Transport Phenomena (Kaohsiung City, Taiwan, November 2-5 2010);лInternational Conference on Fluid dynamics and Thermodynamics (WASET. Amsterdam, Netherlands 13-15 July 2011);

б) Всероссийских:

Байкальские чтения по математическому моделированию в синергетических системах (Улан-УдэЦТомск, 10-23 Июля, 1999); IX Всероссийская научнотехническая конференция Физика и химия высокоэнергетических систем (Томск, 2003); IV, V, VI Всероссийская научная конференция Фундаментальные и прикладные проблемы современной механики (Томск, 5Ц7 октября 2004, 3Ц5 октября 2006, 30 сентября - 2 октября, 2008); XI Всероссийская научно-техническая конференция Физика и химия высокоэнергетических систем (Томск, 2005); VI Всероссийская научная конференция Современная баллистика и смежные вопросы механики (Томск 17 - 19 ноября, 2009); Всероссийская молодёжная научная конференция Актуальные проблемы современной механики сплошных сред (Томск, 16Ц19 октября 2010).

Результаты диссертационной работы докладывались устно на научных семинарах технического факультета Университета Эрланген-Нюрнберг (г.

Эрланген, Германия, в 2003 и 2008 гг.).

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

Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения и списка используемой литературы, включающего 2наименований, содержит 97 рисунка, 5 таблиц - всего 275 страниц.

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

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

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

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

Мелкие частицы вне этого объема V - Vc, не испытывая влияния крупной, оседают с их собственной скоростью uf. Тогда средняя скорость оседания Vc Vc мелких частиц в выделенном объеме V, равна us uc uf 1 .

VV Vc dc Учитывая, что f f (c ), где c - объемная доля крупной V V фракции частиц в суспензии, а скорости оседания частиц пропрциональны квадратам их диаметров, получается:

us dc 1 f (c ) f (c ). (1) uff d Соотношение (1) показывает, что увеличение скорости седиментации мелких частиц в бидисперсной суспензии пропорционально квадрату онтношения диаметров частиц крупной и мелкой фракций.

Вид функции f (c ) можно установить, принимая те или иные гипотезы относительно геометрии ячейки и течения в ней. Если принять ячейку в виде сферы радиусом R, и крупная частица оказывает одинаковое влияние на все мелкие частицы, находящиеся на расстоянии меньшем dc/2 от своей поверхности, то f (c) 7c. Если же принять ячейку в виде цилиндра с радиусом основания R и высотой 2R, а лактивную часть ячейки в виде цилиндра такого же радиуса и высотой dc с выброшенным объемом крупной 1/1 2 1/3 частицы, то f (c ) 3c 2 2c 3.

Во втором параграфе моделируется движение мелкой частицы в поле течения крупной частицы, вызванного ее оседанием, на основании чего выводится формула для скорости оседания мелкой частицы в присутствии крупной. Задача рассматривается в системе координат, связанной с крупной частицей, помещенной в цилиндрическую ячейку радиусом основания R и высотой H=2R. Предполагается, что поле течения вокруг крупной частицы является Стоксовым; мелкие частицы не влияют на поле теченния жидкости (Stf < 1).

Средняя скорость прохождения мелкой частицей ячейки определяется из соотношения up( yj ) H / T( yj ), где T ( yj ) - время прохождения ячейки частицей, имеющей во входном сечении ячейки координату y.

j Среднепотоковая скорость мелких частиц, проходящих через ячейку, имеет R вид: up 2 yjup( yj )dyj, а их средняя скорость седиментации в R2 лабораторной системе координат:

us uc up. (2) Время прохождения мелкой частицы ячейки находится из решения системы кинематических уравнений:

dt 1 dy v(x, y) ; .(3) dx u(x, y) uf dx u(x, y) uf со следующими начальными условиями x R : t 0; y yj. (4) Здесь поле скоростей жидкости опрделяется согласно Г. Стоксу по формулам (Шлихтинг Г. Теория пограничного слоя. - М.: Наука, 1974. - 712 с.):

1 1 3x2 u(x, y) uc 1 3 1 4r r2 4r3rc2 r2 3xy v(x, y) uc 1 ;

4r3rc2 r2 2 2 x y xy r2 ; 1; 1;

rc rc rc rc x R, R.

T( yj ) t(R, y(R)), rc dc 2.

Рис.1. Траектории мелких частиц, D= Анализ размерностей показывает, что решение задачи (2)Ц(4) зависит от двух параметров U uf uc, D R rc f, т.е. us uc Us F(U, D). На рис.1 приведены характерные траектории f мелких частиц, огибающих крупную частицу. Наличие собственной скорости седиментации ведет к тому, что траектории становятся более выгнутыми (штриховые линии), и частицы вынуждены пробегать больший путь при увеличении Uf. При этом увеличивается время прохождения частицей ячейки, что в свою очередь ведет к увеличению Us. Расчеты показывают, что зависимость Us F(U, D) с погрешностью менее одного процента f аппроксимируется формулой Us 1.35 D U 11.2 D2 в диапазоне f значений D > 2 и Uf < 0.3, которая в размерных переменных принимает вид:

us dc 2/3 1/ 11.57c 1.54c (5) uff d Оценка границ применимости модели получается на основе анализа точки останова мелкой частицы в направлении х при y=0: u(x,0) = uf или 3 X 1U . Полагая X D, и учитывая, что для данного типа ячейки f 2 X 1/ D 2 3c, указанное неравенство в размерных переменных запишется в виде:

4/d f 1/ 1 c c (6) dc Если отношение диаметров частиц меньше значения, даваемого правой частью неравенства (6), то мелкая частица проникает в ячейку и может быть увлечена крупной частицей. Увеличение концентрации крупных частиц ведет к снижению размера мелких частиц, увлекаемых крупной.

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

Оценивается время пребывания частицы в ячейке: r 2(1 2), где 1 - время прохождения прямолинейного участка траектории, D X0 U, а 2 - время прохождения половины кругового Рис.2. Схема траектории частицы участка, 1.5ln 2 X Y0. Здесь U - средняя скорость прохождения прямолинейного участка 1U 3 2D 1 Y0 D. Средняя скорость седиментации находится из f D 4 Yсоотношения Us 1 f dY0 U a D, где a [1.5 2,1.5].

D (Y0) 0 r Переходя к размерным переменным, полагая a=1.3 и учитывая, что 1.5D3 1 c, получается соотношение близкое к (5):

us dc 1/ 11.48c (7) uff d Анализ траекторий частиц показывает, что основную роль для эффекта ускорения седиментации играет не та малая доля мелких частиц, лежащих на пути крупной частицы, чьи траектории резко искажаются в ходе ее обтекания, а та основная доля частиц, попадающих в ячейку в периферийной её области и испытывающих относительно слабое увлечение.

Справедливость полученных формул (1), (6) и (7) подтверждается экспериментальными данными по седиментации частиц кварцевого песка суспензии в тарельчатой центрифуге (Ch. Gerhart, J.Dck, Th.Neee.

Grundlagenuntersuchungen zur behinderten Sedimentation polydisperser Suspensionen bei der Hydrostromklassierung. Teil I: Untersuchungen in einer Laborzentrifuge //Aufbercitungs technik. - 1999. - Vol. 40, №7. - P. 328Ц334.) и по всплыванию частиц парафинного масла эмульсии в вертикальной колонне (S.Kumar, Th.W.Pirog, D.Ramkrishna. A new method for estimating hindered creaming/settling velocity of particles in Polydisperse System // Chemical Engineering Science. - 2000. - V.55. - P.1893Ц1904).

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

Рассмотрены два вида ячеек:

цилиндрическая ячейка (рис.3а), для которой справедливо 1.5D3 1 c, и сферическая ячейка Хаппеля радиуса R (рис.3б), для которой D3 1 c. Кроме геометрии, а) б) Рис.3 Влияние силы Саффмэна на траектории ячейки отличаются по мелких частиц. а - цилиндрическая ячейка; б - условиям, заданным на границе:

ячейка Хаппеля.

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

dxp dyp up (xp, yp ); vp(xp, yp ). (8) dt dt В начальной точке траектории частицы задаются условия: xp t0 x0 ;

yp t0 R - для цилиндрической ячейки и yp t0 R2 x0 - для сферической ячейки. Поле скорости частиц находится из решения уравнения баланса сил, действующих на частицу: Fg FA FS FSt 0, где Fg d3S g j 6 - сила тяжести; FA d3l g j 6 - сила Архимеда;

f f 0.FS 1.615ld2 l rot ul ul up rot ul - сила Саффмэна, вызванная f вращением частицы в сдвиговом потоке, FSt 3ld ul up - сила f сопротивления Стокса, действующая со стороны жидкости на частицу.

Система уравнений (8) решается до выполнения одного из условий: а) пересечение траектории мелкой частицы с границей крупной частицы x2 y2 dc 2, кривая 1, рис. 3; б) выход мелкой частицы на предельную p p траекторию, отвечающую условию dyp dxp xp 0 0, кривая 2, рис. 3; в) достижение частицей границы ячейки: yp R или x2 y2 R2, кривая 3, p p рис. 3.

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

us d f 1 1 a b d dc . (9) uf dc f 0.331 0.4Здесь коэффиценты a 4.474c, b 0.972c - для ячейки Хаппеля и 0.449 0.3a 15.971c, b 0.532c - для цилиндрической ячейки получены на основе аппроксимации рассчитанных величин x* R функцией auf b uf с точностью до 1% в диапазоне 2 R rc 10.

Из полученной зависимости (9) следует, что при увеличении размера крупной частицы dc 1.5d средняя скорость оседания мелких частиц us f будет уменьшаться, что противоречит известным экспериментальным данным. Следовательно, силой Саффмэна невозможно объяснить наблюдаемого многократного ускорения седиментации мелких частиц в полидисперсной суспензии.

В пятом параграфе формула вида (7) обобщается на случай полидисперсной суспензии, которая описывается как набором дискретных фракций, имеющих объемную долю j, так и непрерывной функцией распределения частиц по размерм q(d). В первом случае совместное влияние всех крупных фракций на ускоренное оседание мелкой фракции описывается 1/N us зависимостью d 1 1.48 d6 , а во втором случае эта f j j uf d d* j 1/ us зависимость принимает вид: d 1 1.48 q(x)x6dx. Здесь d* - f uf d* наименьший размер частиц крупной фракции, котороые способны увлекать мелкую частицу.

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

В первом параграфе излагается общая феноменологическая теория для разделения полидисперсной суспензии в поле гравитационных и центробежных сил (рис.4) при следующих допущениях: 1) Суспензия состоит из двух несжимаемых фаз: фазы твердых частиц и жидкой фазы; 2) Вся совокупность твердых частиц разбита на N фракций; 3) Размеры частиц намного меньше рассматриваемой области; 4) Перед началом осаждения суспензия равномерно распределена по Рис.4. Схема вращающейся центрифуги всему объему; 5) Массообмен между фазами отсутствует; 6) Все вязкие эффекты смеси опрделяются вязкими эффектами жидкости.

Система уравнений, описывающая разделение суспензии, в системе координат, связанной с центрифугой, согласно работе Берреса и Бюргера (Berres S., Brger R. On gravity and centrifugal settling of polydisperse suspensions forming compressible sediments // Int. J. Solids and Structures. - 2003. ЦV.40. - P. 4965Ц4987), имеет вид:

Уравнение сохранения массы i - ой фракции для твердой фазы i iVi 0, i 1,..., N (10) t Уравнение баланса сохранения смеси:

V 0, (11) N N где V 1 Vf Vj, .

j j j1 j Уравнение сохранения импульса для i - ой фракции твердой фазы:

N DVi s ii i iibi mif , i 1,..., N (12) m ij Dt j Уравнение сохранения импульса для жидкостой фазы:

N DVf f 1 f 1 bf (13) f f m i Dt iE i piI i, f pf I E - тензоры напряжений для жидкой и твердой f фаз; pi и pf - фазовые давления частиц i-ой фракции и жидкости, соответственно, которые определяются через поровое давление p i E pf 1 p, pi i p e ; i и E - вязкие тензоры напряжений f * 0, частиц i-ой фракции и жидкости; e - k 0 /* 1, * эффективное напряжение твердой фазы; bi, bf - удельные массовые силы, действующие на частицы фракции i и жидкость, bi gk r 2 Vi ;

i i di d s bf gk r 2 Vf ; mij e j j 3 j Vi Vj Vi Vj - 2 idi3 d j j сила взаимодействия между частицами фракций i и j; mif i Ui pi - сила взаимодействия между частицей и жидкостью, где i - коэффициент сопротивления связанный с передачей импульса между жидкостью и частицей i-ой фракции; Ui Vi Vf - скорость частиц относительно жидкости; i - объемная доля частиц i-ой фракции; i - плотность вещества частицы i-ой фракции; j - единичный вектор, направленный вдоль оси вращения; r - радиус вектор, напрвленный вдоль радиуса центрифуги; k - единичный вектор, направленный в абсолютной системе координат вертикально вверх.

Анализ размерностей системы (11)Ц(13) показывает, что решение задачи об оседании полидисперсной суспензии в тарельчатой центрифуге определяется следующими безразмерными параметрами: числом Фруда, Fr U0 gL0 - отношение инерционных сил к силам гравитации; числом Россби, Ro U0 0L0 - отношение сил инерции к силам Кориолиса;

седиментационным числом Рейнольдса, Re Ud 0 - отношение сил инерции к вязким силам; числом Фруда системы, F 0 L0 g - отношение центробежной к гравитационной составляющей массовых сил. В качестве характерных масштабных величин выбраны: L0 - характерное расстояние до оси вращения; d - диаметр частицы самой крупной фракции; U0 - скорость оседания одиночной частицы самой крупной фракции в неограниченной неподвижной жидкости; гидростатическое давление столба жидкости высотой L0 и плотностью 0, p0 L00g ; время оседания самой крупной частицей, t0 L0 U0. Здесь 0 - кинематическая вязкость жидкости; g - ускорение свободного падения; 0 - плотность жидкости; 0 - характерная скорость вращения центрифуги; * - объемная доля частиц, при которой они начинают соприкасаться друг с другом.

Оценка безразмерных параметров для водно-кварцевой суспензии с размером частиц не превышающим 25 мкм, помещенной в тарельчатую центрифугу радиусом 0.1 м, которая вращается со скоростью 0=100 сЦ1, дает следующие значения: Fr ~ 10Ц8, Fr/Re ~ 10Ц5, Ro ~ 10Ц5, Fr/Ro ~ 10Ц2. Тогда система уравнений (11)Ц(13) существенно упрощается и принимает вид:

i M fi 2r gk i V ai ,, i 1,..., N (14) t V 0 (15) N p e() (gk 2r)(1) E 1 (16) f j j f j Здесь N N N fi M V i i i k k .

j j j j k j1 k 1 j N 1 ai , V () ie i j j j .

i N e i j j i j 2 Здесь d1 (18f ), i di2 d1, .

j j f 2.5* 1 *, 0 * V .

* 0, Для закрытой тарельчатой центрифуги V 0 уравнения (14) становятся несвязанными с (15) и (16). В случае выполненения условий 2L0 g 1, const, e 0 уравнения (14) принимают вид:

j i iVi 0, i 1,..., N (17) t N 2rd1 где Vi V 1 k .

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

Система уравнений, описывающая изменение объемной доли частиц бидисперсной суспензии (d1 > d2) в соответствие с (17) имеет вид:

r1 r1V1 r2 r2V 0, 0, rmin r rmax (18) t r t r Здесь Vi - скорость движения частиц i-ой фракции относительно стенок центрифуги:

V1 U1(11) U22, V2 U11 U2(12) (19) Ui - скорость движения одиночных частиц i-ой фракции относительно суспензии:

1/U1 USt,1V 1, U2 USt,2V 1 A1 d1 d2 , 1 2.

1 2rdi2 USt,i - стоксовская скорость оседания одиночной частицы в 18f неподвижной жидкости под действием центробежных сил.

В начальный момент времени задается постоянное значение объемной доли частиц: i (0,r) i,0. На левой и правой границах задается равенство нулю потоков частиц: i (t,rmin)Vi(t,rmin ) 0 и i(t,rmax )Vi(t,rmax ) 0, i = 1, 2.

Система уравнений (18) с начальными и граничными условиями решается численно методом конечных объемов с определением потоков частиц на гранях ячеек по методу А.Н.Крайко (распад разрыва в среде, лишенной собственного давления).

Исследование решения системы (18) проводится в безразмерных t r Vi 2 переменных: , , i , где t* rmax V*, V* 0 rmaxd1 18f.

t* rmax V* В задаче присутствуют четыре параметра: 1,0, 2,0, D rmin rmax, d2 d1, значения которых задавались следующими: 0.1, D 0.1. Начальная объемная доля частиц 1,0 и 2,0 варьировалась.

На рис.5 проиллюстрирована динамика полей объемной доли мелкой и крупной фракций частиц для слабоконцентрированной суспензии 1,0 0.005, 2,0 0.045. Крупные частицы оседают достаточно быстро, и уровень их объемной доли в зоне седиментации непрерывно падает со временем, тогда как ее распределение по радиусу в зоне седиментации остается равномерным.

а) б) Рис.5. Распределения концентрации крупных частиц (а) и мелких частиц (б) вдоль радиуса центрифуги в различные моменты времени. 0=0.05, 1,0=0.005, 2,0=0.045.

=4.8810-4 (1); 0.477 (2); 0.942 (3); 1.40 (4); 2.75 (5); 34.47 (6); 77.58 (7); 162.04 (8).

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

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

Динамика профилей скоростей седиментации обеих фракций U1, U2 для данной суспензии показана на рис.6. Скорость крупных частиц, зависит только от суммарной объемной доли твердой фазы, линейно растет с ростом , приближаясь со временем к значениям, даваемым стоксовой формулой, рис.6 а).

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

а) б) Рис.6.Профили скорости седиментации крупных (а) и мелких (б) частиц в различные моменты времени. 0=0.05, 1,0=0.005, 2,0=0.045.

=4.8810-4 (1); 0.477 (2); 0.942 (3); 1.40 (4); 1.85 (5); 2.30 (6); 2.75 (7);

Для сильно-концентрированной суспензии при значениях 0=0.2, 1,0=0.18 и 2,0=0.02, картина седиментации подобна той, которая получена для слабоконцентрированной суспензии, с тем лишь различием, что двухволновая структура фронтов осаждения здесь выражена ярче.

На основе созданной модели оседания бидисперсной суспензии делается проверка правомерности использования формулы для экспериментального определения скорости оседания частиц в тарельчатой центрифуге, которая применялась в работе [Gerhart Ch., Dck J., Neee Th. Grundlagenuntersuchungen zur behinderten Sedimentation polydisperser Suspensionen bei der Hydrostromklassierung. Teil I: Untersuchungen in einer Laborzentrifuge // Aufbereitungstechnik. - 1999. - Vol. 40, №7. - P. 328Ц334.] и является точной, вообще говоря, для случая монодисперсной суспензии: Vi 0.5rp lni t, здесь rp - расстояние от центра до точки забора образца суспензии. В эксперименте удается взять 5-6 образцов суспензии с интервалом 5 сек, до того как частицы полностью осядут. Поэтому скорость седиментации частиц предлагается определять по конечно-разностной формуле:

ln i,k i,k rp (20) Vi,k 0.5 2 tk 1 tk Здесь k - порядковый номер взятого образца суспензии.

Будем называть лэкспериментально измереннными скоростями оседания частиц Ve скорости, полученные по формуле (20), в которой концетрации берутся из решения системы уравнений (18) в моменты времени tk и tk 1, а теоретически определенными скоростями Vt - скорости, найденные из соотношений (19) на момент времени tk 0.5 tk 1 tk 2.

На рис.7 показаны отношения значений скоростей Ve VSt в зависимости от отношения Vt VSt, полученные как для мелких частиц (темные символы), так и для крупных (светлые символы) в широком диапазоне объемной доли.

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

Рис.7. Сравнение теоретических и лэкспериментальных значений скоростей седиментации 1, 2 - 0=0.2, 1,0=0.18, 2,0=0.02;

3 и 4 - 0=0.05, 1,0=0.005, 2,0=0.045;

5 и 6 - 0=0.2, 1,0=0.02, 2,0=0.18;

7 и 8 - 0=0.05, 1,0=0.045, 2,0=0.005.

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

ri riVi 0, i 1,..., N. rmin r rmax (21) t r Здесь N 2rdm Vi V 1 fE i,i k2 fE k,k ; (22) i k 18f k N 1/fE i,i G mj - функция увлечения; - относительная j j ji i j объемная доля частиц j-ой фракции; mj ; - отношение размера j крупной частицы к размеру мелкой, которую способна увлечь первая;

i di dm ; dm - размер самой крупной частицы;

3 G 2.52 3 exp 5, - константа, теоретическое значение которой, подтверждённое экспериментами [ Dueck J., Neesse Th., Minkov L., Kilimnik D., Hararah M. Theoretical and experimental investigation of disturbed settling in a polydisperse suspension // Proceedings of ICMF-2004. Fifth International Conference on Multiphase Flow, Yokohama (Japan), 30 May - June, 2004. Y. Matsumoto, K. Hishida, A. Tomiyama, K. Mishima, S. Hosokawa (editors). Paper No. 106. P.1-8.], равно 10.

Начальные и граничные условия для системы уравнений (21) задаются аналогично предыдущей задаче: i 0,r 0mi,0, iVi t,rmin 0, iVi t,rmax 0.

В качестве масштаба длины выбирается внешний радиус центрифуги rout ; масштаб скорости - Стоксовская скорость оседания самой крупной частицы диаметром dm на расстоянии rmax от оси центрифуги, VSt,m rmax2dm 18f ; масштаб времени - tm rmax VSt,m. Расчеты проводились при следующих значениях исходных данных: rmin 0.008 м;

rmax 0.08 м; 0 0.2 ; dm 56 мкм; tm = 0.85 сек. Функция распределения частиц по размерам mi0 задавалась в соответствие с данными работы [7], [Hararah M.A. Settling of fine particle in dense polydisperse suspensions.

Doctoral dissertation. Erlangen, 2004. 100 p.] рис.8. Число фракций принималось равным 67.

Для сравнения результатов расчета с экспериментальными данными [0] используется число Фруда системы z rmax2 g (отно-шение центробежного ускорения на расстоянии rout от оси вращения к ускорению свободного падения).

На рис. 9 показана зависимость отношения суммарной объемной доли частиц к ее начальному значению в точке наблюдения rp 0.045 м от времени. Видно, что неучет эффекта увлечения мелких частиц крупными ( , толстая сплошная кривая) дает завышенные значения 0 по сравнению с экспериментальными данными (точки), тогда как решение системы уравнений (21) с использованием модели увлечения частиц ( 10, тонкая сплошная кривая) с хорошей точностью согласуется с экспериментом.

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

Рис.8. Исходная функция распределения Рис. 9. Изменение полной суммарной частиц по размерам. объемной доли частиц.

Сравнение рассчитанных относительных объемных долей i i,0 в зависимости от времени для фракций с частицами размером 2 и 3.2 мкм с экспериментальными данными показало, что экспериментальные значения могут превышать рассчитанные значения до 20%. Для этих фракций следует ожидать и наибольшего расхождения измеренных скоростей оседания с теоретическими значениями, рис.10 а).

i/i,i/i, а) б) Рис.10. Сравнение вычисленной относительнойобъемной доли для отдельных фракций мелких и крупных частиц с экспериментальными данными Для фракций с частицами 5 мкм и выше рассчитанные кривые изменения концентраций и качественно, и количественно согласуются с экспериментально измеренными, рис. 10 б).

Численное исследование седиментации частиц полидисперсной суспензии в тарельчатой центрифуге позволил через эволюцию профилей относительной объемной доли частиц различных фракций выявить следующие особенности. 1) Фронт осветления мелких частиц (d<5 мкм) практически не изменяет своего положения за время t tm 60, соответствующего времени проведения эксперимента, рис.11 а); 2) Для мелких и средних фракций в окрестности фронта осветления образуется область сгущения, вызванная противоположно направленным движением частиц, увлекаемых вытесняемой жидкостью, за счет оседания частиц более крупных фракций, рис.11 а), б); 3) Для крупных фракций область сгущения не возникает, рис.11 в); 4) В середине области (в экспериментах rp=0.045 м) радиальные градиенты объемной доли практически равны нулю до тех пор, пока фронт осветления не дойдет до точки наблюдения.

Как следует из системы уравнений (21) скорость седиментации частиц i/i,0 i/i,i/i, а) б) в) Рис.11. Распределение концентраций частиц вдоль r для различных моментов времени.

d=0.56 мкм (а); 10 мкм (б); 36 мкм (в) Для а), б): t/tm=5 (1); 10 (2); 20 (3); 40 (4); 60 (5), для в): t/tm =0.5 (1); 1 (2); 2(3); 3 (4).

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

N r ln(i ) r (1,...,n ) iVi j Vi , где (1,...,n ) , которая 2 t i j1 r r j вследствие выявленного свойства равенства нулю градиентов объемной доли частиц различных фракций упрощается до r ln(i ) Vi , (23) 2 t что позволяет сделать вывод о справедливости использования формулы (20), являющейся конечно-разностным аналогом формулы (23), для экспериментального определения скорости седиментации частиц полидисперсной суспензии в тарельчатой центрифуге.

При экспериментальном определении скорости седиментации частиц на основе измеренной их объемной доли, забор проб с суспензией осуществляется через определенный интервал времени t, величина которого может оказывать влияние на получаемые результаты. На рис. показаны относительные теоретические скорости седиментации частиц, вычисленнные с использованием формулы (22) (сплошные линии), и относительные лэкспериментальные скорости седиментации частиц, определенные с использованием формулы (20) для промежутков времени t tm 1 (светлые точки) и t tm 5 (темные точки).

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

Эволюция скорости седиментации для частиц с диаметром d=3.2 мкм V/VSt,m V/VSt,m а) б) Рис.12. Сравнение скоростей седиментации частиц.

(а) - d = 0.56 мкм; (б) - d = 3.2 мкм 1 - теоретические значения; 2, 3 - лэкспериментальные значения;

2 - t/tm=1; 3 - t/tm=5.

дана на рис.12 б). Здесь видна борьба двух эффектов. Вначале из-за оседания частиц (особенно крупных фракций) вытесняется жидкость и увлекает частицы этой фракции в сторону внутренней стенки, так, что скорости движения частиц становятся отрицательными, а потом, когда концентрация твердой фазы становится достаточно малой, процесс ускорения мелких частиц крупными становится более эффективным по сравнению с увлечением за счет вытеснения воды. На конечной стадии седиментации скорость устанавливается, сравниваясь со стоксовским значением. Очевидно, что расхождение между значениями скорости вычисленной по теоретической и по лизмерительной методике, связаны в основном с численным дифференцированием уравнения (23). Дополнительная ошибка может возникать из-за небольших градиентов концентрационных профилей различных фракций.

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

Определение скорости седиментации частиц суспензии (3), залитой через воронку в предварительно раскрученный корпус тарельчатой центрифуги (2), основано, согласно формуле (20), Рис.13. Схема забора суспензии из тарельчатой центрифуги на измерении в разные моменты времени объемной доли дисперсной фазы суспензии. Забор суспензии из корпуса центрифуги осуществляется с помощью шприца (6) через капилляр (4) и патрубок (5). Число капилляров, расположенных через равные углы, равно шести. Внешний край капилляра находится на достаточном удалении от внешнй стенки корпуса центрифуги, чтобы не оказаться внутри слоя седимента (1).

Объемная доля всей твердой фракции определяется путем взвешивания определенного объема суспензии по формуле (msus Vsus liq) (p liq ).

Объем суспензии измеряется при помощи пикнометра с точностью 0.0005 мл, взвешивание производится на аналитических весах с точностью 0.000005 г. Для получения распределения частиц по размерам используется прибор Mastersizer 2000 фирмы Malvern. Статистическая ошибка, связанная с ограниченным количеством измерений, составила 5%.

Рассматриваются следующие возможные систематические ошибки:

1) Ошибка, связанная со временем заполнения суспензией корпуса 1 центрифуги через воронку: tз t1 t2, t1 Vsus (g1 2tg1 6) - время истечения суспензии через воронку, t2=Ц1arch(R/r0) - время заполнения вращающегося корпуса центрифуги. Для значений Vsus = 150см3, g = 10 м/с2, =/4, R = см, r0 = 1 см, = 78.5 сЦ1 время заполнения составляет 0.1 с, что намного меньше промежутка времени (5 с) через который берется первая проба.

2) Ошибка в определении объемной доли частиц, обусловленная длительностью процесса забора пробы. Относительная ошибка в определении объемной доли составляет d22 (p liq)t (18liq). Для значений р = 2500 кг/м3, liq = 1000 кг/м3, liq = 10Ц3 Пас, = 78.5 сЦ1, t=0.с ошибка для частиц размером 1 мкм составляет 0.005%, а для 20 мкм - 2.2%.

3) Ошибка в определении объемной доли, обусловленная наличинм следа чистой жидкости за внешним краем капилляра, объем которой должен быть намного меньше объема забираемой суспензии. 2r2 Vs2 3. Для объема,суспензии, втягиваемой через один капилляр, Vs,1= 0.4 cм3 и для радиуса капилляра r = 0.35 мм ошибка, вносимая следом, составляет = 0.4 %.

4) Ошибка, обусловленная инерционностью частиц, при втягивании суспензии в капилляр. На основе анализа характерных времен релаксации частиц и времени забора проб получается оценка для минимально допустимого объема втягиваемой через капилляр суспензии Vs,1 Vlim d22rp(p liq)r2t (18liq). Для частиц размером 1 мкм предельно допустимый объем составляет 1.410Ц5 см3, а для частиц размером 20 мкм эта величина равна 5.510Ц3 см3, что намного меньше Vs,1= 0.4 cм3.

5) Ошибка, обусловленная внесением возмущения в положение седиментационного фронта из-за взятия проб, составляет 6Vs,1 2hrp.

При высоте центрифуги h = 1 см, 6V1 =2.4 см3, rр = 7 см величина = 0.8%.

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

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

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

В четвертой главе рассматривается оседание частиц полидисперсной суспензии с учетом увлечения мелких частиц крупными под действием массовых сил в классификаторе упрощенной формы (partition model) и исследуется влияние дополнительной инжекции жидкости на классификаионные характеристики аппарата.

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

Рис.14. Схема упрощенного классификатора Предполагается, что двухфазная смесь втекает в классификационный аппарат (рис.14) слева с постоянной скоростью U, а вытекает справа через верхнее и нижнее выходные отверстия. При этом каждая частица j-й фракции смещается вниз (седиментирует) с некоторой скоростью VS, j, зависящей от размера частиц d этой фракции. Коэффициент турбулентной диффузии j частиц D считается константой. Ширина верхнего отверстия hov, а ширина всего аппарата h.

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

jj Pe 0 (25) j j На входе в классификатор задаются начальные концентрации частиц 1. (26) j На стенках классификатора ставятся условия равенства нулю потоков частиц j Pe 0. (27) j j где 0, j, x x, x h2U D, y h ; Ре hVS, j D ; j 1, N j j j Безразмерный параметр Реj характеризует размер частиц j-ой фракции.

Сепарационная функция, определяющая долю частиц данной фракции, отводимую через нижнее отверстие, определяется как (28) T Pe, , d jj S S где S = hov / (h - hov) - сплит-параметр.

Важнейшими параметрами, характеризующими классификационный процесс, являются: а) параметр Ре[50], дающий диаметр разделения частиц, поступающих на 50% в нижнее отверстие; б) острота разделения Q, которая определяется отношением диаметров частиц, соответствующих значениям сепарационной функции 0.25 и 0.75, Q d25 d75.

Точное решение задачи (25)Ц(27) имеет вид:

Pe j , exp Pe jj exp Pe j Pe n (29) 2Pe exp j 1 exp j j 1 Pe ncos n Pe 2 n sin 2 n j 2 j exp n2 Pe2 4 n Pe 4 n j Для достаточно больших , ограничиваясь первым членом суммы в (29), получается автомодельное решение , , 1 , exp Pe2 4 (30) j j j j где Pe j , exp Pe .

jj exp Pe j Сравнение точного и автомодельного решений, показывает, что автомодельное решение, достаточно близко к точному при > 0.1.

Выражение для характерной длины релаксации частицы j-ой фракции можно найти, положив Pe2 4 1, откуда j r, j Pe2 4 (31) j VS2, j U D Эта формула в размерных переменных имеет вид: . Здесь xr, j h2 4D выделяются следующие характерные времена процесса классификации, а именно: диффузионное время (td = h2/D), время седиментации (ts,j = h/Vs,j) и время пребывания частиц в зоне установления их концентрации (tr,j = xr,j /U ).

Из (31) видно, что длина релаксации уменьшается при увеличении параметра Реj. Поскольку большие значения Реj соответствуют более крупным частицам, концентрация таких частиц устанавливается на более коротком расстоянии от входа в аппарат, чем в случае мелких частиц.

При численном решении задачи (25)Ц(27) получены зависимости диаметра разделения Ре[50] (рис.15) и остроты разделения от (рис.16).

Рис. 15. Влияние значения сплит- Рис. 16. Зависимость остроты разделения параметра S на зависимость Ре[50] от от длины аппарата. Сплошные кривые - длины аппарата S = 4 (1); 9 (2); 19 (3 ) численное решение; штриховые - автомодельное решение S = 4 (1); 9 (2) Из рис.15 видно, что на начальном участке кривой значение диаметра разделения падает Ре[50] при увеличении , а при дальнейшем увеличении значение Ре[50] приходит к постоянной величине, поскольку оседание более мелких частиц уравновешивается диффузией. При увеличении сплитпараметра S, когда размер нижнего слива уменьшается относительно размера верхнего, растет и значение Ре[50].

Вид зависимости остроты разделения от длины аппарата , при различных значениях параметра S, полученной численным способом и на основе автомодельного решения (30), показан на рис.16. Расхождения при малых значениях естественны, поскольку автомодельное решение верно только для достаточно больших значений . При увеличении значение остроты разделения Q падает, за исключением случая коротких аппаратов.

Вид функции Q() определяется характером установления Pe[j25] и Pe[j75].

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

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

Во втором параграфе проводится численное исследование fish-hook эффекта на основе модели упрощенного классификатора, рис.14.

Предполагается, что распределение частиц по размеру подчиняется закону RRSB (РозенЦРаммлерЦСпенлингЦБеннет), а мелкие частицы увлекаются крупными по механизму, описанному в первой главе.

Система уравнений, описывающая изменение относительной объемной доли седиментирующей совокупности полидисперсных частиц в классификационном аппарате, начальные и граничные условия имеют вид аналогичный (25)Ц(27), в которых Pe PemWS, j.

j Здесь WS, j VS, j VSt,m - безразмерная скорость седиментации частицы jой фракции, которая для распределения RRSB q() mexp m m записывается в виде:

m WS, j , H g fE g fE mt tm exp(t)dt j j j 1 m Pem hVSt,m D ; j d d* ; fE (x) exp(t)dt ;

j t x VSt,m - стоксовская скорость оседания частицы размером d* ; d* - средний размер частиц в суспензии; параметр m характеризует остроту распределения частиц по размерам.

.

3 H 1 1 0.6 ; g 2.252 3 exp 5 ;

n , 0, j ,.

j j Расчеты проводились при следующих значениях параметров:

10, m 1.7. Общее число фракций n принималось равным 40, безразмерный диаметр частиц j-ой фракции определялся по правилу j 0.1100K 1, где j номер фракции.

j Параметрами задачи являются: объемная доля частиц, подаваемых в аппарат 0, средний размер частиц суспензии Pem, геометрический сплит- параметр S.

На рис. 17 показано поведение сепарационных кривых, полученных на выходе из классификационного аппарата (=1) при разных значениях начальной объемной доли твердой фазы суспензии для Pem 10, S 9.

Видно, что для частиц с размером больше, чем 0.5, T 3 сепарационная функция 0. монотонно возрастает. Доля 0 0.частиц с размерами, 0.отводимая через нижнее сливное T отверстие, монотонно 0.4 уменьшается с увеличением, что объясняется выносом мелких частиц через нижнее сливное 0.отверстие за счет их увлечения более крупными частицами.

Увеличение начальной 0.1 1 объеминой доли твердой фазы в суспензии ведет к монотонному Рис.17. Влияние объемной доли твердой фазы на сепарационную функцию.

снижению значений =1. =0.02 (1); 0.04 (2); 0.08 (3); 0.16 (4) сепарационной функции в области 0.5<10, что вызвано замедлением скорости седиментации частиц в силу стесненных условий оседания, ведущего к уменьшению относительной концентрации крупных частиц в нижнем выходном отверстии. При этом зерно разделения (значение , при T=0.5) смещается в сторону более крупных частиц - факт хорошо известный из экспериментов и практики.

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

Для характеристики немонотонности сепарационной кривой вводится понятие глубины fish-hook эффекта как разница h = T(0.1) Ц1/(1+S). На рис.18 а) показана зависимость h от начальной объемной доли твердой фазы при разных значениях числа Pem, которая является немонотонной с ярко выраженным максимумом. Такое поведение глубины эффекта объясняется двумя противоборствующими факторами - увлечением мелких частиц крупными с одной стороны, и замедлением процесса оседания вследствие возрастания вязкости и плотности суспензии - с другой. Из рисунка также видно, что увеличение числа Pem, которое соответствует или уменьшению коэффициента диффузии частиц, или увеличению их собственной скорости седиментации, приводит к возрастанию максимума h и смещению его в область низких значений начальной объемной доли. Здесь видно влияние функции распределения частиц по размерам на характеристики сепарации.

а) б) Рис.18. Зависимость глубины увлечения от суммарной объемной доли частиц на входе.

а - расчет; б - эксперимент.

Зависимость h от начальной объемной доли твердой фазы при разных значениях числа Pem, полученная в расчетах, качественно полностью соответствует измерениям, приведенным в [Gerhart Ch. Untersuchungen zum Trennverhalten in Hydrozyclonen niedriger Trennkorngren: Dissertation.

Erlangen: Universitt Erlangen-Nrnberg, 2001.], рис.18 б) В одной из серий экспериментов измерялась глубина fish-hook эффекта в зависимости от концентрации подаваемой в гидроциклон суспензии. Для двух используемых аппаратов с диаметрами рабочего цилиндра 25 и 39 мм устойчиво наблюдался максимум h при =0.04. Значение максимума уменьшалось в 1.9 раз при смене большего гидроциклона на меньший, т.е. при изменении Pem в 1.56 раз. Хотя из рис.18 а) следует гораздо более сильная зависимость от Pem, следует учесть, что изменение размера гидроциклона при удержании давления приводит к сильному изменению поля течения, что не отражено в проведенных расчетах.

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

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

x y hov h Vin L H Рис.19. Схема инжектирования классификатора Система уравнений, описывающая изменение относительной объемной доли j-ой фракции частиц, в безразмерной форме имеет вид:

U( ) jj WS, j W() Pem 0. (32) j слив Верхний Вход слив Нижний Условие на входе классификатор: 1.

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

j j PemWS, j 0 ; Pem WS, j W () 0. (33) j j 0 Поле скоростей описывается функциями:

1, 0 L Hin U 1 PemWin L Hin, L Hin L 0, 0 L Hin W Win, L Hin L Решение задачи (32)Ц(33) определяется следующими параметрами:

Pem hVSt,m D - параметр Пекле (характеризуется высотой классификатора, скоростью Стокса для масштабной частицы и коэффициентом турбулентной диффузии); Win Vin VSt,m - отношение скорости инжекции к скорости Стокса для масштабной частицы; Hin H x* - безразмерная ширина отверстия для дополнительной инжекции жидкости; L L x* - безразмерная длина классификатора; параметр m, отвечающий за остроту распределения частиц по размерам исходной суспензии.

Увеличение скорости инжекции Win ведет к снижению содержания мелких частиц в потоке крупных, выгружаемых через нижний слив классификатора T(0.01) и, Win0=следовательно, качество Win0=разделения повышается (Рис.20).

0.Win0=Кроме того, увеличение скорости Win0=1инжектирования ведет к 0.снижению глубины fish-hook эффекта.

T Расчеты показывают, что 0.возрастание скорости инжекции ведет к увеличению размера зерна 0.разделения, что является нежелательным для эффективного классификационного про0.01 0.1 1 10 1цесса. Размер зерна разделения j зависит от концентрации частиц и в случае более разжиженных Рис.20. Влияние начальной скорости суспензий он уменьшается.

инжекции на поведение сепарационной Увеличение плотности суспензии функции.

=0.04; Hin=0.001; Pem=10; m=1.7. ведет к повышению остроты разделения.

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

для сепарационной функции T ( ) ;

j 1 S exp WS, j ( ) 0.5Win Pem j для зерна разделения Pem 50 050 1 Win ;

2ln S для остроты разделения PemWin 1 PemWin Q Q0 1 1 2 ln S 3 2 ln 3S Здесь 050, Q0 - зерно разделения и острота разделения при отсутствии [ инжекции, соответственно: Q0 ln S 3 ln 3S, WS, j (050]) ln(S) Pem Для малых Win острота разделения растет с увеличением скорости инжекции, а при Win острота разделения стремится к единице.

а) б) Рис.21. Сравнение расчетных и экспериментальных зависимостей.

а - T0(qinj); б - d[50](qinj).

1 - экперимент; 2 - расчет по (32)Ц(33) Для сравнения модели (32)Ц(33) с экспериментальными данными [Mohamed Galal Farghaly Aly. Controlled Wash Water Injection to the Hydrocyclone Underflow. - Erlangen - 2009. - 155 p.] были проведены расчеты со следующими исходными данными: n=51, 0=0.094, S=9, Pem = 10, что соответствовало dm = 21 мкм, h=50 мм, H=10Ц3 мм.

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

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

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

уравнение сохранения массы смеси:

mUm 0, (1) уравнение сохранения массы для k-ой фракции частиц с учетом их диффузии за счет турбулентности:

k m k UM,k Um Dturb,k Dturb,k 0, k=1, 2,Е, N (2) s k m уравнение сохранения количества движения смеси:

N mUmUm p m turb mg kkUM,kUM,k . (3) k 0 Здесь m - плотность смеси, s - плотность вещества частицы, Um - вектор скорости смеси, k - объемная доля k-ой фракции дисперсной фазы; g - ускорение свободного падения. Диффузионная скорость k-ой фракции UM,k выражается через скорость скольжения k-ой фракции относительно s N жидкой фазы Uk,liq из соотношения: UM,k Uk,liq U, а скорость m j1 j j,liq скольжения k-ой фракции дисперсной фазы находится из зависимости:

k s m Uk,liq a, в котороой ускорение частицы a согласно модели fdrag,k s смеси можно найти через градиент давления: a g p. Тензор вязких m напряжений смеси определяется через вязкость суспензии 1.m liq 1 * и тензор скоростей деформации суспензии:

m m Um UT, а тензор турбулентных напряжений смеси m определяется через тензор Рейнольдсовых напряжений: turb mR. Для коэффициента сопротивления частицы используется формула ШиллераНауманна: fdrag, k 1 0.15Re0.687, пригодная для Rek 1000, где k liqdk Rek Uk,liq,* 0.62. Согласно модели, изложенной в первой главе, liq мелкие частицы, увлекаемые крупными, приобретают скорость пропорциональную квадрату кажущегося диаметра частицы, который в случае полидисперсной суспензии будет dk G() fE (dk ), где 1/ N 3 G() 2.51/3 exp 5, fE (dk ) d - функция увлечения j j j d dk j частицы. Поэтому время релаксации частицы в формуле для скорости скольжения частиц относительно жидкости имеет следующий вид:

s k dk G() fE (dk ).

18liq Для определения параметров поля турбулентности привлекается модель Рейнольдсовых напряжений, описываемая следующей системой уравнений:

Уравнение переноса для рейнольдсовых напряжений:

turb mUmR m 1 R P , (4) где P m R Um UT RT тензор, отвечающий за производство - m турбулентных напряжений; - тензор, отвечающий за перераспределение турбулентности за счет быстрых деформаций, медленных деформаций и за счет перераспределения напряжения вблизи стенки; - тензор скорости диссипации, Im. Коэффициент турбулентной вязкости вычисляется по kформуле turb C, в которой кинетическая энергия турбулентных пульсаций определяется через тензор рейнольдсовых напряжений:

k tr(R).

Уравнение переноса для скорости диссипации турбулентной энергии:

turb 1 mUm m C1 2 tr(P) k C2m k, (5) где 1.0 ;1 0.82 ; C1 1.44 ; C2 1.92 ; C 0.09 - параметры модели. Коэффициент турбулентной диффузии частиц Dturb,k определяется sdk по формуле Дика-Неессе mDturb,k turb 11.6 .

k liq В соответствие с допущением об осевой симметрии закрученного Rxx Rxr Rx течения Um Um, Vm,Wm, R Rxr Rrr Rr , i j.

x r Rx Rr R Граничные условия для системы уравнений (1)Ц(5) имеют следующий вид:

на входе в гидроциклон:

Rxx Rrr R 0.1Um, Rxr Rx Rr 0.

k3 C 4, 0.07 Dhyd 2ainbin Dhyd - гидравлический диаметр входного сечения,.

ain bin QF QF Um 0, Vm , Wm Dcain ainbin на нижнем и верхнем выходах из гидроциклона задавалось давление равное атмосферному p pвнеш, для остальных параметров задавались условия x 0, Um, Vm, Wm, Rxx, Rrr, R, Rxr, Rx, Rr, на оси симметрии задавались условия симметрии Vm 0, Wm 0, Rxr 0, Rx 0, Rr 0, 0, Um, p, Rxx, Rrr, R,.

r на стенке гидроциклона задавалось условие равенства нулю скорости смеси 0, Um, Vm, Wm. В системе координат (,,), связанной со стенкой, где - продольная координата, - нормальная координата, - бинормальная координата напряжения Рейнольдса в пристенной области определялись через кинетическую энергию турбулентых пульсаций k.

R 1.098k, R 0.247k, R 0.655k, R 0.255k Скорость диссипации турбулентной энергии в пристенной области определялась на основе гипотезы локального равновесия между 3/C k3/производством турбулентной энергии и ее диссипацией , y здесь - постоянная Кармана, 0.4187, y - расстояние до стенки.

Решение системы уравнений (1)-(5) выполнялось с помощью противопоточной схемы второго порядка точности с привлечением алгоритма PRESTO для расчета давления на гранях ячеек. Согласование между полем давления и полем скорости реализовано на основе алгоритма PISO.

Вся совокупность частиц разбивалась на 12 фракций, объемная которых приведена в Таблице № 1 2 3 4 5 d, мкм 0.66 0.88 1.63 2.73 4.43 6.k103 0.596 0.743 0.972 1.236 1.531 1.6№ 7 8 9 10 11 d, мкм 9.44 12.71 17.74 24.91 36.90 45.k103 2.293 3.388 3.301 2.312 0.617 0.1 На входе в гидроциклон скорости частиц и жидкости совпадают.

Объемный расход суспензии на входе в гидроциклон равнялся 74.9 л/мин.

Плотность вещества частиц s = 2650 кг/м3. Объемная доля частиц составила 0.0189, что соответствует содержанию твердой фазы в суспензии 50 г/литр.

Для анализа полученных результатов все фракции частиц были разбиты на три основные группы: мелкие - с размером менее 5 мкм, среднеразмерные - с размером от 5 мкм до 15 мкм, крупные с размером от 15 мкм и выше. Объемная доля мелких частиц, поступающих в аппарат, составила 0.0051, средних - 0.0074, крупных - 0.0064. На рис.представлены распределения объемной доли мелких частиц поперек гидроциклона в различных его сечениях: а - выходная граница верхнего слива; б - середина цилиндрической части гидроциклона, х=0.165 м; в - выходная граница нижнего слива. Учет эффекта увлечения мелких частиц крупными сильно изменяет ее распределение объемной доли для мелкоразмерных фракций, и практически не сказывается на распределении среднеразмерных и крупных частиц. Видно, что доля мелких частиц (кривая 1) снижается в верхнем сливе (рис. 22 а) и в самом аппарате (рис.22 б), и увеличивается в нижнем сливе, (рис.22 в), по сравнению со случаем без учета рассматриваемого эффекта. К тому же, происходит перераспределение содержания мелких частиц поперек гидроциклона. Доля мелких частиц возрастает по направлению к стенке гидроциклона (рис.22, б и рис. 22 в), что объясняется повышенной скоростью мелких частиц за счет их увлечения крупными.

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

а) б) в) Рис.22. Профили объемной доли мелких частиц.а - верхнее сопло; б - цилиндрическая часть; в - нижнее сопло. 1 - без учета увлечения, 2 - с учетом увлечения.

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

23 б,в).

а) б) в) Рис.23. Радиальная скорость частиц VM,k в цилиндрической части гидроциклона.

а - фракция №4; б - фракция №8; в - фракция №11.

1 - без учета увлечения, 2 - с учетом увлечения.

Если не учитывать эффекта увлечения мелких частиц крупными, то сепарационная кривая монотонно поднимается до 1 при увеличении размера частиц (рис. 24, кривая 1).

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

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

мелкоразмерных частиц (рис. 24, кривая 2). Для сравнения на рис.24 приведены экспериментальные данные из работы [17] списка публикаций по теме диссертации.

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

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

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

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

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

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

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

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

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

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

Проведен анализ воздействия инжекционного потока на характеристики классификационного аппарата.

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

СПИСОК ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ Публикации в изданиях, рекомендованых ВАК 1. Дик.И.Г., Миньков Л.Л., Неессе Т. Гидродинамическая модель ускорения седиментации мелких частиц в бидисперсной суспензии // Теплофизика и аэромеханика. - 2001. - Т.8, №2. - С.283Ц294.

2. Дик И.Г., Миньков Л.Л. Кинематика обтекания шара медленным потоком взвешанных частиц в приложении к седиментации бидисперсной суспензии // Инженерно - физический журнал. - 2002. - Т.75, №4. - С. 3Ц11.

3. Дик И.Г., Миньков Л.Л., Ларионова Н.В., Неессе Т. Седиментация бидисперсной суспензии в тарельчатой центрифуге // Теплофизика и аэромеханика. - 2002. - Т.9, №3. - С.481Ц494.

4. Дик И.Г., Килимник Д.Ю., Миньков Л.Л., Неессе Т. Измерение скорости седиментации мелкодисперсных частиц в тарельчатой центрифуге // Инженерно-физический журнал. - 2003. - Т.76,№4. С.7Ц 17.

5. Neesse Th., Dueck J., Minkov L. Separation of finest particles in hydrocyclones // Minerals Engineering. - 2004.ЦVol.17.Ц P.689Ц696.

6. Dueck J., Neesse Th., Minkov L. Zum Einfluss der Partikelrotation in dichten Suspensionen // AufbereitungsTechnik. - 2006. - Vol.47, No. 1-2. - P. 6Ц14.

7. Дик И.Г., Миньков Л.Л., Пикущак Е.В. О сепарационных кривых проточного классификационного аппарата конечной длины // Инженерно-физический журнал. - 2006. - Т.79, №3. - С.171Ц178.

8. Дик И.Г., Миньков Л.Л., Пикущак Е.В. Моделирование "Fish-hook" эффекта в классификационном аппарате // Инженерно-физический журнал. - 2007. - Т.80, №1. - С.60Ц69.

9. Дик И.Г., Миньков Л.Л., Пикущак Е.В. Влияние функции распределения частиц по размерам в полидисперсной суспензии на сепарационный процесс в классификационном аппарате // Вестник Томского государственного университета. Математика и механика. - 2007. - №1.ЦС.63Ц71.

10. Миньков Л.Л., Пикущак Е.В., Дик И.Г. Моделирование седиментации частиц полидисперсной суспензии в тарельчатой центрифуге // Теплофизика и аэромеханика. - 2009. - Т.16,№1. - С. 79Ц88.

11. Дик И.Г, Пикущак Е.В., Миньков Л.Л. Моделирование изменения характеристик разделения классификатора путем инжекции воды в аппарат // Теплофизика и аэромеханика. - 2009. - Т.16, №2. - С. 261Ц 273.

12. Dueck J., Pikushchak E., Minkov L., Farghaly M., Neesse Th.. Mechanism of hydrocyclone separation with water injection // Minerals Engineering. - 2010. - V.23. - P.289Ц294.

13. Hararah М.А., Endres E., Dueck J., Minkov L., Neesse Th. Flow conditions in the air core of the hydrocycloner // Minerals Engineering. - 2010. - Vol.23.Ц P. 295-300.

14. Крохина А. В., Дик И. Г., Неессе T., Миньков Л. Л., Павлихин Г. П.

Исследование гидродинамики гидроциклона с дополнительным двухструйным инжектором // ТОХТ. - 2011. - Т.45, №2. - С.227Ц235.

15. Миньков Л.Л., Крохина А.В., Дик И.Г. Гидродинамические механизмы влияния инжекции на классификационные характеристики гидроциклона // ИФЖ. - 2011. - Т.84, №4. - С.747Ц758.

16. Minkov L., Dueck J. CFD-modeling of a flow in a hydrocyclone with an additional water injector // Компьютерные исследования и моделирование.

Ц 2011. - Т. 3, № 1. - С. 63-76.

17. Миньков Л.Л., Крохина А.В., Дик И. Г.Расходные характеристики гидроциклона со встроенным инжектором // Теплофизика и аэромеханика - 2011. - Т.18, №3. - С.413Ц426.

18. Pikushchak E., Dueck J., Minkov L. Modeling of Cross Flow>

Year 7. - Issue. 78. - 2011. - P.141Ц143.

19. Дик И.Г., Миньков Л.Л. Нестоксовская седиментация в приложении к анализу взаимодействия частиц в суспензии //ИФЖЦ2012. - Т.85, №1. - С.18Ц26.

20. Dueck J., Krokhina A.V., Minkov L.L. Controlling Characteristics of Hydrocyclone via Additional Water Injection // Theoretical Foundations of Chemical Engineering. - 2012. - Vol. 46, No. 3. - P. 296Ц306.

Публикации в других научных изданиях 21. Дик И.Г., Миньков Л.Л. Гидродинамическая модель седиментации бидисперсной суспензии // Сборник трудов Всероссийской науч. конф.

"Байкальские чтения по математическому моделированию в синергетических систе-мах",10-23 Июля 1999 Улан-Удэ-Томск. - Изд-во ТГУ. - С.146Ц149.

22. Dueck J., Neee Th., Gerhart Ch., Minkov L. Behaviour of finest particles in hydrocyclones. // Science & Technology of Filtration and Separations for the 21st Century, Editors Shiao-Hung Chiang and Samuel E. Lee, American Filtration & Separation Society, Vol.15, Pittsburgh, 2001.

23. Dueck J., Neee Th., Gerhart Ch., Minkov L. Residence time model for the settling of a polydisperse suspension. // Proceedings of the ICMF-2001.

Advances in Filtration and Separation Technology. Proceedings of the ICMF2001 (4th International Conference on Multiphase flow, Editor E.

Michaelides, Tulane University), May 27 to June1, 2001, New Orleans LA, USA.

24. Дик И.Г., Миньков Л.Л., Неессе Т. Моделирование седиментации бидисперсной суспензии в тарельчатой центрифуге. // Труды Международной конференции "Байкальские чтения-II по моделированию процессов в синергетических системах" 18-23 июля 2002. Улан-Удэ, Томск. Изд-во Том. ун-та, 2002. - С.193Ц197.

25. Дик И.Г., Миньков Л.Л. Седиментация монодисперсной суспензии в тарельчатой центрифуге. // Труды Международной конференции "Байкальские чтения-II по моделированию процессов в синергетических системах" 18-23 июля 2002. Улан-Удэ, Томск. Изд-во Том. ун-та, 2002. - С. 188Ц192.

26. Пикущак Е.В., Дик И.Г., Миньков Л.Л. Модель fish-hook эффекта в классификационном аппарате // Доклады IX Всероссийской научнотехнической конференции Физика и химия высокоэнергетических систем - Томск 2003. - С. 39Ц40.

27. Dueck J., Neesse Th., Minkov L.,Kilimnik D., Hararah M. Theoretical and experimental investigation of disturbed settling in a polydisperse suspension.

// ICMF-2004. Fifth International Conference on Multiphase Flow.

Yokohama, Japan, May 30-June 4, 2004. Paper No.106, P.1Ц28. Миньков Л.Л., Дик И.Г. Влияние силы Саффмэна на процесс увлечения мелких частиц при оседании крупных частиц в суспензии // Фундаментальные и прикладные проблемы современной механики.

Доклады IV Всероссийской научной конференции. Томск, 5-7 октября 2004 г. Изд-во Том.ун-та. - 2004. - С.123Ц124.

29. Миньков Л.Л., Дик И.Г. Роль силы Саффмэна при совместном оседании мелких и крупных частиц в суспензии // Известия вузов.

Физика. - 2004. Ц№10. - С.81Ц88.

30. Дик И.Г., Миньков Л.Л., Пикущак Е.В. Исследование основных сепарационных характеристик в классификационном аппарате для упрощенной постановки задачи // Доклады XI Всероссийской научнотехнической конференции Физика и химия высокоэнергетических систем - Томск 2005 - С. 21Ц22.

31. Dueck J., Minkov L.L. Collective effects by settling of polydisperse dense suspension // Eurasian Physical Technical Journal. - 2005. - Vol.2, No.1(3).

Ц P.47Ц63.

32. Миньков Л.Л., Пикущак Е.В., Дик И.Г. Обоснование метода измерения скорости седиментации частиц отдельных фракций в плотной полидисперсной суспензии // Материалы V Всерос. науч. конф.

"Фундам. и прикл. проблемы современной механики". Томск 3-октября 2006 г. - 2006. - С. 394 - 396.

33. Hararah M.A., Dueck J., Neesse Th., Minkov L. Investigation of disturbed settling in polydisperse suspensions // Eurasian Physical Technical Journal. - 2006. - Vol.3, No.2(6). - P.27Ц37.

34. Пикущак Е.В., Дик И.Г., Миньков Л.Л. Упрощенная модель классификатора с дополнительным впрыском жидкости на выходе // Известия ВУЗов. Физика. Ц2008. - Т.51, №8/2. - С. 201Ц206.

35. Dueck J.G., Minkov L.L., Pikushchak E.V., Neesse Th. Origin of the "Fishhook" Effect // HEAT 2008, Proc. of the 5th Intern. Confer. on Transp.

Phenomena in Multiphase Systems, June 30-July 3, 2008, Bialystok, Poland.

V.2. ЦP.205Ц212.

36. Миньков Л.Л., Дик И.Г, Джалал М., Пикущак Е.В. Особенности течения в гидроциклоне со встроенным инжектором // Материалы VI Всерос. науч. конф. "Фундам. и прикл. проблемы современной механики". Томск 30 сентября - 2 октября 2008 г. - Томск; Изд. Том. гос.

ун-та. 2008 - С.381 - 382.

37. 28. Dueck J., Pikushchak E., Minkov L., M.Galal, Th.Neesse. Simulation of Water Injection in Hydrocyclones // Proceedings of Physical Separation '09, Falmouth, UK, June 16-17, 2009. - P. 1Ц13.

38. Dueck J., Krokhina A., Minkov L.L., Neesse T. Hydrodynamics of a cyclone with wash water injection // Proceedings of 7th World Conference on Experimental Heat Transfer, Fluid Mechanics and Thermodynamics 28 June - 03 July 2009, Krakow, Poland (ExHFT-7), P.1953Ц1960.

39. Пикущак Е.В., Дик И.Г., Миньков Л.Л.. Влияние дополнительной инжекции на классификационный процесс // Материалы VI Всерос.

науч. конф. "Современная баллистика и смежные вопросы механики".

Томск 17 - 19 ноября 2009. г. Томск; Изд. Том. гос. ун-та. Ц2009. - С.2Ц 286.

40. Minkov L., Pikushchak E., Galal M., Dueck J. On the improving the>

41. Dueck J., Krokhina A., Minkov L., Neesse Th. Experimental and numerical investigation of hydrodynamic characteristics of the hydrocyclone with water injection // Separation processes. Summaries 2. 19 th International Congress of Chemical and Process Engineering CHISA 2010. 28 August - 1 September 2010, Prague, Czech Republic. P3.289. - P.800 - 801.

42. Dueck J., Minkov L., Neesse Th.. Effect of disturbed settling on convective particle transport in a Hydrocyclone // The 21st International Symposium on Transport Phenomena. 2Ц5 November 2010, Kaohsiung City, Taiwan. - 2010. - P. 1 - 7.

43. Миньков Л.Л., Крохина А.В., Дик И.Г. Влияние инжектирования воды на характеристики гидроциклона // Труды Томск. ун-та. - Т.276. - Сер.физ-мат. - Томск: Изд-во Том. ун-та, 2010. - С.9 - 18.

44. Minkov L., Farghaly M., Dueck J. Hydrodynamics of Hydrocyclone with Built-in Water Injector //Eurasian Physical Technical Journal. - 2011.ЦV.8, No.1(15). - P.29Ц41.

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