Сегодня я заварил себе чай и задумался
Сегодня утром я задумался, пока размешивал два кубика сахара в чашке с только что заваренным чаем. Задумался о форме жидкости, которую она принимает при вращении. Безусловно, все представляют себе что будет, если очень быстро начать размешивать сахар в чашке с чаем. Мне захотелось рассмотреть этот банальный и привычный процесс подробнее и попытаться рассказать Вам немного интересного из физики окружающих нас в быту явлений.
Идея эксперимента
Давайте представим, что мы имеем некоторую цилиндрическую тару, в которой находится некоторая жидкость. Вращаться жидкость можно заставить, как минимум, двумя очевидными способами: размешать её каким-нибудь предметом или начать вращать цилиндрическую тару, что, благодаря силам трения между жидкостью и поверхностью сосуда, приведет к вращению жидкости, увлекаемой содержащим её вращающимся сосудам.
Физическая модель
Остановимся на втором варианте. Итак, у нас есть вращающийся с постоянной циклической частотой сосуд, в котором при динамическом равновесии с постоянной циклической частотой вращается жидкость в том же направлении.
Вырежем из всей жидкости элементарный бесконечно малый объем около поверхности и рассмотрим какие силы на него действуют. В силу симметрии задачи, будем ориентироваться на цилиндрические координаты, что заметно упростит расчеты.
Качественный расчет формы поверхности
Запишем второй закон Ньютона для элементарного кусочка объема жидкости: ma = ΣF
К примеру, после размешивания ложкой сахара в чашке только что заваренного чая, жидкость вращается вокруг оси симметрии, отсюда наш элементарный кусочек объема имеет центростремительное ускорение. Поэтому спроецируем наш закон Ньютона на ось, совпадающую с радиусом-вектором от элементарного объема до оси симметрии. Не будем учитывать вязкость и поверхностное натяжение. Сила, сообщающая центростремительное ускорение (в правой части нашего закона движения) возникнет из-за разности давлений столбов жидкости, что можно увидеть на увеличенной части первого рисунка.
Таким образом, у нас получится следующее выражение:
Детальный расчет формы поверхности
Теперь мы получили вполне ясную зависимость для формы поверхности и с уверенностью можем сказать, что это параболоид. Но нам неизвестна постоянная величина C . Давайте её определим для полного понимания физики процесса.
Так как объем жидкости не меняется (мы считаем, что не пролили ни капли, пока размешивали наш чай ツ), то запишем объемы до вращения и во время вращения с постоянной циклической частотой.
До вращения:
Во время вращения:
Некоторые заметки
Хотелось бы обратить внимание на то, что форма поверхности зависит от частоты вращения, ускорения свободного падения, геометрических параметров сосуда, первоначального объема жидкости, но не зависит от плотности жидкости. Это выражение мне показалось довольно интересным, так как с его помощью можно легко смоделировать примерное расположение жидкости внутри вращающегося вокруг своей оси симметрии цилиндрического сосуда. Для этого можно воспользоваться MathCAD’ом и построить несколько графиков.
Графическое представление результатов расчета
Возьмем вполне реальные параметры системы, соизмеримые с размерами чашки или стакана.
Далее перепишем нашу функцию для её отображения в MathCAD.
Для 2D отображения сечения:
Для 3D отображения поверхности:
В качестве изменяющегося параметра будем менять циклическую частоту вращения w0. Результаты можно наблюдать на рисунках ниже:
Выводы
Видно, что если циклическая частота превысит значение 39.5 (рад/с), то мы увидим дно вращающегося цилиндрического сосуда, и, начиная с этой частоты, жидкость будет плавно «переходить» на стенки сосуда, всё сильнее оголяя дно. Очевидно, что при очень больших частотах вся жидкость растечется по стенкам сосуда. Теперь мы знаем все параметры такой жидкости. Зная её уравнение, не составит большого труда рассчитать толщину слоя жидкости на стенке сосуда на определенной высоте при определенной частоте.
Замечания к эксперименту
Отдельно хотелось бы подчеркнуть те противоречащие друг другу допущения, которые были приняты при рассмотрении задачи:
1. Считалось что, жидкость вращается благодаря вращению сосуда, который её содержит. Это может быть только при учете внутреннего трения, вязкости и поверхностного натяжения.
2. Но при выводе формы поверхности эти явления не учитываются для того, чтобы упростить решение и показать только качественный результаты моделирования. Т.е. решение немного противоречит описываемой изначально модели. Учет всех явлений, включая нелинейность процесса при высоких частотах, настолько бы усложнил задачу, что её вряд ли можно было бы решить аналитически и показать примерную и понятную модель для человека, который не связан с математикой/физикой.
3. Цель состоялась в том, чтобы показать лишь очень приближенное и самое простое решение, включающее в себя ряд допущений.
Эту статью я писал на habr примерно 3 года назад:
О форме вращающейся жидкости
И по-хорошему, надо бы усложнить эту модель с помощью рассмотрения уравнения Эйлера для движения жидкости