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

Расчёты и пояснения

Оглавление · Предметный указатель

GWEB_C0097

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

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

Что это даёт на практике

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

Что делает вихрь

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

Привычная запись центробежной инерционной силы:

Fц = m·vθ²/r, а соответствующее ускорение aц = vθ²/r.

Здесь m — масса частицы в килограммах, v_θ — её местная окружная скорость в метрах в секунду, r — радиус траектории в метрах. Сила получается в ньютонах, ускорение — в м/с². Это описание во вращающейся системе отсчёта; физический смысл прост: чем быстрее поворот и меньше его радиус, тем труднее частице повторять изгиб газового потока.

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

Небольшой расчёт: сто ускорений свободного падения

Возьмём условную окружную скорость частицы 10 м/с и радиус 0,10 м:

aц = 10² / 0,10 = 1000 м/с² ≈ 102g, где g ≈ 9,81 м/с².

Для сферической частицы m = ρp·π·d³/6. Здесь ρ_p — плотность частицы в кг/м³, d — диаметр в метрах, π ≈ 3,1416; делитель 6 получается из геометрической формулы объёма шара через диаметр. Возьмём d = 10 мкм и ρ_p = 2000 кг/м³. Последовательность на обычном калькуляторе:

1. Переводим размер: 10 мкм = 10·10⁻⁶ м = 1·10⁻⁵ м.
2. Возводим в куб: d³ = 1·10⁻¹⁵ м³.
3. Считаем объём: πd³/6 ≈ 0,5236·10⁻¹⁵ м³.
4. Умножаем на плотность: m ≈ 2000·0,5236·10⁻¹⁵ = 1,047·10⁻¹² кг.
5. Умножаем на ускорение: Fц ≈ 1,047·10⁻¹²·1000 = 1,047·10⁻⁹ Н.
6. Для сравнения, вес: mg ≈ 1,047·10⁻¹²·9,81 = 1,027·10⁻¹¹ Н. Отношение сил ≈ 102.

Это учебные исходные данные, а не режим выбранного циклона. И сто g здесь ещё не «стократная очистка»: газ сопротивляется относительному движению частицы. Чтобы узнать, успеет ли она отделиться, нужно перейти от силы к скорости, а затем к пути и времени.

Скорость витания: почему мелочь так долго остаётся в газе

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

Для небольшой сферы при медленном обтекании действует закон Стокса. Из баланса веса с учётом выталкивания и сопротивления получается:

wос = (ρp − ρg)·g·d² / (18μ).

ρ_p — плотность частицы, ρ_g — плотность газа, обе в кг/м³; d — диаметр в метрах; μ — динамическая вязкость газа в Па·с. Насыпная плотность топлива сюда не подходит: она включает пустоты между частицами. Для пористого углеродистого остатка нужна плотность, согласованная с массой и внешним размером именно частицы.

Число 18 здесь не подгоночный коэффициент для золы. В стоксовом приближении сила сопротивления сферы Fс = 3πμdw, где F_с измеряется в ньютонах, w — относительная скорость в м/с. Коэффициент 3π следует из решения задачи медленного вязкого обтекания сферы. Приравнивая это сопротивление силе (ρp − ρg)·πd³g/6, сокращаем π и получаем 18 = 3·6 в знаменателе скорости. Для иной формы или иного режима менять «18 на глаз» вместо выбора подходящей модели нельзя.

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

Reₚ = ρg·w·d / μ, где w — скорость частицы относительно газа, м/с. Re_p безразмерно: единицы в числителе и знаменателе сокращаются. Индекс p означает частицу; это не число Рейнольдса всего аппарата.

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

Десять и двадцать микрометров

Зададим условные свойства газа: ρ_g = 0,50 кг/м³, μ = 3,0·10⁻⁵ Па·с. Пусть плотность сферических частиц равна 2000 кг/м³. Это выбранные числа для сравнения, а не табличные свойства любой золы при любой температуре.

Для десятимикронной частицы распишем вычисление:

1. Разность плотностей: 2000 − 0,50 = 1999,5 кг/м³.
2. Квадрат диаметра: (1·10⁻⁵)² = 1·10⁻¹⁰ м².
3. Числитель скорости: 1999,5·9,81·10⁻¹⁰ = 1,96151·10⁻⁶.
4. Знаменатель: 18·3,0·10⁻⁵ = 5,4·10⁻⁴.
5. Делим: wос = 1,96151·10⁻⁶ / (5,4·10⁻⁴) ≈ 0,003632 м/с = 3,632 мм/с.
6. Проверяем Re_p: 0,50·0,003632·10⁻⁵ / (3,0·10⁻⁵) ≈ 0,000605. Условие выполнено.

