Фундаментальными управляющими уравнениями являются уравнение неразрывности для несжимаемой жидкости и уравнения движения Навье–Стокса. Для изотермической, несжимаемой ньютоновской жидкости с постоянными плотностью (\rho) и динамической вязкостью (\mu) они имеют вид

[ \nabla \cdot \mathbf{v}=0 ]

[ \rho\left(\frac{\partial \mathbf{v}}{\partial t} +\mathbf{v}\cdot\nabla\mathbf{v}\right)

-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}. ]

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

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

Система управляющих уравнений

Сохранение массы: уравнение неразрывности

Общее уравнение сохранения массы имеет вид

[ \frac{\partial \rho}{\partial t}+\nabla\cdot(\rho\mathbf{v})=0. ]

Для несжимаемой жидкости с постоянной плотностью

[ \boxed{\nabla\cdot\mathbf{v}=0} ]

Это условие означает, что элемент жидкости не расширяется и не сжимается при движении. В декартовых координатах

[ \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} =0, ]

где (u), (v) и (w) — компоненты скорости в направлениях (x), (y) и (z).

Сохранение импульса: уравнения Навье–Стокса

Баланс линейного импульса имеет вид

[ \boxed{ \rho\frac{D\mathbf{v}}{Dt}

-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g} } ]

где материальная производная равна

[ \frac{D\mathbf{v}}{Dt}

\frac{\partial\mathbf{v}}{\partial t} + (\mathbf{v}\cdot\nabla)\mathbf{v}. ]

Таким образом, уравнение импульса представляет собой баланс между:

Ускорением жидкости: (\rho D\mathbf{v}/Dt)
Силами давления: (-\nabla p)
Вязкими силами: (\mu\nabla^2\mathbf{v})
Объёмными силами: (\rho\mathbf{g}), обычно силой тяжести

Конвективный член ((\mathbf{v}\cdot\nabla)\mathbf{v}) имеет важное значение в изгибах, сужениях, расширениях, смесителях и других областях, где скорость изменяется в пространстве.

Компонентная форма в декартовых координатах

Для скорости (\mathbf{v}=(u,v,w)) уравнения импульса имеют вид

[ \rho\left( \frac{\partial u}{\partial t} +u\frac{\partial u}{\partial x} +v\frac{\partial u}{\partial y} +w\frac{\partial u}{\partial z} \right)

-\frac{\partial p}{\partial x} +\mu\nabla^2u+\rho g_x, ]

[ \rho\left( \frac{\partial v}{\partial t} +u\frac{\partial v}{\partial x} +v\frac{\partial v}{\partial y} +w\frac{\partial v}{\partial z} \right)

-\frac{\partial p}{\partial y} +\mu\nabla^2v+\rho g_y, ]

[ \rho\left( \frac{\partial w}{\partial t} +u\frac{\partial w}{\partial x} +v\frac{\partial w}{\partial y} +w\frac{\partial w}{\partial z} \right)

-\frac{\partial p}{\partial z} +\mu\nabla^2w+\rho g_z. ]

Эти три уравнения в сочетании с уравнением неразрывности дают четыре уравнения для четырёх основных полей: (u), (v), (w) и (p).

Почему ньютоновское поведение имеет значение

Закон вязкости Ньютона

В ньютоновской жидкости вязкие напряжения изменяются линейно со скоростью деформации. В тензорной форме

[ \boldsymbol{\tau}

2\mu\mathbf{D}, ]

для несжимаемого течения, где

[ \mathbf{D}

\frac{1}{2} \left[ \nabla\mathbf{v} + (\nabla\mathbf{v})^{T} \right] ]

— тензор скорости деформации.

Для постоянной вязкости и несжимаемого течения дивергенция вязкого напряжения сводится к

[ \nabla\cdot\boldsymbol{\tau}

\mu\nabla^2\mathbf{v}. ]

Именно это упрощение приводит к вязкому члену в стандартном уравнении Навье–Стокса.

Физическая интерпретация

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

Условия, необходимые для замыкания модели

Условие прилипания на стенке

На неподвижной твёрдой стенке

[ \mathbf{v}=0. ]

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

Для движущейся стенки скорость жидкости на стенке равна скорости стенки.

Граничные условия на входе

На входе может быть задан профиль скорости

[ \mathbf{v}=\mathbf{v}_{\text{in}}, ]

или массовый расход

