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

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

Допущения приближения псевдоустановившегося состояния

Тепловые профили устанавливаются быстрее, чем перемещается граница раздела

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

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

Скрытая теплота преобладает в балансе подвижной границы

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

Таким образом, баланс на границе раздела основывается на соотношении

[ \text{отведённое тепло}=\text{скрытая теплота, выделенная новой твёрдой фазой}. ]

Построение модели термических сопротивлений

Определите геометрию трубы и твёрдой фазы

Пусть

(R_i) — внутренний радиус трубы,
(R_t) — внешний радиус трубы,
(R_s(t)) — радиус подвижной границы раздела твёрдой и жидкой фаз,
(k_w) — теплопроводность стенки трубы,
(k_s) — теплопроводность слоя твёрдой фазы,
(h_i) — коэффициент теплоотдачи со стороны внутреннего хладагента,
(T_c) — температура хладагента в объёме,
(\rho_s) — плотность твёрдой фазы,
(L) — скрытая теплота затвердевания или кристаллизации.

Для осадка, растущего наружу от поверхности трубы, твёрдая фаза занимает область между (R_t) и (R_s), а её толщина равна

[ \delta(t)=R_s(t)-R_t. ]

Выразите сопротивление на единицу длины трубы

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

[ \mathcal{R}_i'=\frac{1}{h_i,2\pi R_i}, ]

[ \mathcal{R}_w'=\frac{\ln(R_t/R_i)}{2\pi k_w}, ]

[ \mathcal{R}_s'=\frac{\ln(R_s/R_t)}{2\pi k_s}. ]

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

Таким образом, суммарное сопротивление на единицу длины равно

[ \mathcal{R}_{\text{tot}}'

\frac{1}{h_i,2\pi R_i} + \frac{\ln(R_t/R_i)}{2\pi k_w} + \frac{\ln(R_s/R_t)}{2\pi k_s}. ]

Рассчитайте мгновенную скорость отвода тепла

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

[ q'(R_s)

\frac{T_m-T_c}{\mathcal{R}_{\text{tot}}'}. ]

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

Вывод скорости роста твёрдой фазы

Запишите баланс скрытой теплоты

Когда граница раздела перемещается наружу на величину (dR_s), объём вновь образовавшейся твёрдой фазы на единицу длины трубы равен

[ dV'=2\pi R_s,dR_s. ]

Соответствующая выделившаяся скрытая теплота равна

[ dQ_{\text{latent}}'

\rho_s L,2\pi R_s,dR_s. ]

Разделяя на время и приравнивая эту скорость к отводимому теплу, получаем

[ q'(R_s)

2\pi R_s\rho_s L\frac{dR_s}{dt}. ]

Подстановка выражения для сопротивления даёт основное уравнение роста:

[ \boxed{ \frac{dR_s}{dt}

\frac{T_m-T_c} { 2\pi R_s\rho_s L \left[ \frac{1}{h_i,2\pi R_i} + \frac{\ln(R_t/R_i)}{2\pi k_w} + \frac{\ln(R_s/R_t)}{2\pi k_s} \right] } } ]

Это уравнение задаёт мгновенную радиальную скорость роста слоя твёрдой фазы.

Интерпретируйте снижение скорости роста

По мере увеличения (R_s) сопротивление твёрдой фазы

[ \frac{\ln(R_s/R_t)}{2\pi k_s} ]

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

Также важен множитель (2\pi R_s) в члене, описывающем скрытую теплоту: при увеличении длины окружности границы раздела заданное радиальное приращение создаёт больший объём твёрдой фазы на единицу длины трубы.

Полезные упрощения для инженерных оценок

Используйте плоское приближение для тонкого слоя

Если толщина твёрдой фазы мала по сравнению с радиусом трубы, (\delta \ll R_t), цилиндрическое сопротивление можно приближённо заменить плоским сопротивлением. Уравнение роста принимает вид

[ \frac{d\delta}{dt}

\frac{T_m-T_c} {\rho_s L \left[ \frac{1}{h_i} + \frac{t_w}{k_w} + \frac{\delta}{k_s} \right]}, ]

где (t_w) — толщина стенки трубы.

Если сопротивления плёнки хладагента и стенки объединить в постоянное сопротивление (R_0), то

[ \frac{d\delta}{dt}

\frac{T_m-T_c} {\rho_s L\left(R_0+\delta/k_s\right)}. ]

Интегрируйте уравнение для оценки зависимости толщины от времени

При постоянных свойствах и разности температур интегрирование даёт

[ t

\frac{\rho_s L}{T_m-T_c} \left[ R_0\delta+\frac{\delta^2}{2k_s} \right]. ]

Этот результат показывает два предельных режима.

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

[ \delta \propto \sqrt{t}. ]

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

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

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

Применение модели к трубчатой технологической установке

Оцените эффективный коэффициент теплопередачи

Модель сопротивлений позволяет определить, какая часть системы ограничивает рост. Внутренняя плёнка хладагента описывается величиной (1/h_i), тогда как растущая твёрдая фаза в плоском приближении вносит вклад (\delta/k_s).

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

[ U=\frac{1}{R_{\text{total}}}. ]

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

Оцените рабочие условия хладагента

Температура и расход хладагента влияют на движущую силу теплопередачи и внутренний коэффициент (h_i). Снижение (T_c), как правило, увеличивает движущую силу (T_m-T_c), тогда как увеличение скорости хладагента может уменьшить сопротивление внутренней плёнки.

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

Учитывайте влияние температуры и состава

В системах кристаллизации (T_m) может изменяться в зависимости от концентрации, давления или оттеснения растворённого вещества на границе раздела. Если жидкость постепенно концентрируется, температура фазового перехода и доступная движущая сила могут изменяться в ходе процесса.

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

Свяжите модель с осевым поведением в трубе

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

[ \Delta T(z,t)=T_m(z,t)-T_c(z,t), ]

а уравнение роста необходимо решать как по длине трубы, так и во времени.

Это даёт локальную толщину (R_s(z,t)), что позволяет выявить области, где вероятны чрезмерное осаждение, закупорка или неравномерная кристаллизация.

Понимание компромиссов

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

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

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

Постоянные свойства могут быть нереалистичными

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

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

Состояние поверхности может изменять поведение границы раздела

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

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

Рост может снижать полезную производительность процесса

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

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

Выбор подхода в зависимости от цели

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

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

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

Сводная таблица:

Ключевой элемент	Инженерный вывод
Основное допущение	Тепловые профили устанавливаются быстрее, чем перемещается граница раздела твёрдой и жидкой фаз.
Движущая сила	Температура фазового перехода минус температура хладагента в объёме, ΔT = Tm − Tc.
Путь теплопередачи	Внутренняя плёнка хладагента, стенка трубы и растущий слой твёрдой фазы.
Баланс роста	Отведённое тепло равно скрытой теплоте, выделенной вновь образовавшейся твёрдой фазой.
Цилиндрический рост	Для толстых осадков используйте текущий радиус границы раздела и логарифмическое сопротивление твёрдой фазы.
Оценка для тонкого слоя	Сначала рост происходит почти линейно, а затем при увеличении сопротивления твёрдой фазы приближается к зависимости δ ∝ √t.
Основные факторы	Температура хладагента, скорость хладагента, теплофизические свойства, скрытая теплота и толщина осадка.
Ограничения модели	Нестационарный нагрев за счёт явной теплоты, переменные свойства, влияние состава, загрязнение и пористые осадки могут потребовать моделирования более высокого уровня точности.

Прогнозируйте и совершенствуйте свой трубчатый процесс

Нужно управлять замораживанием, кристаллизацией или осаждением твёрдой фазы в пилотной установке? Свяжитесь с LABPARK, чтобы обсудить вашу задачу с нашей командой.

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

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

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