Для проверки силы сопротивления в тех же условиях: F_с = 3·3,1416·3,0·10⁻⁵·10⁻⁵·0,003632 ≈ 1,027·10⁻¹¹ Н — практически тот же вес с учётом малого выталкивания.

• При d = 10 мкм итог: wос ≈ 3,63 мм/с; Reₚ ≈ 0,000605.
• При d = 20 мкм — wос ≈ 0,0145 м/с, или 14,5 мм/с. Reₚ ≈ 0,00484.

Удвоили диаметр — получили вчетверо большую скорость. Оба результата укладываются в выбранное приближение. Если у десятимикронной частицы уменьшить условную плотность до 500 кг/м³, скорость будет уже около 0,907 мм/с. Так становится понятно, почему две одинаковые по внешнему размеру частицы, плотная минеральная и лёгкая углеродистая, отделяются по-разному.

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

Успеет ли частица добраться до стенки

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

wᵣ ≈ (ρp − ρg)·d²·aц / (18μ).

Это приближённая скорость смещения относительно газа к периферии, м/с; все плотности, d и μ имеют те же единицы, a_ц берётся в м/с². Разность плотностей учитывает также действие градиента давления вращающегося газа. Время прохождения небольшого радиального пути можно оценить как t ≈ Δr/wᵣ: Δr — путь в метрах, t — время в секундах. Это обычное отношение пути к скорости при её приблизительном постоянстве на данном отрезке.

Оставим свойства предыдущего примера и a_ц = 1000 м/с². Пусть нужно пройти последние 5 мм, а доступное время условно составляет 0,030 с:

Диаметр частицы	Оценка wr	Rep	Время на 5 мм
5 мкм	0,0926 м/с	0,00771	0,0540 с
10 мкм	0,370 м/с	0,0617	0,0135 с

Проверим вторую строку вручную: w_r = 1999,5·10⁻¹⁰·1000 / (5,4·10⁻⁴) ≈ 0,3703 м/с. Затем Re_p = 0,50·0,3703·10⁻⁵ / (3,0·10⁻⁵) ≈ 0,0617. Наконец, t = 0,005 / 0,3703 ≈ 0,01350 с. У пяти микрометров квадрат диаметра вчетверо меньше, поэтому скорость вчетверо меньше, а время на тот же путь вчетверо больше. Re_p уменьшается ещё и за счёт самого диаметра, всего в восемь раз.

В этой модели первая частица не успевает, вторая успевает. Предполагается, что существенное скольжение направлено радиально; если есть ещё заметное осевое или окружное отставание, Re_p проверяют по полной относительной скорости.

А теперь полезная ловушка. Для частицы 20 мкм та же подстановка дала бы 1,48 м/с и Re_p около 0,494. Наше условие медленного обтекания уже не выполнено. В гравитационном примере формула подходила, в более сильном поле для той же частицы требует пересмотра. Поэтому «размер маленький, значит Стокс» — недостаточное обоснование.

Реальный циклон содержит целое распределение траекторий и времён пребывания. До стенки ещё нужно добраться, после этого — не сорваться обратно и выйти через сборник. Оценка выше объясняет роль размера и времени, но характеристику аппарата получают для его геометрии и рабочего потока. В частности, d₅₀ — размер, для которого фракционная эффективность составляет 50% в заданных условиях. Это не самый крупный размер пыли, которая вообще способна проскочить.

Расход меняет сразу несколько вещей

Средняя скорость газа на входе равна uвх = Qраб/Aвх, где Q_раб — рабочий объёмный расход в м³/с, A_вх — площадь входного сечения в м²; u_вх получается в м/с. Например, для условных 0,10 м³/с и 0,0050 м²: u_вх = 0,10 / 0,0050 = 20 м/с. Это только пример деления расхода на площадь, а не рекомендуемая скорость. Это другая скорость, чем v_θ частицы в вихре.

Для Q_раб нужен весь газовый поток, включая водяной пар, при температуре и абсолютном давлении в самом циклоне. Если исходно известен приведённый сухой расход, сначала восстанавливают влажную основу, затем пересчитывают объём к рабочему состоянию. Иначе можно подобрать вход под один объём, а пропускать через него другой.

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

Сквозной пример: сколько пыли останется в газе

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

• Q = 1000 м³/ч — полный рабочий влажный расход при неизменных условиях на входе и выходе;
• cвх = 10 г/м³ — концентрация твёрдых частиц на той же объёмной основе;
• массовое распределение пыли и фракционная эффективность — в таблице ниже.

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