[ \dot{m}

\int_A \rho,\mathbf{v}\cdot\mathbf{n},dA. ]

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

Граничные условия на выходе

Типичные условия на выходе включают заданное давление

[ p=p_{\text{out}}, ]

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

Начальное условие

Для нестационарного режима необходимо задать начальное поле скорости:

[ \mathbf{v}(\mathbf{x},0)=\mathbf{v}_0(\mathbf{x}). ]

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

Распространённые упрощения в модулях пилотных установок

Установившееся течение

Если рабочие условия существенно не изменяются со временем

[ \frac{\partial\mathbf{v}}{\partial t}=0. ]

Уравнение импульса принимает вид

[ \rho(\mathbf{v}\cdot\nabla)\mathbf{v}

-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}. ]

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

Полностью развитое внутреннее течение

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

Например, для одномерного течения между параллельными пластинами

[ 0=-\frac{dp}{dx}+\mu\frac{d^2u}{dy^2}. ]

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

Пренебрежимо малая инерция

При достаточно малом числе Рейнольдса инерционные эффекты малы по сравнению с вязкими. Уравнение импульса принимает форму уравнения течения Стокса:

[ 0=-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}. ]

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

Гравитация и гидростатическое давление

Если скорость пренебрежимо мала, уравнение импульса сводится к

[ \nabla p=\rho\mathbf{g}. ]

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

Применение уравнений к технологическим операциям

Трубопроводные сети и насосы

Уравнения описывают распределение скоростей и градиенты давления во всей трубопроводной системе. Затем для оценки потерь давления и определения требований к насосу используются интегральные или упрощённые формы.

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

Теплообменники

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

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

[ Q=\dot{m}C_p\Delta T, ]

а характеристики теплообменника часто выражаются соотношением

[ Q=UA\Delta T_{\mathrm{lm}}F. ]

Это дополнительные тепловые соотношения, а не замена уравнениям течения.

Реакторы и смесительное оборудование

В реакторах с мешалкой и проточных реакторах уравнения Навье–Стокса предсказывают циркуляцию, рециркуляцию и связанные со смешением поля скоростей.

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

Аппараты с насадкой и колонны

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

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

Двухжидкостные или многофазные модули

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

[ \mathbf{v}_1=\mathbf{v}_2 ]

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

[ \boldsymbol{\tau}_1\cdot\mathbf{n}

\boldsymbol{\tau}_2\cdot\mathbf{n} ]

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

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

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

Изотермичность не означает, что энергия всегда несущественна

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

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

Несжимаемость имеет определенную область применимости

Жидкости часто моделируются как несжимаемые, поскольку изменения плотности малы. Газы также можно рассматривать как несжимаемые, если скорость значительно ниже скорости звука, а вызванные давлением изменения плотности несущественны.

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

Постоянная вязкость может быть плохим приближением

Стандартное уравнение предполагает постоянную (\mu). Это разумно для многих жидкостей в умеренных диапазонах температур и скоростей сдвига, но не для жидкостей с сильной температурной зависимостью вязкости или неньютоновских материалов.

Если вязкость изменяется в пространстве, вязкий член, как правило, следует сохранять в форме дивергенции напряжений, а не записывать просто как (\mu\nabla^2\mathbf{v}).

Предположение о ламинарности не требуется уравнениями Навье–Стокса

Уравнения Навье–Стокса применимы к ньютоновскому течению независимо от того, является оно ламинарным или турбулентным. Однако прямое разрешение турбулентных масштабов требует значительных вычислительных ресурсов, поэтому для оборудования пилотного масштаба часто необходимы модели турбулентности или экспериментально подтверждённые корреляции.

«Вязкая диссипация» относится к балансу энергии

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

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

Правильный выбор в зависимости от цели

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

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

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

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

Элемент модели	Управляющее соотношение	Основное назначение
Сохранение массы	∇ · v = 0	Обеспечивает неразрывность несжимаемого течения
Сохранение импульса	ρ(∂v/∂t + v · ∇v) = −∇p + μ∇²v + ρg	Предсказывает эффекты скорости, давления, вязкости и объёмных сил
Ньютоновская вязкость	τ = 2μD	Связывает вязкое напряжение со скоростью деформации
Граничные условия	Условия на входе, выходе, стенке и начальные условия	Замыкает математическую модель
Дополнительная физика	Уравнения энергии, переноса компонентов, реакции или многофазного течения	Расширяет модель для описания теплопередачи, реакций или нескольких фаз

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