Скорость свободного движения плохообтекаемых тел в
жидкости
Д.т.н., профессор М.А. Михалев*,
ФГБОУ ВПО Санкт-Петербургский государственный политехнический университет
Ключевые слова: физическое моделирование; гидравлическая крупность; числа подобия; критерии подобия; критериальное уравнение
Явление осаждения (или всплывания) различных тел в жидкости (или газе) широко используется в технических устройствах (отстойники, камеры флотации и аэрации, гидро- и пневмотранспорт и т.п.). Скорость осаждения или всплывания тел в жидкости или газе, называемая гидравлической крупностью, представляется важной характеристикой явления. Однако существующие трудности математического моделирования позволяют получить решение задачи осаждения в жидкой или газообразной среде плохообтекаемых тел только при очень малых числах Рейнольдса (ползущее, или стоксово течение).
Что же касается методов физического моделирования явления, то оно в основном базируется на результатах моделирования обтекания закрепленных тел набегающем внешним потоком, согласно известному постулату относительности движения. Но, несмотря на кажущуюся аналогию движущегося и обтекаемого с той же скоростью закрепленного тела, с позиций методов физического моделирования - это два различных явления: в первом скорость движения заранее не известна, во втором - она изначально задана [1,2]. Сведения о свободном движении тел в жидкости, содержащиеся в литературных источниках, посвященных физическому моделированию явлений, в основном базируются на результатах изучения воздействия внешнего потока, набегающего на закрепленное тело [3 - 8]. Как правило, приводятся графики зависимости коэффициента сопротивления (числа Эйлера) от числа Рейнольдса. Но в задаче о гидравлической крупности ни число Рейнольдса, ни коэффициент сопротивления не являются критериями подобия, так как содержат скорость движения тела, которая до проведения опытов не известна. Иными словами, для свободно движущегося тела упомянутая зависимость представлена не в критериальной форме, поэтому воспользоваться ею можно только методом последовательных приближений.
При изучении сложных явлений, так или иначе связанных с гидравлической крупностью частиц, возникает необходимость правильного учета влияния на эту скорость (или на коэффициент сопротивления) дополнительных факторов [9 - 11]. Например, если тело, кроме сил тяжести и архимедовой, испытывает воздействие еще какой-то силы [9], то это приведет к изменению характера обтекания тела жидкостью, величины коэффициента сопротивления, в конечном итоге - скорости движения тела. Турбулентность внешнего потока и относительное содержание в нем частиц (несвободное движение тел в жидкости) могут оказать существенное влияние на величину гидравлической крупности [10]. Убедительные сведения на этот счет содержатся в [12]. Любые исследования на физических моделях (или основанные на методах физического моделирования явлений) должны заканчиваться обработкой экспериментальных данных обязательно в критериальной форме [11]. Ниже на примере гидравлической крупности будет показано, каким образом эти требования можно реализовать на практике.
Рассмотрим задачу о гидравлической крупности подробнее, основываясь на методах физического моделирования гидравлических явлений [1].
В основу положим принцип относительности движения, как это принято в механике: процесс обтекания тела одинаков, если закрепленное тело обтекается потоком с той же скоростью, с которой движется тело самостоятельно. Собственно, на этом принципе основывается физическое моделирование моделей летательных аппаратов в аэродинамических трубах. Будем учитывать пока только процессы во внешней среде, возникшие вследствие перемещения в ней тела (движение самого тела обязательно учтем ниже).
Отметим, что при стационарном равномерном движении тела в ньютоновской жидкости (несжимаемой или сжимаемой, но вдали от области, где сжимаемость среды начинает оказывать сопротивление движению - при малых числах Маха) явление характеризуется тремя числами подобия: Рейнольдса, Фруда и Эйлера. Среди них нет ни одного критерия, так как во все числа в качестве характерной входит заранее неизвестная скорость осаждения (или всплывания) ®0:
Яе ; Ег = р ; Ей=-Ар
V ' Р 8а (Р-Р)' р&о
Здесь d - характерный линейный размер тела, в качестве которого чаще всего принимается объем шара, равный объему тела; V и р - вязкость и плотность жидкости; р - плотность вещества, из которого состоит плохо обтекаемое тело; Ар - перепад давления вдоль контура тела, возникающий вследствие его обтекания внешним потоком жидкости; 8 - ускорение силы тяжести. Число подобия, обозначенное символом Егр , называют «плотностным числом Фруда», поскольку в него в качестве активной силы входит разность между силой тяжести и архимедовой силой. Эта разность определяет поведение тела в жидкости: если р >р , то тело тонет в жидкости, а если р <р , то оно всплывает.
Необходимо остановиться на величине перепада давления, в которую входят две составляющие. Первая, называемая сопротивлением формы тела, возникает вследствие особенности распределения давления по контуру тела: если тело имеет симметричную форму относительно продольной оси своего движения, то это разность между давлением в лобовой и кормовой областях тела. Вторая определяется силой трения, которую испытывает жидкость со стороны контура тела при его обтекании (отсюда и ее название - сопротивление трения). Этим объясняются большие трудности при определении перепада давления. На практике они преодолеваются путем замены перепада давления отношением силы сопротивления к площади приложения ее на контуре тела. Если тело закреплено, то при проведении опытов сила сопротивления легко определяется. Для свободно падающего тела силу сопротивления ¥Т можно вычислить, если задать ее формулу:
2 г"
¥„ = ; ^=СТ Р<°1 , (1)
где СТ - коэффициент сопротивления; Б - площадь приложения силы сопротивления. Запишем число Эйлера и примем Ар = ¥Т/Б в соответствии со второй формулой в (1), получим:
АР = СТ ра>1
Еи. (2)
Таким образом, если задать силу сопротивления формулой (1), то число Эйлера с точностью до постоянной станет равным коэффициенту сопротивления, входящему в виде сомножителя в эту формулу.
Воспользуемся далее условием равномерного движения тела в жидкости. Оно заранее известно, ради него не нужно ставить эксперимент, достаточно приравнять силу тяжести (с учетом архимедовой силы), приложенную к телу, силе сопротивления:
С 2
¥8 = ¥т; Г8(р1 -р)=-2-р®1 Б , (3)
где V - объем тела. Примем отношение У/Б равным d , а постоянную, в общем случае зависящую от формы тела, введем в коэффициент сопротивления. После преобразований из (3) найдем:
gd (р-рр 1 Ст
—=Еи . (4)
®0 Р 2 ' Егр
Михалев М.А. Скорость свободного движения плохообтекаемых тел в жидкости
Таким образом, из условия равномерного движения тела в жидкости получается, что число Эйлера,равное половине коэффициента сопротивления, обратно пропорционально плотностному числу Фруда. Следует отметить, что величину коэффициента сопротивления после определения гидравлической крупности в результате проведения экспериментов находят из левого уравнения, содержащегося в (4). Получается, что явление определяется не тремя, а только двумя числами подобия. Имея в виду дальнейшие преобразования, оставим для анализа числа подобия Рейнольдса и плотностное число Фруда. Используя метод комбинаций чисел подобия [1], найдем следующий критерий подобия:
Яе2 gd3(р -р)
— = 2 =АГ , (5) ^гр V р
который называется числом (в данном случае - критерием) Архимеда. В связи с тем, что критерий заменяет только одно число подобия, критериальное уравнение можно представить в двух видах:
Яе = /(Аг); Frр=^(Ar) . (6)
Приведем формулы, соответствующие первому и второму критериальному уравнению в (6), полученные в результате осаждения частиц кварца в воде. В области стоксова течения, которому соответствуют следующие условия: Яе < 1; Аг <16; Бгр < 0,0625; d < 0,01см , имеем:
_ Аг Аг
Яе=—; = —. (7)
18 р 182
В области квадратичного сопротивления, существование которой определяют условия: Яе> 100; Аг> 1,6 • 104; Бгр >0,625; d> 0,1см , получено:
Яе = 1,2Аг0,5; ?гр = 1,22 . (8)
В области, переходной от стоксова течения до квадратичного сопротивления, найдены следующие зависимости:
и Аг Аг
Яе =-0Т; Ргр=-лТТ . (9)
18 + 0,61Аг'5 Р (18 + 0,61Аг0'5)2
Помимо непосредственного использования рассмотренных зависимостей для определения величины гидравлической крупности, отметим еще один метод, с помощью которого на практике решается обратная задача. Он заключается в определении размеров мелких частиц в области стоксова движения после измерении в эксперименте их гидравлической крупности. Далее расчетом из соответствующей формулы определяется линейный размер частиц. Распространим этот метод и на другие области осаждения плохообтекаемых тел, имея в виду возможность его использования не только для определения размеров тел, но и для нахождения вязкости и плотности жидкости, а также плотности вещества, из которого состоят тела. Решим эту задачу, используя методы физического моделирования явлений.
Задача формулируется так, что в ней изначально не заданы ни скорость осаждения, ни линейный размер тела. Но после того как будет выполнен эксперимент, станет известна гидравлическая крупность. Однако среди чисел подобия критерии не появятся, поскольку характерный линейный размер тела не задан. Опуская часть приведенных выше обоснований вывода критериальных уравнений (о числе Эйлера, об использовании интегральных условий движения и т. п.), получим критерий подобия в рассматриваемой задаче. Он равен произведению числа Рейнольдса на число Фруда, так как в такую комбинацию чисел подобия не входит характерный линейный размер тела. Полученная комбинация называется числом (в данном случае это критерий) Келегана, обозначается Ке .
Заметим, что корень кубический из числа Келегана есть безразмерная гидравлическая крупность, которую обозначим символом Ш . Имеем:
КеЕгр = р
V(Р! -р)
=Ке; ^Кё=Ш=
Р
Vg (Р! -Р)
(10)
Отметим также, что корень кубический из числа Архимеда есть безразмерный линейный размер тела, который обозначим символом 1 :
3
Аг = ^ (Р -Р) ; ^ХТ = 1 =
V2 Р
g (Р1 -Р)
V1 Р
(11)
В связи с тем, что и число Рейнольдса, и число Фруда являются функциями числа Архимеда, то естественно, что их произведение (число Келегана) будет тоже функцией этого числа. Но, как об этом было сказано выше, число Келегана в поставленной задаче является критерием, поэтомукритериальное уравнение будет таким:
Аг = ^(Ке) .
(12)
В таком виде связь между числами подобия Архимеда и Келегана в переходной области выражается в виде уравнения четвертой степени, решение которого можно получить только методом последовательного приближения. В практических приложениях гораздо проще рассматривать обратную связь между этими числами. Иными словами, если связь между числами подобия вида (2) доказана, то справедлива и обратная зависимость:
Ке = ^(Аг); Ш = /(1) .
(13)
Обратимся к приведенным выше формулам с тем, чтобы получить в явном виде зависимости (13).
В области стоксова течения, в которой 1 <2,52; Ш <0,35 , получим, используя формулы (7):
Ке=ЯеБгр =
Аг2 _ 12
—г; а =— 183 18
(14)
В области квадратичного сопротивления, существующей при условиях 1 > 25,2; ш > 6 , найдем, обращаясь к формулам (8):
Ке=1,23Аг0'5; Ш=1,210'5
(15)
В области, переходной от стоксова течения до квадратичного сопротивления, получим, используя формулы (9):
Ке =
Аг
2
(18 + °,61Аг0,5)3 '
ш =-
12
18 + 0,611
1,5
(16)
В рамках сформулированной выше задачи приведенные формулы позволяют сравнительно легко по величине Ш , найденной в результате проведения опыта, определить величину 1 , а из нее расчетом - 1 . Как уже было отмечено ранее, исключение представляют расчеты по формулам (16), которые можно выполнить только методом последовательного приближения. Решение обратной задачи - поиск по величине 1 величины Ш - реализуется без преодоления
каких-либо трудностей. Следует только отметить, что по величинам Ш и 1 , предварительно найденным в результате выполненных опытов или расчетов, нужно установить область движения, границы которых указаны выше перед соответствующими формулами.
Михалев М.А. Скорость свободного движения плохообтекаемых тел в жидкости
Предложенная здесь методика расчетов скорости осаждения может быть использована для нахождения других параметров, определяющих явление. Например, правые формулы в зависимостях (7) и (9) позволяют по результатам опыта, в котором была установлена величина скорости осаждения и рассчитано число Фруда, найти число Архимеда, а из него - коэффициент кинематической вязкости среды (в области квадратичного сопротивления явление не зависит от вязкости среды, где происходит движение тела). Из левых формул зависимостей (7), (8) и (9) аналогичным образом находится число Рейнольдса, по его величине - область движения и число Архимеда, а из последнего уравнения - либо плотность вещества, из которого состоит тело, либо плотность среды.
Выводы
В заключение следует обратить внимание на обобщающий характер предложенных зависимостей, учитывающих влияние на процесс осаждения (или всплывания) плохообтекаемых тел основных сил: вязкого трения, тяжести, Архимеда. Если явления относятся к одному классу, то совсем не важно, какие в них вязкость и плотность среды, а также плотность вещества, из которого состоит тело, лишь бы были одинаковы числа подобия, тогда такие явления будут описываться одними и теми же зависимостями. К таким явлениям можно отнести: осаждение тяжелых частиц в водной и воздушной средах, капель воды в воздухе, всплывание пузырьков воздуха в воде и т. п.
Литература
1. Михалев М. А. Физическое моделирование гидравлических явлений: Учебное пособие / СПбГПУ. СПб., 2008. 443 с.
2. Михалев М. А. Теория подобия и размерностей: Учебное пособие. СПб.: СПбГТУ, 2001. 68 с.
3. Леви И. И. Моделирование гидравлических явлений. Л.: Энергия. 1967. 235 с.
4. Зегжда А. П. Теория подобия и методика расчета гидротехнических моделей. Л.-М.: Госстроойиздат. 1938. 162 с.
5. Шарп Дж. Гидравлическое моделирование. М.: Мир. 1984. 280 с.
6. Graf W. H. Hydraulics of Sediment Transport. New York: McGraw-Hill Book Company, 1970. 513 p.
7. Лойцянский Л. Г. Механика жидкости и газа. М.: Наука, 1973. 904 с.
8. Шлихтинг Г. Теория пограничного слоя. М.: Наука. 1969. 737 с.
9. Ватин Н. И., Гиргидов А. А., Стрелец К. И. Численное моделирование трехмерного поля скорости в циклоне // Инженерно-строительный журнал. 2011. №5. С. 5-9.
10. Маркелова Е. А. Технология применения активной примеси при строительстве и эксплуатации подводных трубопроводов // Инженерно-строительный журнал. 2011. №2. С. 55-62.
11. Ватин Н. И., Стрелец К. И, Китаин М. Б. Определение характеристик сварочных искр для расчета их удаления в циклоне // Инженерно-строительный журнал. 2011. №5. С. 25-30.
12. Соу С. Гидродинамика многофазных систем. М.: Мир, 1971. 536 с.
* Михаил Андреевич Михалев, Санкт-Петербург, Россия Тел. раб.: +7(812)535-46-10; эл. почта: mikhalev@cef.spbstu.ru