Аналитический метод заключается в решении задачи о поле течения Стокса вокруг вращающегося объекта, расчете поверхностного напряжения сдвига и интегрировании этого напряжения для получения крутящего момента. Для погруженной сферы радиуса (R), вращающейся с угловой скоростью (\omega) в ньютоновской жидкости с вязкостью (\mu), решение для низких чисел Рейнольдса дает момент сопротивления:

[ G=-8\pi\mu\omega R^3 ]

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

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

Почему применим метод течения Стокса

Пренебрежимо малая инерция жидкости

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

[ -\nabla p+\mu\nabla^2\mathbf{v}=0, \qquad \nabla\cdot\mathbf{v}=0 ]

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

Линейная зависимость между крутящим моментом и скоростью

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

[ |G|=8\pi\mu\omega R^3 ]

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

Решение поля вращательного течения

Применение граничных условий

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

[ v_r=0, \qquad v_\theta=0, \qquad v_\phi(R,\theta)=\omega R\sin\theta ]

Последнее условие представляет тангенциальную скорость поверхности сферы. Вдали от сферы жидкость неподвижна:

[ v_\phi(r,\theta)\rightarrow 0 \quad \text{при} \quad r\rightarrow\infty ]

Определение азимутальной скорости

Решение уравнений Стокса при этих условиях дает:

[ v_\phi(r,\theta)

\frac{\omega R^3\sin\theta}{r^2} ]

Скорость затухает как (1/r^2), что отражает локализованный характер вращательного возмущения при ползущем течении.

Интерпретация распределения течения

Скорость жидкости максимальна на экваторе сферы, где (\sin\theta=1), и обращается в ноль на оси вращения. Все поле течения создается за счет вязкого сцепления между вращающейся поверхностью и окружающей жидкостью.

Расчет поверхностного напряжения сдвига

Вычисление тангенциального напряжения

Для осесимметричного азимутального течения соответствующая компонента напряжения в сферических координатах равна:

[ \tau_{r\phi}

\mu \left( \frac{\partial v_\phi}{\partial r}

\frac{v_\phi}{r} \right) ]

Подстановка поля скоростей и вычисление при (r=R) дает:

[ \tau_{r\phi}(R,\theta)

-3\mu\omega\sin\theta ]

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

Учет плеча крутящего момента

Элемент поверхности под полярным углом (\theta) находится на перпендикулярном расстоянии от оси вращения, равном:

[ R\sin\theta ]

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

Интегрирование крутящего момента

Интегрирование по сфере

Элемент площади поверхности равен:

[ dA=R^2\sin\theta,d\theta,d\phi ]

Таким образом, полный крутящий момент равен:

[ G

\int_0^{2\pi} \int_0^\pi \tau_{r\phi}(R,\theta) (R\sin\theta) R^2\sin\theta ,d\theta,d\phi ]

Подстановка напряжения сдвига дает:

[ G

-3\mu\omega R^3 \int_0^{2\pi}d\phi \int_0^\pi\sin^3\theta,d\theta ]

Используя:

[ \int_0^{2\pi}d\phi=2\pi, \qquad \int_0^\pi\sin^3\theta,d\theta=\frac{4}{3} ]

получаем:

[ \boxed{G=-8\pi\mu\omega R^3} ]

Необходимый приложенный крутящий момент имеет ту же величину и противоположный знак:

[ \boxed{G_{\text{applied}}=8\pi\mu\omega R^3} ]

Связь крутящего момента с мощностью

Вязкая мощность, необходимая для поддержания вращения, равна:

[ P=|G|\omega ]

Для вращающейся сферы:

[ \boxed{P=8\pi\mu\omega^2R^3} ]

Это дает прямую оценку механической мощности, рассеиваемой в жидкости.

Расширение метода на типовые процессы

Использование коэффициентов вращательного сопротивления

Для других погруженных тел та же процедура выражается через коэффициент вращательного сопротивления:

[ G=-\zeta_r\omega ]

Для сферы:

[ \zeta_r=8\pi\mu R^3 ]

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

Применение к мешалкам

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

Различие локальных и суммарных величин

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

Понимание ограничений

Ограничение низкого числа Рейнольдса

Аналитический результат предполагает, что инерцией жидкости можно пренебречь. Если число Рейнольдса перестает быть малым, конвективная инерция изменяет поля скоростей и напряжений, и зависимость (G\propto\mu\omega R^3) может стать неточной.

Результат для сферы не универсален

Коэффициент (8\pi) применим конкретно к изолированной сфере, вращающейся в безграничной ньютоновской жидкости. Его нельзя напрямую переносить на импеллеры, цилиндры, частицы вблизи стенок или ограниченные объемы сосудов.

Предположение о ньютоновском поведении

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

Влияние границ может быть значительным

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

Правильный выбор для вашей цели

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

Если ваша основная задача — оценка крутящего момента для изолированной вращающейся сферы: используйте (G=-8\pi\mu\omega R^3), где знак представляет сопротивление приложенному вращению.
Если ваша основная задача — выбор размера низкоскоростной мешалки: определите специфическое для геометрии мешалки вращательное сопротивление, а не применяйте коэффициент сферы напрямую.
Если ваша основная задача — прогнозирование локального сдвига: используйте поле скоростей и распределение поверхностных напряжений перед интегрированием для получения полного крутящего момента.
Если ваша основная задача — работа вне режима ползущего течения: учитывайте инерцию жидкости и используйте соответствующие корреляции для более высоких чисел Рейнольдса, численные модели или эксперименты.

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

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

Ключевой аспект	Результат или применение
Управляющая модель	Линейные уравнения Стокса для ползущего течения
Крутящий момент сферы	`G = -8πμωR³`
Вращательное сопротивление	`ζᵣ = 8πμR³` для изолированной сферы
Требуемая мощность	`P = 8πμω²R³`
Метод расчета	Решить течение, вычислить поверхностное напряжение и проинтегрировать его момент
Основные ограничения	Форма, стенки, инерция и неньютоновские жидкости требуют поправок или других моделей

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