Перепад давления слоя: почему вентилятор крутится, а газа больше не становится

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

GWEB_C0058

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

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

На щите было 0,23 кПа, стало 0,45 кПа. Перепад почти удвоился. Можно уже начинать искать шлак, мелочь и виноватую смену? Пока рано: если одновременно расход вырос в полтора раза, для совершенно неизменного слоя получится почти такая же картина. Ниже это можно пересчитать на обычном калькуляторе.

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

Сначала найти оба конца измерения

Обозначим давление перед рассматриваемым участком p₁, после него p₂. Здесь «перед» и «после» относятся к направлению движения газа. Тогда положительный перепад по ходу потока:

Газ проходит через слой между отборами p₁ и p₂. Они подключены к измерению ΔP, а отдельный расходомер даёт Q с учётом температуры и абсолютного давления. — Два отбора ограничивают участок слоя; расход измеряют отдельно и приводят к его условиям. Схема не включает колосник и очистку.

ΔP = p₁ − p₂

Все три величины выражены в паскалях, Па. Для примера p₁ = −1000 Па, p₂ = −3000 Па относительно атмосферы. Получаем ΔP = −1000 − (−3000) = 2000 Па = 2 кПа. Оба отбора находятся под разрежением, но сопротивление участка положительное. 1 кПа = 1000 Па; 1 мбар = 100 Па. Миллиметр водяного столба соответствует приблизительно 9,81 Па, поэтому 2000 Па составляют около 204 мм вод. ст., а не 2000 мм.

На схеме стоит нарисовать отборы и подписать всё, что находится между ними. Если туда попали слой, колосник и отводящий патрубок, прибор показывает их совместный перепад. Приписывать весь прирост слою уже нельзя. Отдельный перепад очистки позволяет увидеть, где действительно изменилось сопротивление. Эта простая мысль есть и в старых руководствах по газостанциям: давление измеряли по участкам именно для поиска места затруднённого прохода газа. [1]

Для связи измерения с расчётом ниже предполагается, что отборы передают статическое давление, а изменения скорости на входе и выходе и гидростатическая поправка малы либо учтены отдельно. Если один отбор оказался в струе, а другой после расширения, разность показаний содержит ещё и последствия изменения потока. Табличка «ΔP слоя» на экране геометрию патрубков не исправляет.

Что делает газ между кусками топлива

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

Разделю расчёт на три действия:

ΔP = ΔP_в + ΔP_и

ΔP_в = 150 × μ × L × (1 − ε)² × u / (ε³ × d²)

ΔP_и = 1,75 × ρ × L × (1 − ε) × u² / (ε³ × d)

Здесь:

ΔP_в и ΔP_и — вязкая и инерционная составляющие перепада, Па. Индексы «в» и «и» служат подписями, а не дополнительными множителями.
L — макроскопическая длина участка в направлении среднего потока, м. Для вертикального осевого потока это высота засыпки между границами расчёта, а не суммарная длина извилистого пути между частицами.
ε — межчастичная порозность, доля объёма слоя, свободная от частиц. Величина безразмерная: ε = 0,40 означает 40% свободного объёма. В формулу подставляют 0,40, а не 40. Мелкие внутренние поры самого древесного угля не следует автоматически прибавлять к промежуткам между кусками.
d — диаметр частиц, м; в рассматриваемом простом варианте частицы одинаковые и сферические. 10 мм нужно записать как 0,010 м.
u — поверхностная, или приведённая к пустому сечению, скорость газа, м/с. Она равна расходу газа в условиях этого участка, делённому на полную площадь его поперечного сечения.
μ — динамическая вязкость газа, Па·с, то же самое кг/(м·с). Она характеризует вязкое сопротивление газа и зависит от температуры и состава. Это не кинематическая вязкость, измеряемая в м²/с.
ρ — плотность газа в этом участке, кг/м³. Плотность при нормальных условиях сюда вместо рабочей не подставляют.
150 и 1,75 — безразмерные эмпирические коэффициенты исходной корреляции. Они получены согласованием зависимости с опытными данными. Это не физические константы, обязательные для любой щепы, золы и укладки. В руководстве разработчиков Pronghorn этот предел прямо оговорён; запись двух слагаемых также приведена в материалах эксперимента NASA с насыпным слоем. [2, 3]

