Обтекание шара. Гидравлическая крупность

Зависимость коэффициента сопротивления шара от числа Рейнольдса имеет сложный вид (см. рис. 10.2). В первом приближении она может быть описана формулой [1]

 

или                                                        (8.12)

 

которая действительна при Re<105. В этой формуле Re=υd/ν (d - диаметр шара.) При очень малых числах Рейнольдса из уравнения (10.12) следует:

 

                                                   (8.13)

 

Подставляя это выражение в формулу (10.11), получим формулу Стокса:

 

Fд=3πμυd,                                                    (8.14)

 

При очень больших числах Рейнольдса

 

.                                                     (8.15)

 

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

 

,                                       (8.16)

где ρтв - плотность твердого тела;

  ρж - плотность жидкости;

  Сд - коэффициент сопротивления шара.

С учетом выражения (8.12) формула (8.16) принимает вид:

 

,                                      (8.17)

 

где ; ω – в см/с; d – в см; ν – в см2

Определение гидравлической крупности (скорости витания) весьма важно для расчетов гидро- и пневмотранспортирования, движения наносов и др.

 

 

Значения гидравлической крупности ω для частиц разного диаметра при их падении в неподвижной воде даны в табл. 8.2 (при температуре воды 20°С),

 

Примеры

  Пример 8.1. Плоская пластинка с размерами L =1 м и l =3 м (размер, перпендикулярный чертежу) и абсолютной эквивалентностью kэ=0,1 мм обдувается в ребро потоком воздуха со скоростью υ=50 м/с. Температура воздуха 15°С. Определить силу трения воздуха о пластинку.

Решение. Коэффициент сопротивления трения для турбулентного пограничного слоя определяем по формуле (8.7):

 

Сf=0,03(kэ/L+83/ReL)0,2

 

 

Кинематический коэффициент вязкости воздуха ν=1,45·10-5 м2/с (см. приложение 5).

Число Рейнольдса в рассматриваемом случае

 

.

 

Коэффициент сопротивления трения

 

.

 

Сила трения воздуха по двум сторонам пластинки при ρ= 1,2 кг/м3 (см. приложение 5)

Н.

 

Пример 8.2. Вычислить силу давления ветра, которую испытывает 1 м2 лобовой площади дымовой трубы (ω =1 м2). Коэффициент сопротивления такой трубы Сд=0,67 определен путем испытания модели. Наибольшая скорость ветра υ=50 м/с. Температура воздуха 15°С.

Решение. Плотность воздуха ρ= 1,21 иг/м3, давление ветра находим по формуле

Fд=Cд ωρν2/2=0,67·1·1,21·502/2≈1·103 Н

Пример 8.3. Осевая сила, с которой поток действует на круглую прямую трубу диаметром d=0,3 м, по динамометру R=7·102 Н (рис. 8.4). Определить

давление ρ1 на входе в трубу, если вода вытекает из трубы в атмосферу.

    Рис. 8.4                                                                                                    Рис. 8.5

 

Решение. Составим уравнение количества движения в проекции на направление движения для сечений 1—1 и 2—2:

 

ρQυ1 – ρQυ2 = p1ω – F

 

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

υ1= υ2; p1ω =R,

откуда

Пример 8.4. Определить избыточное давление на входе в диффузор с условием, чтобы сила, действующая на диффузор в направлении течения, равнялась нулю, если Q=0,01м3/с; d1=0,03 м; d2=0,1 м; α=60° (рис. 8,5).

Решение. Запишем уравнение количества движения в проекции на направление движения в виде

 

ρQυ1 – ρQυ2 = p1ω1 – p2ω2 + F

 

По условию задачи F=0. Выразим давление на выходе из диффузора через искомое давление p1, используя уравнение Бернулли:

 

P2=P1 + ρυ12/2 – ρυ22/2 –ΔPпот.

 

Найдем скорости на подходе к диффузору и на выходе из него:

 

υ1 =Q/ω1=1,27·0,01/0,032=14,1 м/с;

 

υ2 =Q/ω1=1,27·0,01/0,12 = 1,27 м/с.

 

Потери давления в диффузоре

 

ΔPпот=ζρυ22/2

 

Где

 

ζ=Кп.р.12-1)=kп.р.(d22/d12-1)

 

По табл. 4.3 при α=60° Кп.р.= 0,95. Тогда

 

ζ=0,95(0,12/0,032-1)2=95.

 

При плотности воды ρ=998,2 кг/м3, подставляя численные значения, получим:

 

ΔPпот=95·998,2·1,272/2=0,765·105 Па.

 