Фракция по размеру	Массовая доля во входной пыли	Заданная эффективность
Более 50 мкм	40%	95%
От 10 до 50 мкм включительно	35%	70%
Менее 10 мкм	25%	20%

Сумма массовых долей: 40 + 35 + 25 = 100%. Здесь считаются доли массы, а не числа частиц.

Шаг 1. Входная масса за час. Массовый расход пыли M_вх = Q·c_вх:

1000 м³/ч · 10 г/м³ = 10 000 г/ч = 10 кг/ч.

Шаг 2. Разбиваем вход по фракциям. Для каждой M_i = M_вх·f_i, где f_i — её массовая доля в единицах. Перевод процентов в долю — деление на 100: 40% = 0,40. Получаем 10·0,40 = 4,0 кг/ч; 10·0,35 = 3,5 кг/ч; 10·0,25 = 2,5 кг/ч.

Шаг 3. Уловленное и проскок. M_ул,i = M_i·η_i, где η_i — заданная фракционная эффективность в долях. M_пр,i = M_i − M_ул,i. Для крупной фракции: 4,0·0,95 = 3,8 кг/ч уловлено; 4,0 − 3,8 = 0,2 кг/ч прошло дальше.

Фракция	Вход, кг/ч	Уловлено, кг/ч	Прошло дальше, кг/ч
Более 50 мкм	4,0	4,0·0,95 = 3,80	0,20
10–50 мкм	3,5	3,5·0,70 = 2,45	1,05
Менее 10 мкм	2,5	2,5·0,20 = 0,50	2,00
Итого	10,0	6,75	3,25

Проверка баланса: 6,75 + 3,25 = 10,0 кг/ч. Все массовые величины здесь относятся к твёрдым частицам.

Шаг 4. Общая эффективность. η_общ = M_ул/M_вх = 6,75/10,0 = 0,675, то есть 67,5% по массе. Это взвешенный результат для данного распределения. Обычное среднее трёх процентов было бы неправильным: фракции приносят неодинаковую массу.

Шаг 5. Выходная концентрация. При принятом неизменном расходе газа c_вых = M_пр/Q:

3,25 кг/ч · 1000 г/кг / 1000 м³/ч = 3,25 г/м³.

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

Что получилось по существу? Крупную пыль аппарат снимает хорошо, но 2,0 из оставшихся 3,25 кг/ч приходится на мелкую фракцию. Именно с ней в основном предстоит работать следующей ступени. Для сравнения с требованиями потребителя нужна ещё их допустимая концентрация на согласованной основе. А расчёт 6,75 кг/ч уловленного остатка уже задаёт нагрузку на сбор и выгрузку — этот материал не должен бесконечно храниться в конусе.

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

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

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

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

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

Какие данные стоит получить до большого аппарата

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

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

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

Так исследования и расчёт переходят в опытно-конструкторскую работу, а затем — в масштабирование. Если эту середину пропустить, роль стенда получает уже большая установка у заказчика. И обнаруженная липкая пробка становится не небольшим изменением корпуса в своей мастерской, а переделкой горячего тракта с ограниченным доступом и производством, которое ждёт газ. Такая доводка может занять годы; её масштаб плохо сочетается с первоначальной идеей сэкономить на испытаниях.

К оглавлениюALVERJOB

Источники и материалы

Источники и связи

• US EPA / Public Health Service. Control Techniques for Particulate Air Pollutants, AP-51, §§4.3.1–4.3.2, с. 4-15–4-17: инерционное отделение и силовой масштаб.
• University of Washington, M. Benjamin. Particle Settling / Sedimentation, с. 7–9; US EPA, APTI Course 413, гл. 3: сопротивление, осаждение, Re и границы приближений. Численные примеры рассчитаны для явно заданных условных величин; это не испытания циклона.
• US EPA. Cyclones, EPA-452/F-03-005, с. 1, 3–4: свойства частиц, сопротивление, липкие материалы, подсос и вторичный унос.
• R. Zwart и др. Tar removal from low-temperature gasifiers, ECN-E—10-008, §§2.2.1, 3.1, 3.2.1, с. 18, 29–31: разделение ступеней и отложения при холодном пуске.
• Huffman, Halligan, Peterson. Conversion of Cattle Feedlot Manure to Ethylene and Ammonia Synthesis Gas, EPA-600/2-78-026, с. 40–41: нагретый циклон и отдельное последующее отделение конденсируемых веществ в исследованной установке.
• Связи: C0089 — рабочая и приведённая основы расхода; C0095 — состав потока в точке отбора; C0098/C0101 — тонкая фильтрация; C0100 — мокрая очистка.