Связь расхода и скорости выглядит проще:

u = Q / A

Q здесь — местный объёмный расход, м³/с; A — полная площадь проходного сечения аппарата, м². Именно полная, без вычитания площади частиц. Если разделить Q ещё и на ε, получится другая скорость, относящаяся к свободному объёму. Повторно подставлять её в приведённую запись Эргуна нельзя: влияние порозности там уже учтено.

Размерности тоже можно проверить без специальной программы. У вязкой части μLu/d² имеет единицы (Па·с) × м × (м/с) / м² = Па. У инерционной части ρLu²/d получаем (кг/м³) × м × (м²/с²) / м = кг/(м·с²) = Па. Порозность и оба числовых коэффициента единиц не добавляют.

Где заканчивается формула и начинается настоящий слой

Для расчётного примера ниже специально взяты одинаковые шарообразные частицы, равномерная укладка, постоянное сечение, один газ и неизменные свойства по длине. Перепад мал относительно абсолютного давления. В таком случае можно пользоваться одним значением μ, ρ, u и просто умножать удельное сопротивление на L.

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

С размером несферических частиц отдельная история. Длинная тонкая щепка и округлый кусок того же ситового класса укладываются по-разному. Меняются свободный объём, ориентация и площадь поверхности, которую обтекает газ. Для некоторых засыпок используют эквивалентный диаметр по наружной поверхности и объёму частиц, поправки на форму либо коэффициенты, определённые испытанием. Простое среднее между длиной, шириной и толщиной щепки не становится от этого проверенным гидравлическим диаметром. Исследования Nemec и Levec как раз сравнивали форму частиц и способы укладки; их результаты не дают основания объявить один коэффициент универсальным для топлива. [4]

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

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

Считаем один слой, а затем меняем одну причину

Рассмотрю учебный однородный участок: числа не описывают авторскую установку и не задают размеры рекомендуемого аппарата.

L = 0,80 м; A = 0,20 м²;
d = 0,010 м; ε = 0,40;
ρ = 0,40 кг/м³; μ = 0,000030 Па·с;
Q = 360 м³/ч именно при условиях слоя;
абсолютное давление условно около 100 кПа; выбранные свойства газа остаются постоянными во всём примере.

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

Сначала переводим расход: 360 / 3600 = 0,100 м³/с. Число 3600 — количество секунд в одном часе. Затем u = 0,100 / 0,20 = 0,500 м/с.

Подготовим повторяющиеся действия: 1 − ε = 0,60;(1 − ε)² = 0,36;ε³ = 0,064;d² = 0,000100 м²;u² = 0,250 м²/с².

Вязкая часть:

ΔP_в = (150 × 0,000030 × 0,80 × 0,36 × 0,500) / (0,064 × 0,000100)

ΔP_в = 0,000648 / 0,0000064 = 101,25 Па

Инерционная часть:

ΔP_и = (1,75 × 0,40 × 0,80 × 0,60 × 0,250) / (0,064 × 0,010)

ΔP_и = 0,084 / 0,00064 = 131,25 Па

Складываем: 101,25 + 131,25 = 232,50 Па = 0,2325 кПа. Это около 0,23% от заданных 100 кПа абсолютного давления. Пренебрежение изменением плотности только из-за этого небольшого перепада в нашем примере согласовано с исходным упрощением.

Теперь оставим слой прежним, а расход увеличим до 540 м³/ч. Скорость станет 540 / 3600 / 0,20 = 0,750 м/с, то есть вырастет в 1,5 раза. Вязкая часть также увеличится в 1,5 раза; инерционная — в 1,5² = 2,25 раза:

ΔP = 101,25 × 1,5 + 131,25 × 2,25 = 151,875 + 295,3125 = 447,1875 Па

Округлённо 447,2 Па. Вот и почти удвоившийся перепад из начала текста. Слой в этом расчёте вообще не меняли.