Тогда

 

.

 

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

 

 

и находим

P1= - 0,44·105 Па= - 44 кПа.

 

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

 Пример 8.5. Определить расход воздуха, поступающего в каждое отверстие квадратного сечения в промышленном здании (рис. 10.6). Вентиляция осуществляется за счет динамического воздействия ветра (ветрового давления). Скорость ветра υ=5 м/с; температура воздуха 20°С; площади отверстий ω1 =15 м2, ω2 =30 м2; ω3=10 м2.

Рис. 10.6

 

Решение. Ветровое давление на поверхность здания (на единицу площади) определяем по формуле (8.5).

 

 

Где С – в данном случае аэродинамический коэффициент сопротивления, зависящий от характера обтекания ветром рассматриваемой поверхности.

Ветровые коэффициенты принимаем соответственно: С1=0,5, С2= - 0,3,      С3= - 0,1. Плотность воздуха ρ=1,22 кг/м3.    

Давление ветра на наветренную поверхность у первого отверстия

 

 

Давление ветра на наветренную поверхность у первого отверстия

 

 

Давление у третьего отверстия

 

 

Предположим, что общий баланс воздуха в помещении имеет вид

 

 

где Q1 - расход воздуха, поступающего в помещение через первое отверстие; Q2 и

Q3 - расходы воздуха, уходящего из помещения через второе и третье отверстия соответственно.

Если давление в помещении обозначить через pп, то расходы воздуха в каждом отверстии (разность в подкоренном выражении всегда положительна):

 

 

 

 

 

 

Коэффициенты расхода отверстий μ1, μ2, μ3 в общем случае зависят от числа Рейнольдса. В первом приближении принимаем квадратичный законн истечения, тогда μ123. Следовательно,

 

 

 

Полученное трансцендентное уравнение решаем графически относительно давления в помещении рп. Для этого представим уравнение в виде

 

.

 

Задаваясь различными значениями рп, вычисляем соответствующие f(pп):

рп = 1

 

рп = 0

 

 

рп = - 2

 

Построим график зависимости f(рп) от рп (рис. 10.7). Решение уравнения находим при f(рп)=1, рп= - 1,8. Результаты расчета давления в помещении показывают, что так как рпвз, через третье отверстие воздух будет поступать в помещение.

          Рис. 8.7                                                            Рис. 8.8

Находим числа Рейнольдса и значения коэффициентов расхода μ при движении воздуха в отверстиях [ν=15·10-6 м2/с (см. приложение 4)]. Число Рейнольдса

 

 

Эквивалентный диаметр квадратного отверстия d=a=√ω находим по приложению 17. Тогда

 

.

 

По рис. 7.2 находим: μ1=0,60. Аналогичным способом находим: μ2=0,60;

Μ3=0,60.

Расход воздуха, поступающего в помещение:

через первое отверстие

=34,4 м3/с;

через третье отверстие

=4,3 м3/с.

Расход воздуха, уходящего из помещения через второе отверстие

=39 м3

Проверяем баланс воздуха в помещении:

 

 

Приток воздуха =38,7 м3/с; удаление воздуха =39 м3/с. Погрешность расчета составляет:

 

 

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

Пример 8.6. Струя, вытекающая из коноидального насадка диаметром

d=0,15 м, должна воздействовать на небольшую преграду с силой R=2·104 Н.

Определить расход воды Q и давление перед насадком p, если преграда делит струю на две части, отклоняемые на угол φ=60° (рис. 10.8).

Решение. Силовое воздействие струи в направлении ее оси определяем по

формуле (10.1):

 

F=ρων2(1-cosφ),

 

Откуда

 

Для коноидального насадка коэффициент сжатия струи ε= 1, поэтому площадь сжатого сечения струи ω равна площади выходного насадка ω0:

 

ω— ω0 = πd2/4=0,785·10-2 м2.

 

Принимая ρ = 998,2 кг/м3, находим необходимую скорость струи:

 

м/с

 

Расход, соответствующий этой скорости истечения,

 

=υω=48·1,77·10-2=0,85 м3

 

Расход связан с перепадом давления зависимостью

 

.

 

Принимая для коноидальното насадка μ0=0,98, имеем;

 

 

Таким образом, подавая к насадку расход воды  =0,85 м3/с под давле-

нием P = 1180 кПа, обеспечиваем необходимое динамическое воздействие на

преграду R=2·104 Н.

 


Понравилась статья? Добавь ее в закладку (CTRL+D) и не забудь поделиться с друзьями:  



double arrow
Сейчас читают про: