Расчёт кристаллизатора: выход и управление пересыщением
Расчёт кристаллизатора — прежде всего упражнение по балансу массы: в разобранном ниже примере один проход охлаждения с 60 до 25 °C извлекает 48,8% растворённого сульфата меди в виде кристаллов CuSO4·5H2O. Выход задают кривая растворимости и баланс; размер кристалла задают управление пересыщением и время пребывания. Ниже приведены методика расчёта, рабочие окна и режимы отказа, которые отделяют предприятие, выполняющее цель по извлечению, от того, которое никак её не достигает.
Ключевые выводы
- Выход кристаллов определяется балансом массы и кривой растворимости, а не выбором оборудования.
- Пересыщение, выраженное через ΔC или относительное пересыщение σ, — движущая сила и нуклеации, и роста.
- σ = 0,2–2% — практическое стационарное окно для большинства систем неорганических солей.
- В кристаллизаторах MSMPR средний размер слабо зависит от времени пребывания; управление мелочью — более сильный рычаг по размеру.
- Большинство провалов по выходу связано с гидратными переходами и неизмеренной кинетикой; и то и другое закрывается в пилотных испытаниях.
Основы: растворимость, пересыщение и выход
Кривая растворимости — граница, разделяющая ненасыщенный и пересыщенный раствор в равновесии: для каждой температуры она даёт концентрацию C* растворённого вещества (в расчёте на единицу свободной воды, то есть воды, остающейся растворителем, без воды, связанной в кристаллической решётке), при которой раствор ровно насыщен. Для сульфата меди пентагидрата C* растёт примерно с 0,207 кг безводного CuSO4 на кг свободной воды при 20 °C до примерно 0,400 при 60 °C. Для сульфата натрия кривая прерывается при 32,4 °C, где стабильная твёрдая фаза меняется с мирабилита (Na2SO4·10H2O) на тенардит (безводный Na2SO4) — переход, который уже ставил в тупик не одну технологическую схему.
Пересыщение — термодинамическое состояние, существующее, когда фактическая концентрация C превышает равновесное значение C* при той же температуре. Это единственная движущая сила и для нуклеации, и для роста. Оно всегда переходно в том смысле, что система его расходует: каждый грамм соли, осаждённый на поверхности кристалла, снижает C обратно к C*. Практическое проектирование поэтому состоит в том, чтобы создавать пересыщение с контролируемой скоростью и уничтожать его через поверхность кристаллов.
Ширина метастабильной зоны (MSZW) — интервал концентрации или температуры от кривой растворимости до точки начала спонтанной (первичной) нуклеации. Ниже MSZW чистый раствор может оставаться пересыщенным сколь угодно долго без нуклеации; внутри неё нуклеация становится вероятной, а затем бурной. Растворы с затравкой можно вести при ΔC значительно ниже предела MSZW — именно так работает промышленное оборудование: метастабильная зона это рабочая комната, а порог нуклеации — её стена.
Наконец, гранулометрический состав (CSD) — статистика популяции, описывающая распределение кристаллической массы по классам крупности. Выход говорит, сколько тонн вы продаёте; CSD говорит, будет ли продукт центрифугируемым и сыпучим материалом или мелким шламом, который кольматирует ткань и вносит маточный раствор в продукт. Весь расчёт ведётся на сухую безводную основу — килограммы безводного вещества на килограмм свободной воды, — потому что гидратные кристаллы забирают и воду, и соль, а небрежный учёт кристаллизационной воды — самая частая арифметическая ошибка в этой области.
Расчёт выхода: метод баланса массы
Выход стадии кристаллизации следует из одного баланса растворённого вещества и одного баланса воды. Поскольку гидрат забирает и то и другое, балансы связаны и решаются совместно.
Шаг 1 — зафиксировать основу. Выберите удобное количество исходного раствора (здесь 1 000 кг) и представьте его на сухой основе: масса свободной воды и масса растворённого безводного вещества.
Шаг 2 — взять растворимость при обеих температурах. По кривой прочитайте C1 (температура питания) и C2 (конечная температура) в кг безводного вещества на кг свободной воды. Для примера: CuSO4 при 60 °C C1 ≈ 0,400; при 25 °C C2 ≈ 0,230.
Шаг 3 — вычислить гидратное отношение. R = M(гидрат)/M(безводное). Для CuSO4·5H2O: 249,7 против 159,6 г/моль, R = 1,5646. Один килограмм продукта содержит 1/R = 0,6392 кг безводного CuSO4 и 0,3608 кг кристаллизационной воды.
Шаг 4 — решить связанный баланс. Для прохода охлаждения без испарения (V = 0) выход гидратных кристаллов составляет Y = S1·(C1 − C2) / [1/R − C2·(1 − 1/R)], где S1 — начальная масса свободной воды. Знаменатель учитывает воду, которую гидрат выводит из оборота. Если проход дополнительно испаряет V кг свободной воды, замените числитель на C1·S1 − C2·(S1 − V).
Шаг 5 — перевести в извлечение. Извлечение растворённой соли = (Y/R)/(C1·S1) в процентах.
Шаг 6 — проверить оба баланса. Просуммируйте безводное вещество и свободную воду по продукту и маточнику; каждый должен сойтись с питанием. Если не сходится водный баланс, неверен член кристаллизационной воды.
Разобранный пример. Основа: 1 000 кг очищенного раствора при 60 °C с 28,57% масс. безводного CuSO4 (C1 = 0,400; S1 = 714,3 кг). Охлаждение до 25 °C при C2 = 0,230, без испарения и уноса, продукт — чистый CuSO4·5H2O. Числитель: 714,3 × (0,400 − 0,230) = 121,4. Знаменатель: 0,6392 − 0,230 × 0,3608 = 0,5562. Y = 121,4 / 0,5562 = 218,3 кг кристаллов CuSO4·5H2O. Извлечение = (218,3 × 0,6392) / 285,7 = 139,5 / 285,7 = 48,8% растворённого сульфата меди.
| Поток | Свободная вода (кг) | Безводный CuSO4 (кг) | Итого (кг) |
|---|---|---|---|
| Исходный раствор, 60 °C | 714,3 | 285,7 | 1 000,0 |
| Кристаллы продукта | 78,8 (гидратная вода) | 139,5 | 218,3 |
| Маточник, 25 °C | 635,5 | 146,2 | 781,7 |
Оба баланса сходятся (714,3 = 78,8 + 635,5; 285,7 = 139,5 + 146,2; 218,3 + 781,7 = 1 000,0), а концентрация маточника воспроизводит C2: 146,2 / 635,5 = 0,230. Чувствительность: доведение того же раствора до 20 °C (C2 ≈ 0,207) даёт Y = 714,3 × 0,193 / 0,5645 = 244,3 кг, извлечение 54,7% — примерно на шесть процентных пунктов больше за 5 K дополнительного охлаждения. Числа, на которых стоит расчёт, хороши ровно настолько, насколько хороша кривая растворимости за ними, а на заводских растворах — с железом, кислотой и магнием — кривая смещается; поэтому литературную растворимость мы считаем начальной оценкой, а не проектным входом.
Как создаётся и измеряется пересыщение
Поскольку выход берётся из баланса, единственное, что остаётся регулировать, — это скорость и путь создания и расходования пересыщения. Абсолютное пересыщение ΔC = C − C*, относительное σ = ΔC/C*. Микропример: при 50 °C C* для CuSO4 ≈ 0,336 кг/кг; раствор 0,340 имеет ΔC = 0,004 кг/кг и σ = 0,004/0,336 = 1,2%. Промышленная практика для неорганических солей держит σ в полосе 0,2–2% в стационарном режиме: достаточно для приемлемых скоростей роста и достаточно мало, чтобы не касаться стены нуклеации.
| Метод | Как создаётся ΔC | Типичное σ в стационаре | Типичная гидрометаллургическая задача |
|---|---|---|---|
| Поверхностное охлаждение | Растворимость падает с температурой в теплообменнике | 0,5–2% | CuSO4·5H2O, MgSO4·7H2O, квасцы |
| Испарение | Растворитель удаляется при почти постоянной температуре | 0,2–1% | Na2SO4, MgCl2, никелевые соли с плоской кривой |
| Вакуумное (адиабатическое) охлаждение | Вскипание уносит явную теплоту и растворитель | 0,5–2% | Мирабилит, бура |
| Высаливание | Смешиваемый антирастворитель снижает C* | 1–5%, кратковременно | Спирто-рассольные системы в рафинировании Ni–Co |
| Реакция | Вещество образуется in situ быстрее, чем может осаждаться | 1–5%, переходно | Гипс, Mg(OH)2, осаждение никеля |
Для испарительного пути оборудованием выбора обычно служит испаритель-кристаллизатор с принудительной циркуляцией, где контур циркуляции держит поверхность теплообмена обеднённой пересыщением, а ΔC развивается в зоне кипения. Для чувствительных к размеру продуктов на охлаждении кристаллизатор Oslo держит кипящий слой крупных кристаллов в слабо пересыщенном растворе, обменивая производительность на качество CSD.
MSZW заслуживает ещё одного абзаца, потому что прямо превращается в политику управления. Измеренные значения для обычных сульфатных систем лежат примерно в 5–25 K охлаждения при скоростях 5–20 K/ч без затравки (величины порядка литературных, зависят от системы и примесей). Затравка может вдвое сократить используемый запас. Суждение, которое даёт только практика: мы вводим затравку примерно на трети измеренной MSZW, а не в точке термодинамического насыщения, и никогда не позволяем пусковым переходным процессам определять характер нуклеации прогона. Аппарат, плохо зародышеобразованный в первые тридцать минут, весь остаток кампании производит мелочь.
Баланс популяций MSMPR: нуклеация против роста
Большинство непрерывных промышленных кристаллизаторов хорошо описывается идеализацией MSMPR. В её предположениях баланс популяций даёт числовую плотность кристаллов n(L) = n⁰·exp(−L/Gτ), где n⁰ — плотность нуклеации, G — линейная скорость роста, τ — время пребывания продукта.
Поведением управляют два кинетических закона. Рост подчиняется G = kg·σ^g с g обычно 1–2; нуклеация (преимущественно вторичная, контактная на существующих кристаллах) — B⁰ = kb·σ^b с b обычно 2–4, то есть нуклеация является более крутой функцией пересыщения. Это неравенство показателей — центральный факт управления пересыщением: когда σ растёт, новые кристаллы появляются быстрее, чем растут существующие, и та же масса соли распределяется по большему числу более мелких кристаллов.
Классическое следствие для размеров: средний по массе размер L43 = 3,67·G·τ. Микропример: при G = 1,0 × 10⁻⁸ м/с и τ = 6 ч, L43 = 3,67 × 1,0 × 10⁻⁸ × 21 600 = 0,79 мм. Удвоит ли время пребывания до 12 ч этот размер? Нет, и баланс популяций объясняет почему: при фиксированной плотности суспензии баланс массы требует B⁰·G³·τ⁴ = const; подстановка кинетических законов даёт σ ∝ τ^(−4/(b+3g)), а при b = 2g получается L43 ∝ τ^0,2. Удвоение τ повышает L43 лишь в 2^0,2 = 1,15 раза, до ≈ 0,91 мм. Это второе правило проектирования по размеру: крупные аппараты покупают удивительно мало; сильные рычаги — снижение нуклеации и удаление мелочи.
| Параметр | Обозн. | Типичный диапазон (неорг. соли) | Эффект при росте |
|---|---|---|---|
| Время пребывания продукта | τ | 2–12 ч | Больше L43, но слабо (∝ τ^0,2 при b = 2g); больше капитала на тонну |
| Линейная скорость роста | G | 10⁻⁹–10⁻⁷ м/с | Больше производительности при том же σ |
| Скорость нуклеации | B⁰ | 10⁵–10⁸ шт./(м³·с) | Больше мелочи, меньше L43 |
| Доля удаляемой мелочи | Rf | 5–20% циркулирующей магмы | Растворённая мелочь возвращается ростом; CSD крупнее и уже |
| Плотность суспензии | MT | 50–400 кг/м³ | Больше площади на единицу ΔC; умеряет σ, но может усилить контактную нуклеацию |
Использовать ли эти рычаги в непрерывной или периодической схеме, зависит от масштаба, изменчивости питания и требований по размеру. Для целей проектирования практическая ценность MSMPR — в пропорциональном мышлении: он говорит, какая ручка двигает CSD и с какой силой, ещё до подгонки пилотных данных.
Самые частые ошибки проектирования
Они повторяются на предприятиях по солям; каждый пункт дан как симптом — причина — решение.
- Выход рассчитан только по тепловому балансу. Симптом: извлечение на заводе не дотягивает до проекта на 5–15 процентных пунктов. Причина: использована растворимость в чистой воде либо проигнорированы кристаллизационная вода и гидратные переходы. Решение: пересчитать баланс по кривой растворимости на реальном составе раствора, на сухой основе, с явным учётом воды кристаллизации.
- Игнорирование гидратных переходов. Симптом: завод сульфата натрия выдаёт шламообразную мелочь ниже 32,4 °C или выход исчезает в тёплый сезон. Причина: переключение мирабилит/тенардит меняет растворимость, R и потребность в воде скачком. Решение: осознанно выбрать температурную полосу — холодная кристаллизация мирабилита ради высокого выхода (высокий R) либо горячий тенардит ради сыпучего продукта — и проектировать тепловую схему под выбранную фазу.
- Слишком низкая продувка мелочи. Симптом: средний размер продукта сползает неделя за неделей, ткань центрифуги кольматируется. Причина: контур удаления мелочи занижен или его уставка никогда не настраивалась по реальной нуклеации. Решение: измерить CSD, затем настроить продувку так, чтобы контур снимал достаточную площадь и переносил рост на выжившие кристаллы.
- Обрастание поверхностей теплообмена из-за форсирования ΔC. Симптом: тепловая нагрузка падает за дни, циклы чистки укорачиваются. Причина: пересыщение создаётся у стенки, а не в объёме. Решение: снизить движущую силу, поднять скорость циркуляции и перенести испарение от теплообменника в зону кипения.
- Дрейф размера после изменения питания. Симптом: CSD мельчает после модернизации выше по схеме, повысившей содержание металла. Причина: большая производительность при той же геометрии требует большего ΔC, а B⁰ реагирует круче G, поэтому побеждает нуклеация. Решение: повторить проверку баланса популяций при новом содержании и добавить затравку либо мощность разрушения мелочи.
- Отсутствие измерения MSZW до детального проектирования. Симптом: аппарат бурно нуклеирует при пуске, недели продукта вне спецификации. Причина: приняты литературные данные MSZW, а примеси заводского раствора её сузили. Решение: измерить MSZW на реальном растворе, с затравкой и без, при представительных скоростях охлаждения.
И снова суждение из практики: большинство этих отказов — не ошибки расчёта, а ошибки учёта: вода гидрата, оборотные потоки, растворимость, сдвинутая примесями. Любая экспертиза проекта должна начинаться и заканчиваться балансом массы, который ничего не прощает и ничего не забывает.
Стоимость и энергия: испарять или охлаждать
Энергия разделяет два пути сильнее, чем капитал. Испарение должно подвести скрытую теплоту воды: около 2 260 кДж/кг при 100 °C и около 2 407 кДж/кг при 40 °C под вакуумом. Однокорпусной испаритель потребляет примерно 1,1–1,3 кг пара на кг испарённой воды; трёхкорпусная схема — 0,35–0,45 кг/кг; установка механической рекомпрессии пара (MVR) — 20–45 кВт·ч на тонну испарённой воды при эффективном ΔT 4–8 K, повышении температуры кипения (BPE) 3–8 K для насыщенных сульфатных растворов и КПД компрессора 75–80%. Для контекста: 30 кВт·ч/т — это лишь 4–5% теплосодержания той же задачи на пару, и именно поэтому MVR доминирует в новых испарительных проектах при наличии электроэнергии.
Кристаллизация охлаждением термодинамически намного дешевле. Нагрузка — явная теплота плюс теплота кристаллизации; для CuSO4·5H2O последняя ≈ 78 кДж/моль, или ≈ 0,31 МДж на кг гидрата. В примере: явная нагрузка ≈ 1 000 кг × 3,5 кДж/(кг·K) × 35 K = 122,5 МДж, теплота кристаллизации ≈ 218,3 × 0,31 = 67,7 МДж, итого холодильник видит ≈ 190 МДж на тонну питания — примерно 53 кВт·ч, сопоставимо на тонну с нагрузкой MVR, но в виде низкопотенциального холода, а не работы вала. Экономика выбора поэтому держится на наклоне растворимости: крутые кривые (сульфат меди) вознаграждают охлаждение, плоские заставляют испарять.
Оговорка по BPE: поскольку она сжимает эффективный ΔT и для теплообменника, и для компрессора MVR, её недооценка на концентрированном растворе — классическая причина заниженной испарительной мощности. Ни одна из этих цифр не является коммерческим предложением; это диапазоны проектного контекста.
Каких выходов ждать для меди, никеля и сульфата натрия
Три иллюстративных случая по одной и той же методике и все явно помечены как иллюстративные (растворимости порядка литературных, чистые растворы, без оборота).
Сульфат меди (охлаждение): как разобрано выше, 1 000 кг раствора 28,57% масс., охлаждённые с 60 до 25 °C, дают 218,3 кг CuSO4·5H2O и 48,8% извлечения; углубление до 20 °C поднимает его до 54,7%. Крутая кривая делает сульфат меди одной из немногих гидрометаллургических солей, где одно охлаждение извлекает половину растворённой ценности металла.
Сульфат никеля (испарение): основа 1 000 кг очищенного раствора 30% масс. NiSO4 (сухая основа), испарить 250 кг свободной воды, закончить насыщением при 30 °C с иллюстративным C2 = 0,50 кг/кг. При R = 262,85/154,75 = 1,699 для гексагидрата: Y = (0,4286 × 700 − 0,50 × 450) / (0,5887 − 0,50 × 0,4112) = 75 / 0,3832 = 195,7 кг NiSO4·6H2O, то есть 38,4% растворённого никеля. Умеренно плоская кривая — причина, по которой заводы сульфата никеля испарительные.
Сульфат натрия (вакуумное охлаждение): основа 1 000 кг раствора 25% масс. Na2SO4, охлаждение с 40 до 20 °C без чистого испарения. При R = 322,20/142,04 = 2,269 для мирабилита: Y = 750 × (0,3333 − 0,195) / (0,4409 − 0,195 × 0,5591) = 103,7 / 0,3318 = 312,7 кг Na2SO4·10H2O, извлечение 55,1%: высокое гидратное отношение (более половины кристалла — вода) поднимает выход. Любое мгновенное испарение в вакуумном охладителе добавляет извлечение сверх этого пола.
Обратите внимание на разброс: 38–55% извлечения за один проход по одной методике — вот почему расчёт кристаллизатора всегда специфичен для соли.
Частые вопросы
Чем отличаются абсолютное и относительное пересыщение?
ΔC — абсолютная движущая сила: разность фактической концентрации и равновесной растворимости C* при той же температуре, в кг/кг свободной воды. Относительное пересыщение σ = ΔC/C* выражает эту разность долей растворимости, и именно с ним коррелируют скорости нуклеации и роста. ΔC используйте для учёта и выбора размеров оборудования, σ — для сравнения систем и установления пределов регулирования.
Увеличивает ли затравка выход кристаллов?
Нет. Выход задан балансом массы и кривой растворимости, поэтому затравка не может изменить, сколько соли уйдёт из раствора при равновесии. Она меняет как этот выход проявляется: расходует доступное пересыщение ростом на существующих поверхностях вместо спонтанной нуклеации, давая меньше и крупнее кристаллов с более узким CSD. Количество, крупность и момент ввода затравки — среди самых чувствительных рабочих переменных.