Во втором опыте вернём 360 м³/ч, но уменьшим ε с 0,40 до 0,35. Это именно изменение структуры при прежней высоте, а не расчёт осадки одной и той же массы. Для реальной осадки изменится и L, что пришлось бы учесть отдельно. Новый свободный объём меньше, а оставшиеся проходы теснее. Вязкая часть меняется в отношении:

[(0,65² / 0,35³) / (0,60² / 0,40³)] = 1,75186265

Инерционная часть меняется в отношении:

[(0,65 / 0,35³) / (0,60 / 0,40³)] = 1,61710398

Следовательно, ΔP_в = 101,25 × 1,75186265 = 177,38 Па; ΔP_и = 131,25 × 1,61710398 = 212,24 Па. Сумма — 389,62 Па. Порозность уменьшилась всего на пять процентных пунктов, а сопротивление при прежнем расходе выросло примерно на 68%.

Если вместо порозности изменить только d с 10 до 7 мм, вязкую часть нужно умножить на (10/7)², а инерционную на 10/7. Получится 206,63 + 187,50 = 394,13 Па. Это отдельный условный опыт, в котором порозность удержана прежней. Для настоящей мелочи оба изменения могут произойти одновременно.

Учебный вариант
Исходный слой
Местный расход, м³/ч
360
Порозность
0,40
Диаметр, мм
10
Перепад, Па
232,50

Учебный вариант
Только больше расход
Местный расход, м³/ч
540
Порозность
0,40
Диаметр, мм
10
Перепад, Па
447,19

Учебный вариант
Только меньше порозность
Местный расход, м³/ч
360
Порозность
0,35
Диаметр, мм
10
Перепад, Па
389,62

Учебный вариант
Только мельче частицы
Местный расход, м³/ч
360
Порозность
0,40
Диаметр, мм
7
Перепад, Па
394,13

Учебный вариант
Только половина расхода
Местный расход, м³/ч
180
Порозность
0,40
Диаметр, мм
10
Перепад, Па
83,44

Последнюю строку тоже легко проверить: 101,25 × 0,5 + 131,25 × 0,5² = 83,4375 Па. Зависимость здесь не чисто квадратичная: сохраняется и линейное слагаемое. Поэтому привычное деление перепада на квадрат расхода не даст одинакового числа на всех нагрузках даже у неизменного слоя.

А вентилятор согласен с новым расходом?

Пока мы задавали расход извне, словно кто-то обязан его обеспечить. Настоящий вентилятор работает на пересечении своей характеристики и характеристики сети. При фиксированных оборотах увеличение сопротивления обычно смещает рабочую точку к меньшему расходу; величина изменения зависит от формы обеих кривых. Если автоматика поддерживает расход повышением оборотов, уменьшения расхода сначала может не быть, зато изменятся обороты, потребляемая мощность и запас по давлению. Общие принципы такой проверки изложены в руководстве DOE по вентиляторным системам. [5]

Продолжим исключительно учебный расчёт. Обозначим q = Q / Q₀, где Q₀ = 360 м³/ч, а q — безразмерное отношение. Все характеристики относятся к одному местному объёмному расходу; охлаждения и подсосов между элементами в этой модели нет. Перепад исходного слоя уже известен:

ΔP_сл = 101,25q + 131,25q², Па

Зададим сопротивление остальной сети ΔP_ост = 200q² Па и характеристику вентилятора при неизменных оборотах ΔP_вент = 700 − 267,5q² Па. Числа 200, 700 и 267,5 имеют единицу Па, поскольку q безразмерно. Они выбраны для наглядного примера, не взяты из каталога. Индексы «сл», «ост», «вент» означают слой, остальную сеть и вентилятор. Давления всех характеристик согласованы между собой; различиями скоростных напоров в граничных сечениях пренебрегаем.

При q = 1 сеть требует 101,25 + 131,25 + 200 = 432,5 Па. Вентилятор развивает 700 − 267,5 = 432,5 Па. Равенство выполнено, расход равен 360 м³/ч.

После уменьшения порозности до 0,35 характеристика слоя станет 177,376q + 212,245q² Па. При прежнем расходе сеть уже требует 589,621 Па, а вентилятор всё ещё может дать 432,5 Па. Удержаться в прежней точке система не сможет.

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

700 − 267,5q² = 177,376q + 212,245q² + 200q²

Переносим слагаемые: 679,745q² + 177,376q − 700 = 0. По формуле положительного корня квадратного уравнения:

q = [−177,376 + √(177,376² + 4 × 679,745 × 700)] / (2 × 679,745) ≈ 0,89267

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

Новый расход Q = 0,89267 × 360 = 321,36 м³/ч. Перепад слоя получается около 327,47 Па, остальной сети — 159,37 Па; всего 486,84 Па. Проверка по вентилятору даёт те же 486,84 Па.

Значит, при ухудшении слоя мы увидим одновременно меньший расход и больший перепад слоя. Но измеренные 327 Па будут заметно ниже рассчитанных ранее 390 Па при принудительно сохранённых 360 м³/ч. Без расходомера степень изменения слоя легко недооценить: часть прироста сопротивления уже «оплачена» потерей производительности.

Для реальной установки кривую вентилятора пересчитывают на фактическую плотность у его входа и проверяют допустимую область работы, привод и ограничения изготовителя. После охладителя его объёмный расход может сильно отличаться от расхода горячего газа в слое. Поэтому брать оба значения из одной графы «м³/ч» без условий приведения нельзя. [5]

«Улитка» и Рутс по-разному отвечают на забитый проход

До сих пор мы рассматривали динамический вентилятор. У обычной центробежной «улитки» расход при данных оборотах определяется её характеристикой и сопротивлением сети. Форму кривой берут для конкретного колеса, а не для слова «улитка». Бококанальная вихревая воздуходувка устроена иначе; смешивать эти названия не стоит. В пособии изготовителя Greenheck показано и изменение рабочей точки, и необходимость оставаться в устойчивой области характеристики. [9]

Объёмная роторная воздуходувка типа Рутса переносит газ порциями между входом и выходом. При неизменных оборотах её расходная характеристика обычно жёстче: увеличение противодавления вызывает сравнительно меньшее уменьшение подачи. Поэтому при переменном сопротивлении дутья такой принцип привлекателен. Но «сравнительно постоянный расход» не означает «точно столько же воздуха при любых условиях». Подача зависит от оборотов, внутренних перетечек, состояния машины и условий на всасывании. AERZEN прямо отмечает увеличение перетечек с перепадом давления. [10]

Важное различие видно даже на предыдущем учебном примере. Динамический вентилятор после уплотнения слоя установился на 321,36 м³/ч при 486,84 Па. Если вместо него мысленно поставить идеальный источник прежнего расхода 360 м³/ч, то q останется равным единице, а требуемое давление станет:

ΔP_сети = 177,376 × 1 + 212,245 × 1² + 200 × 1² = 589,621 Па.

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

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

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

Измерять расход и регулировать его можно с обоими типами машин. Для процесса следует заранее выбрать, что поддерживается: массовая подача, объёмный расход при определённых условиях приведения или местный объёмный расход. При изменении температуры и давления одинаковые рабочие м³/ч содержат разную массу воздуха. Для отношения «воздух/топливо» это существенно; при переменной влажности воздуха уточняют и основу измерения. Обороты привода полезны в журнале, но расходомером сами по себе не становятся.

Воздуховод тоже умеет запасать подачу

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

Полезно оценить запас количества воздуха отдельно от скорости его выхода. Возьмём два условных варианта подводящего тракта: внутренний объём V = 0,20 м³ и V = 0,020 м³. Во всём учитываемом объёме зададим одинаковую температуру 60 °C; начальное абсолютное давление P₁ = 130 кПа, конечное P₂ = 110 кПа. Оба состояния считаем равновесными при одной температуре. Это учебные исходные числа, не разрешённые давления и не уставки защиты.

Количество освободившегося воздуха выразим объёмом при отдельно заданных опорных условиях: P_оп = 100 кПа, T_оп = 20 °C. Это наши условия сравнения, не обозначение «нормальных кубометров». Из уравнения идеального газа следует: при одинаковом составе количество вещества пропорционально PV/T. [12] Следовательно, избыточному количеству между двумя состояниями соответствует:

ΔV_оп = V × (P₁ − P₂) / P_оп × T_оп / T.

Здесь ΔV_оп — объём этого количества воздуха при опорных условиях, м³; V — геометрический внутренний объём тракта, м³; P₁ и P₂ — абсолютные давления в начале и конце сравнения; P_оп — абсолютное опорное давление. Все давления в отношении записываются в одинаковых единицах, здесь в кПа. T и T_оп — абсолютные температуры, К. Индекс «оп» означает опорные условия. Температуру в градусах Цельсия переводят в кельвины прибавлением 273,15. Газ считается идеальным; для значительного отклонения от идеальности понадобятся коэффициенты сжимаемости, здесь они приняты равными единице.

Сначала переводим температуры и считаем отношения:

T = 60 + 273,15 = 333,15 К;T_оп = 20 + 273,15 = 293,15 К.

(P₁ − P₂) / P_оп = (130 − 110) / 100 = 0,20.

T_оп / T = 293,15 / 333,15 ≈ 0,87993.

Для большого объёма:

ΔV_оп = 0,20 × 0,20 × 0,87993 ≈ 0,03520 м³ = 35,20 л.

Для малого:

ΔV_оп = 0,020 × 0,20 × 0,87993 ≈ 0,003520 м³ = 3,520 л.

Проверка проста: объём тракта уменьшили в десять раз, запас количества при тех же двух давлениях и температуре уменьшился в десять раз. Внутри после снижения давления газ остаётся; посчитана разность его количества, а не полное опорожнение трубы.

Для масштаба зададим обычную подачу 120 м³/ч именно при наших опорных условиях. Это 120 / 3600 = 0,03333 м³/с, или 33,33 л/с. Полученные 35,20 л соответствуют примерно 35,20 / 33,33 = 1,06 секунды такой подачи, а 3,520 л — 0,106 секунды. Это эквивалент количества, не длительность выброса. Пиковый расход и движение топлива по этому расчёту определить нельзя: нужны нестационарная модель прохода, давления в реакторе, теплообмен и динамика самого препятствия. При быстром расширении температура меняется, поэтому изотермическая оценка не заменяет расчёт переходного процесса.

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

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

Давление ограничивает система защиты, а не терпение двигателя

Нельзя считать, что Рутс безопасно повышает давление «пока хватает двигателя». Ограничения воздуходувки включают перепад, температуру и обороты; допуска корпуса, труб, соединений и реактора могут оказаться строже. В руководстве ROOTS TRI-RAM предел работы определяется первым достигнутым ограничением. Там же предусмотрены предохранительный клапан, разгрузка, контроль давления и температуры; использование предохранительного клапана как постоянного регулятора подачи прямо отвергается. [11]

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

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

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

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

Рабочий режим, форсаж и малая нагрузка

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

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

На малой нагрузке перепад падает, что видно и в последней строке таблицы. Тепловому процессу от этого не обязательно легче. Для нисходящих газогенераторов Imbert FAO описывает ограничение снизу по сохранению температуры горячей зоны и разложению смол, а сверху — по сопротивлению восстановительной зоны. Численные границы из того руководства относятся к его конструкциям и топливу; переносить их на другой аппарат не следует. [6]

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

Когда засорился путь к прибору

Отбор давления и импульсная линия являются частью измерения. Пыль и смола могут перекрыть сообщение с процессом; конденсат создаёт жидкостный столб, которого в расчёте сухой линии не было. Производитель датчиков Yokogawa отдельно рассматривает закупорку импульсных линий и изменение передаваемых колебаний давления. [7]

Если в одной линии появился некомпенсированный вертикальный столб воды высотой h = 0,020 м, его гидростатический вклад по модулю составит:

ΔP_ж = ρ_ж × g × h = 1000 × 9,80665 × 0,020 = 196,13 Па

ρ_ж — принятая для примера плотность воды, 1000 кг/м³; g — стандартное ускорение свободного падения, 9,80665 м/с²; h — разность высот, м. Индекс «ж» означает жидкость. В реальном конденсате плотность может быть другой. Знак ошибки зависит от того, в каком плече и как расположен столб. Речь именно о некомпенсированной разности столбов: две одинаковые и правильно учтённые жидкостные линии — другой случай.

Всего два сантиметра воды дают около 196 Па на фоне наших расчётных 232,5 Па. Поэтому спор о качестве топлива по такому прибору может надолго отвлечь от довольно мокрой причины.

Частичная закупорка часто замедляет отклик; полная способна оставить прибор связанным с запертым объёмом вместо текущего давления. При этом показание не обязано стать нулевым или постоянным: запертый газ меняет давление при изменении температуры. Полезно сопоставлять реакцию обоих давлений и перепада с уже произошедшим штатным изменением расхода, а при наличии независимых измерений — сверять их разность с каналом ΔP. Исчезновение пульсаций само по себе не доказывает закупорку, особенно при смене режима. Это ограничение подтверждено и экспериментами на воздушной линии. [8]

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

Что выяснить на опытном участке до большого аппарата

Самая дешёвая часть нашего примера — нажать на калькуляторе 0,35 в третьей степени. Получить представительное значение порозности для работающего слоя гораздо труднее. У одной и той же исходной щепы после нагрева меняются форма, прочность и количество мелочи. Измеренная плотность углеродного материала сама по себе ещё не даёт межчастичную порозность: внутренние поры куска и проходы между кусками необходимо различать. В расчётной ведомости полезно прямо разделить измеренные величины, справочные свойства газа и пока предположенные параметры засыпки.

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

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

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

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

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

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

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

Продолжить чтение

GWEB_C0075 — совместный разбор давления, температур, расхода и анализа газа на учебных журналах.
GWEB_C0069 — каналы в топливном слое и неравномерное распределение потока.
GWEB_C0018 — зола, спекание и изменение структуры слоя.

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

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

Источники и границы использования

Д. Б. Гинзбург. «Газификация твёрдого топлива», 1958, глава X, печатные с. 92–96. Разделение измерений по участкам и связь режима с несколькими показателями. Старые способы подключения и нормативные числа не перенесены.
NASA, материалы Packed Bed Reactor Experiment, слайд 8. Однофазная запись Эргуна и разделение вязкого и инерционного вкладов. Условия космического эксперимента не перенесены на газогенератор.
Novak и др. Pronghorn Theory Manual, INL/EXT-18-44453, §9.1.1.1, с. 99–100. Происхождение коэффициентов и отсутствие оснований считать их универсальными.
Nemec, Levec. Flow through packed bed reactors: 1. Single-phase flow, 2005. Использованы доступные аннотация и описание исследования формы и укладки. Полный текст не получен; частные поправки не воспроизводятся.
DOE. Improving Fan System Performance, разделы о регулировании расхода и оценке рабочей точки. Использованы принципы согласования сети и вентилятора и поправки на плотность; все кривые в примере придуманы специально для расчёта.
FAO. Wood gas as engine fuel, §2.5.2. Применимость обсуждения верхней и нижней нагрузки ограничена рассматриваемыми нисходящими газогенераторами.
Yokogawa. Impulse Line Blocking Diagnosis. Возможные нарушения передачи давления; гидростатический пример рассчитан отдельно.
Uehara и др. Diagnostics of Impulse Line Blockage with Multi-sensing Differential Pressure Transmitter at Air Line, 2009. Ограничения диагностики по пульсациям при переходных режимах.

Greenheck. Fan Curves, FA/100-99, с. 2–3: семейство характеристик, рабочая точка и область устойчивости. Конкретный вентилятор для газогенератора здесь не выбран.
AERZEN. How to select the most effective blower technology, PDF с. 2: объёмный принцип, регулирование оборотами и увеличение перетечек с перепадом давления. Экономические сравнения для очистных сооружений на газогенератор не перенесены.
ROOTS. TRI-RAM Installation, Operation and Maintenance, 19684-ISRB-2008-0124, печатные с. 2, 5, 7 и рис. 3 на с. 10. Первичное руководство производителя размещено у дистрибьютора. Использованы различие ограничений, защитные функции и запрет подменять регулятор предохранительным клапаном; численные уставки и монтажная схема не копируются.
NASA Glenn. Equation of State: уравнение идеального газа. Расчёт количества воздуха в двух объёмах — авторский учебный пример.

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