Mehanika fluida 1
Факултет
техничких
наука
Нови
Сад
МЕХАНИК
A
ФЛУИДА
1
Маша
Букуров
октобар
, 2011.
1
Кинематика
флуида
и
напонско
стање
У
механици
флуида
решавају
се
динамичке
једначине
кретања
за
разна
струјања
која
су
присутна
у
техничким
,
физичким
,
биолошким
и
другим
областима
.
При
томе
упознају
се
све
спољашње
и
унутрашње
карактеристике
,
односно
инерцијске
,
вискозне
,
притисне
,
турбулентне
,
спољашње
,
еластичне
и
остале
силе
(
величина
,
карактер
,
повезаност
,
дејства
,
утицај
и
последице
).
Кинематика
флуида
,
која
се
бави
проучавањем
геометријских
и
кинематских
карактеристика
струјања
,
само
привидно
је
без
нарочитог
значаја
у
решавању
струјних
проблема
.
Њен
значај
лежи
у
чињеници
да
је
до
данас
решен
ограничен
број
стварних
струјних
проблема
(
реалан
флуид
),
па
је
свака
информација
о
струјању
,
без
обзира
са
које
стране
долази
(
нпр
.
аналогија
са
оптиком
,
електротехником
,
магнетизмом
,
еластомехаником
),
драгоцена
.
Струјна
слика
са
распоредом
брзина
формира
се
непосредно
под
дејством
наведених
динамичких
величина
.
Флуид
је
непрекидна
средина
,
врло
покретна
и
деформабилна
.
Непрекидност
се
задржава
и
при
изузетно
високим
деформацијама
које
су
често
присутне
у
струјном
пољу
.
Уместо
деформација
разматрају
се
умереније
изражене
величине
–
брзине
деформација
које
су
повезане
са
напонским
стањем
.
За
покретан
флуид
потребно
је
да
се
одреде
следеће
величине
:
( ,
,
),
,
,
x
y
z
v
v v v
p
T
ρ
=
Свеукупност
ових
величина
у
посматраном
простору
и
времену
описује
струјно
поље
.
Струјно
поље
је
стационарно
,
ако
су
све
горе
наведене
величине
само
функције
положаја
.
Када
су
ове
величине
променљиве
функције
и
времена
,
поље
је
нестационарно
.
Постоје
два
начина
за
проучавање
кретања
флуида
:
o
Лагранжов
,
када
се
за
сваки
флуидни
делић
одређује
путања
,
брзина
и
убрзање
;
o
Ојлеров
,
који
се
састоји
у
праћењу
брзина
,
убрзања
и
других
својстава
флуидних
делића
у
једној
непокретној
тачки
поља
.

3
O
во
описивање
одговарајуће
je
за
величине
које
су
везане
за
одређене
флуидне
делиће
,
нпр
.
распростирање
полутаната
.
Сви
закони
одржања
(
масе
,
количине
кретања
и
енергије
)
такође
се
физички
јасније
представљају
Лагранжовим
методом
.
Међутим
,
Лагранжова
метода
захтева
велики
број
једначина
за
једно
струјно
поље
,
због
великог
броја
флуидних
делића
,
те
је
стога
неприхватљива
за
већину
проблема
.
Такође
,
било
каква
мерења
тешко
су
изводљива
,
јер
би
мерни
апарати
морали
да
прате
флуидне
делиће
.
Ојлеров
метод
Овим
методом
разматра
се
промена
струјних
величина
у
једном
сталном
месту
простора
,
док
делићи
пролазе
кроз
ово
место
.
Струјне
величине
могу
да
се
мере
непокретним
мерним
уређајима
.
Међусобна
веза
између
Лагранжовог
и
Ојлеровог
метода
дата
је
у
следећој
једначини
за
својство
флуидног
делића
ϕ
:
d
d
d
d
d
d
d
d
grad
x
y
z
x
y
z
t
t
x
t
y
t
z
t
v
v
v
v
t
x
y
z
t
ϕ
ϕ
ϕ
ϕ
ϕ
ϕ
ϕ
ϕ
ϕ
ϕ
ϕ
∂
∂
∂
∂
=
+
+
+
=
∂
∂
∂
∂
∂
∂
∂
∂
∂
=
+
+
+
=
+
∂
∂
∂
∂
∂
Материјална
(
супстанцијална
)
промена
величине
ϕ
(
тотални
извод
)
(
која
може
да
се
одреди
Лагранжовом
методом
)
једнака
је
збиру
њене
локалне
и
конвективне
промене
(
Ојлеров
метод
).
Локална
промена
представља
промену
својства
ϕ
у
једној
тачки
флуидног
простора
током
времена
,
а
конвективна
промена
је
промена
својства
флуида
због
тога
што
делић
флуида
прелази
из
једне
у
другу
тачку
струјног
простора
,
тј
.
конвективна
промена
описује
како
поље
брзина
утиче
на
величину
ϕ
(
grad
v
ϕ
).
На
пример
,
својство
ϕ
може
да
буде
температура
Т
.
Тада
тотални
извод
температуре
може
да
се
прикаже
у
виду
локалне
временске
промене
температуре
и
конвективне
промене
која
се
односи
на
утицај
поља
брзина
на
ту
промену
.
d
grad
d
T
T
v
T
t
t
∂
=
+
∂
Кинематика
флуида
обрадиће
се
кроз
упознавање
са
битним
струјним
својствима
у
скаларном
,
векторском
и
тензорском
пољу
.
4
Скаларно
поље
Скаларно
поље
је
простор
у
коме
је
дефинисана
скаларна
функција
положаја
U
=
U
(
x,y,z
).
Скаларна
функција
положаја
може
да
описује
простор
,
површину
или
линију
,
па
су
скаларна
поља
тада
просторна
,
површинска
или
линијска
.
Линије
или
површине
на
којима
скаларна
функција
задржава
исту
вредност
су
еквискаларне
.
Градијент
скаларне
функције
Слика
2.
Градијент
скаларне
функције
–
правац
најоштрије
промене
Ради
описивања
правца
промене
скаларне
величине
,
својства
општег
за
све
врсте
скаларних
поља
,
уводи
се
диференцијално
-
векторски
оператор
набла
i
j
k
x
y
z
∂
∂
∂
∇ =
+
+
∂
∂
∂
који
се
примењује
на
скаларну
функцију
U=U
(
x,y,z
)
и
одређује
градијент
скаларне
функције
U
(grad
U
):
grad
U
U
U
U
i
j
k
x
y
z
∂
∂
∂
=
+
+
∂
∂
∂
.
grad
U
(
вектор
)
има
правац
и
смер
нормале
у
произвољној
тачки
површине
U=U
(
x,y,z
)
,
а
пошто
је
нормала
најкраће
растојање
између
две
вредности
скалара
U
1
и
U
2
,
то
grad
U
представља
правац
најјаче
промене
скалара
U
.
Смер
grad
U
је
смер
пораста
промене
скалара
(
слика
2
).
За
операције
са
grad
важе
иста
правила
као
и
за
диференцирање
.
Са
познатим
grad
U
може
да
се
одреди
промена
скаларне
функције
U
у
било
ком
правцу
r
.
Тотални
диференцијал
(
)
0
d
d
d
d
grad ,
d
U
U
U
U
x
y
z
U r
r
x
y
z
∂
∂
∂
=
+
+
=
∂
∂
∂

6
Путања
представља
"
траг
"
који
за
собом
оставља
један
флуидни
делић
.
При
стационарном
кретању
струјнице
су
уједно
и
путање
.
Једначина
струјнице
одређена
је
векторском
једначином
:
[ , d ]
0
v
s
=
Развијањем
векторског
производа
добија
се
:
d
d
d
x
y
z
x
y
z
v
v
v
=
=
Увођење
струјнице
одговара
Ојлеровом
начину
посматрања
.
Струјна
цев
Слика
4.
Струјна
цев
Без
неких
значајнијих
промена
карактеристика
струјнице
,
могуће
је
проширити
струјницу
на
струјну
цев
(
слика
4
),
тј
.
на
скуп
свих
струјница
које
пролазе
кроз
све
тачке
затворене
криве
линије
.
У
проточном
пресеку
струјне
цеви
брзине
не
морају
бити
једнаке
,
док
се
притисци
сматрају
непроменљивим
,
кад
год
су
изводнице
струјне
цеви
приближно
паралелне
.
Дивергенција
вектора
брзине
Примена
оператора
∇
преко
скаларног
производа
на
векторску
функцију
доводи
до
скаларне
величине
–
дивергенције
те
векторске
функције
:
(
)
div
y
x
z
x
y
z
v
v
v
v
v
i
j
k
v i
v j
v k
x
y
z
x
y
z
∂
∂
∂
∂
∂
∂
∇ ⋅ =
=
+
+
⋅
+
+
=
+
+
∂
∂
∂
∂
∂
∂
7
Проток
Дивергенција
брзине
може
да
се
доведе
у
везу
са
запреминским
протоком
кроз
произвољну
флуидну
запремину
V
омеђену
површином
А
(
слика
5
).
За
стационарно
струјање
,
ако
тачкама
елементарног
проточног
пресека
припадају
брзине
v
,
мора
бити
:
d
d
( , d )
m
Q
v
A
ρ
ρ
=
=
тј
.
( ,d )
A
m
v
A
ρ
=
∫
где
су
:
m
-
масени
проток
флуида
кроз
површину
А
[kg/s]
Q
-
запремински
проток
флуида
кроз
површину
А
[m
3
/s].
За
течности
густина
се
не
мења
,
па
је
1
2
0
1
1
2
2
0
0
( ,d )
d
d
d
A
A
A
A
m
Q
v
A
v A
v A
v A
ρ
=
=
= −
+
±
∫
∫
∫
∫
где
су
:
A
-
површина
–
омотач
запремине
флуида
V
:
0
1
2
A
A
A
A
=
+
+
v
-
тј
.
0
1
2
,
,
v
v v
-
брзине
нормалне
на
одговарајуће
делове
површине
А
.
Слика
5.
div
v
као
мера
запреминског
протока
кроз
површину
А
Површински
интеграл
може
да
се
замени
запреминским
,
те
следи
div d
V
Q
v V
=
∫
Ако
већа
количина
флуида
истиче
из
запремине
V
него
што
утиче
,
Q
је
веће
од
нуле
,
односно
div
0
v
>
,
тј
.
из
посматране
запремине
извире
течност
која
се
додаје
струји
.
Дивергенција
брзине
div
v
представља
,
дакле
,
разлику
количине
флуида
која
истече
из
јединице
запремине

9
0
y
x
z
x
y
z
v
v
v
v
v
v
t
x
y
z
x
y
z
ρ
ρ
ρ
ρ
ρ
∂
∂
∂
∂
∂
∂
∂
+
+
+
+
+
+
=
∂
∂
∂
∂
∂
∂
∂
или
(
)
grad
div
0
v
v
t
ρ
ρ
ρ
∂
+
+
=
∂
или
div (
)
0
v
t
ρ
ρ
∂
+
=
∂
.
Када
је
струјање
стационарно
,
а
флуид
нестишљив
,
једначина
континуитета
представља
се
са
div
0
v
=
односно
d
0
Q
=
.
Ако
се
у
флуидној
запремини
налази
извор
или
понор
издашности
±
ε
,
једначина
континуитета
је
div (
)
v
t
ρ
ρ
ρε
∂
+
=
∂
где
је
d
V
ε
ε
±
=
-
специфична
издашност
извора
или
понора
[m
3
/s/m
3
].
Једначина
континуитета
представља
закон
одржања
масе
флуида
и
уз
карактеристичну
и
динамичку
једначину
(
законе
кретања
)
служи
за
одређивање
непознатих
струјних
параметара
:
ρ
,
p
и
v
.
Ротор
брзине
Ако
се
оператор
набла
∇
примени
на
векторску
функцију
брзине
преко
векторског
производа
,
добија
се
сложена
векторско
-
диференцијална
величина
која
се
назива
ротор
брзине
rot
v
.
Математички
израз
у
Декартовим
координатама
је
10
[
]
rot
,
y
y
x
x
z
z
x
y
z
i
j
k
v
v
v
v
v
v
v
v
i
j
k
x
y
z
y
z
z
x
x
y
v
v
v
∂
∂
∂
∂
∂
∂
∂
∂
∂
= ∇
=
=
−
+
−
+
−
∂
∂
∂
∂
∂
∂
∂
∂
∂
У
физичком
смислу
ротор
брзине
представља
меру
обртања
флуидног
делића
око
сопствене
тежишне
осе
.
Вртлог
Флуидни
делићи
које
пролазе
кроз
сужење
подлежу
и
деформацији
и
ротацији
као
што
се
може
видети
на
слици
6
.
Слика
6.
Деформација
и
ротација
флуидних
делића
који
пролазе
кроз
сужење
Претпоставиће
се
,
као
што
је
приказано
на
слици
7
.
да
се
елементарни
флуидни
делић
,
приказан
у
равни
у
виду
правоугаоника
ABCD
са
странцама
d
x
и
d
y
у
тренутку
t
налази
се
са
својим
центром
у
тачки
О
.
У
тренутку
t
+d
t
овај
флуидни
делић
заузима
ново
место
у
простору
О
’
деформишући
се
при
томе
у
четвороугао
A’B’C’D’.
Страница
AB
крећући
се
у
правцу
x
-
осе
и
ротирајући
се
за
d
ε
1
заузима
ново
место
у
простору
A’B’.
Страница
AD
креће
се
у
y
правцу
и
ротира
се
за
d
ε
2
.
1
2
1
2
1
2
d
d d
d
d d
d
d
d
d
d
d
d
d
y
x
y
x
v
v
x t
y t
x
y
v
v
t
t
x
x
y
y
ε
ε
ε
ε
θ
θ
∂
∂
=
= −
∂
∂
∂
∂
=
=
=
= −
∂
∂
Угаоне
брзине
страница
AB
и
AD
су
ω
1
и
ω
2
:
1
2
1
2
d
d
d
d
y
x
v
v
t
x
t
y
θ
θ
ω
ω
∂
∂
=
=
=
= −
∂
∂

12
rot
2
v
ω
=
назива
се
вртлогом
.
Вртлог
је
позитиван
уколико
је
обртање
флуидних
делића
у
смеру
супротном
од
кретања
казаљке
на
часовнику
.
Ротор
брзине
показује
колико
је
и
какво
обртање
флуидног
делића
у
ма
којој
тачки
струјног
поља
око
сопствене
тежишне
осе
.
Слика
8.
Вртложно
струјање
;
а
)
заробљен
–
вештачки
вртлог
б
)
слободан
вртлог
Као
што
је
приказано
на
слици
8
,
цилиндричан
суд
који
садржи
течност
обрће
се
око
вертикалне
осе
одређеном
угаоном
брзином
.
Када
се
сва
течност
обрће
константном
угаоном
брзином
ω
=
const.
око
неке
осе
(
тачке
),
као
да
се
ради
о
чврстом
телу
,
а
не
флуиду
;
сви
њени
делићи
обрћу
се
истом
угаоном
брзином
око
те
исте
осе
(
тачке
).
Ово
је
тзв
.
заробљен
–
вештачки
вртлог
(
слично
релативном
мировању
у
статици
флуида
).
На
слици
лево
приказан
је
случај
тзв
.
заробљеног
вртлога
,
при
коме
дрвце
које
плива
на
слободној
површини
ротира
заједно
са
судом
.
Слика
десно
приказује
слободан
вртлог
,
када
дрвце
не
ротира
.
Овај
случај
среће
се
кад
год
флуид
истиче
кроз
мали
отвор
на
дну
суда
.
Иако
се
флуид
обрће
,
његови
делићи
заузимају
исти
правац
без
ротирања
.
Урагани
,
речни
вирови
,
торнада
примери
су
природних
вртлога
.
Иако
је
структура
ових
вртлога
комплексна
,
основни
елементи
који
чине
ове
вртлоге
су
заробљен
вртлог
у
центру
и
слободан
вртлога
на
периферији
.
Многи
од
вртлога
у
природи
су
овог
типа
.
13
Циркулација
Линијски
интеграл
по
затвореној
кривој
С
(
слика
9
)
скаларног
производа
вектора
брзине
v
и
вектора
елемента
затворене
криве
линије
d
s
означава
се
са
Г
(
гама
)
и
назива
се
циркулација
.
( ,d )
cos d
(rot ,d )
(2 ,d )
C
C
A
A
v
s
v
s
v
A
A
α
ω
Γ =
=
=
=
∫
∫
∫
∫
.
Коначна
вредност
циркулације
означава
присуство
"
протока
вртлога
2( , d )
A
A
ω
∑
"
кроз
целу
површину
А
.
При
томе
,
исто
као
и
за
вртлог
,
уколико
се
при
обилажењу
контуре
С
у
правцу
супротном
од
кретања
казаљке
на
сату
добија
вредност
већа
од
нуле
,
циркулација
је
позитивна
и
обрнуто
.
Слика
9.
Циркулација
Г
Слика
10.
Независност
Г
од
облика
и
величине
затворене
контуре
Ако
се
само
у
једној
тачки
В
струјног
поља
(
слика
10
)
флуидни
делићи

15
d
d
d
x
y
z
x
y
z
ω
ω
ω
=
=
За
вртложна
влакна
и
цеви
важи
закон
о
неуништивости
вртлога
и
представља
се
једначином
(
за
.
const
ω =
по
пресеку
А
)
1 1
2
2
A
A
ω
ω
=
где
су
површине
А
1
и
А
2
нормалне
на
ω
1
и
ω
2
.
Слика
11.
Вртложно
поље
Слика
12.
Неуништивост
идеалног
вртлога
Вртложне
линије
не
почињу
и
не
завршавају
се
у
вртложном
пољу
.
Неуништивост
вртлога
важи
за
идеализовани
случај
вртложења
без
трења
(
идеалан
флуид
–
слика
12
).
У
реалном
флуиду
вртлози
,
који
се
најчешће
принудно
формирају
,
нестају
током
времена
.
Њихова
енергија
претвара
се
у
топлотну
и
предаје
околном
флуиду
.
Неуништивост
вртлога
објашњава
дуго
време
постојања
циклон
,
торнада
,
урагана
,
одн
.
вртложних
формација
у
атмосфери
.
Сложено
поље
.
Код
сложеног
поља
важи
div
0
v
≠
и
rot
0
v
≠
.
У
оваквом
пољу
има
и
вртлога
и
извора
и
понора
.
Компликовано
је
за
проучавање
.
16
Векторско
поље
убрзања
Поље
убрзања
значајно
је
јер
репрезентује
инерцијске
силе
.
У
свакој
тачки
овог
поља
одређен
је
вектор
убрзања
флуидних
делића
који
се
у
датим
тренуцима
налазе
у
њој
.
Инерцијске
силе
,
као
и
остале
величине
у
механици
флуида
,
изражавају
се
по
јединици
масе
флуида
,
те
су
представљене
убрзањем
и
имају
такву
димензију
.
d
( , )
d
x
y
z
v
v
v
v
v
v
a
v
v
v
v
v
t
t
x
y
z
t
∂
∂
∂
∂
∂
=
=
+
+
+
=
+
∇
∂
∂
∂
∂
∂
где
су
:
d
d
v
t
-
супстанцијални
или
тотални
извод
брзине
по
времену
;
v
t
∂
∂
-
локални
извод
по
времену
,
одн
.
како
се
мења
брзина
v
у
уоченој
тачки
струјног
поља
са
временом
:
( , )
v
v
∇
-
назива
се
конвективни
извод
брзине
v
и
представља
промену
брзине
v
у
једној
тачки
поља
при
њеном
померању
у
неком
правцу
.
Тензорско
поље
У
савременој
физици
проучавају
се
појмови
који
не
могу
да
се
опишу
уобичајеним
изразима
више
математике
.
У
ову
групу
спадају
и
тензори
.
Тензори
су
скупови
више
величина
којима
се
дефинише
неко
физичко
својство
тела
у
простору
.
Простор
где
су
тензори
дефинисани
назива
се
тензорским
пољем
.
У
таквим
пољима
као
и
у
скаларним
и
векторским
,
тензори
су
функција
положаја
.
У
наставку
помињу
се
тензори
брзине
деформисања
флуидних
делића
и
напони
који
су
с
тим
повезани
.
Тензор
брзина
деформисања
делића
Потпуна
слика
о
кретању
делића
течности
добија
се
једино
када
се
,
поред
кретања
,
познаје
и
начин
на
који
се
делићи
деформишу
.
С
обзиром
на
то
што
гасови
и
течности
представљају
веома
покретљиве
средине
,
делићи
флуида
се
јако
деформишу
,
тако
да
је
важнија
брзина
деформисања
од
самих
деформација
.
Овом
одликом
флуиди
се
битно
разликују
од
чврстих
тела
где
су
деформације
веома
мале
и
споре
.

18
div
y
x
z
v
v
v
v
x
y
z
∂
∂
∂
=
+
+
∂
∂
∂
што
представља
брзину
релативне
запреминске
промене
делића
(
слика
14
),
јер
се
под
дејством
нормалних
напона
првобитна
запремина
шири
или
скупља
.
Слика
14.
div
v
-
брзина
релативне
запреминске
промене
делића
Из
раније
поменуте
једначине
div d
V
Q
v V
=
∫
div
v
представља
и
елементарни
проток
флуида
.
Тензор
напона
–
напонско
стање
Ако
се
из
бестежинског
вискозно
-
струјног
флуидног
простора
издвоји
флуидни
елемент
,
његово
равнотежно
стање
дато
је
одређеном
конфигурацијом
тангенцијалних
и
нормалних
напона
,
који
представљају
дејство
околне
вискозне
средине
на
издвојени
флуидни
делић
.
Нормални
напони
усмерени
су
у
правцу
спољашњих
нормала
граничних
површина
издвојеног
делића
и
представљају
дејство
привлачних
–
вискозних
нормалних
сила
околине
.
Тангенцијални
напони
распоређени
су
у
равнима
–
омотачима
издвојеног
делића
,
а
представљају
дејствовање
смичућих
сила
околине
због
издвајања
њеног
дела
.
Први
индекс
означава
осу
која
је
нормална
на
површиницу
,
а
други
правац
осе
по
којој
се
напон
испољава
.
Напони
σ
xx
,
σ
yy
и
σ
zz
су
19
нормални
напони
услед
чисте
вискозности
,
а
τ
yx
=
τ
xy
,
τ
zx
=
τ
xz
и
τ
zy
=
τ
yz
су
тангенцијални
напони
.
Нормални
и
тангенцијални
напони
имају
позитиван
знак
кад
смер
нормале
на
површину
и
смер
деловања
напона
имају
исти
знак
у
односу
на
оријентацију
координата
.
У
противном
,
напони
су
негативни
(
слика
15
).
Слика
15.
Распоред
напона
на
издвојеном
флуидном
делић
у
Тензор
напона
вискозних
сила
представљен
је
са
:
xx
yx
zx
xy
yy
zy
xz
yz
zz
S
σ
τ
τ
τ
σ
τ
τ
τ
σ
≡
.
Веза
напона
са
брзинама
деформисања
Њутн
је
успоставио
везу
између
напонског
и
деформацијског
стања
у
виду
следећих
једначина

21
Веза
између
напона
и
брзина
деформисања
користи
се
за
извођење
Навије
-
Стоксове
једначине
.
Притисак
Притисак
у
једној
тачки
струјног
простора
има
константну
вредност
,
без
обзира
на
оријентацију
површине
на
коју
дејствује
.
Он
је
последица
дејства
спољашњих
поремећаја
(
клип
који
притискује
флуид
)
или
запреминских
спољашњих
сила
(
привлачна
сила
Земљине
теже
–
тежина
флуида
).
Притисак
дејствује
у
истом
правцу
као
и
нормални
–
вискозни
напон
,
али
је
супротно
оријентисан
–
усмерен
ка
омотачу
издвојене
флуидне
запремине
.
Разлика
интензитета
притиска
и
нормалног
–
вискозног
напона
одређује
општи
нормални
напон
,
који
је
усмерен
ка
површини
на
којој
се
мери
(
не
важи
правило
са
индексима
).
Три
нормална
(
општа
)
напона
дата
су
са
:
;
;
xx
xx
yy
yy
zz
zz
p
p
p
p
p
p
σ
σ
σ
= −
= −
= −
Сабирање
ова
три
нормална
напона
,
водећи
рачуна
да
је
0
xx
yy
zz
σ
σ
σ
+
+
=
добија
се
(
)
1
3
xx
yy
zz
p
p
p
p
=
+
+
Ова
једначина
служи
за
одређивање
вредности
притиска
у
једној
тачки
струјног
простора
.
Притисак
је
непроменљив
у
једној
тачки
и
одређује
се
преко
нормалних
напона
p
xx
,
p
yy
,
p
zz
који
се
мере
појединачно
.
За
флуид
,
у
стању
мировања
,
важи
0
xx
yy
zz
σ
σ
σ
=
=
=
,
односно
xx
yy
zz
p
p
p
p
=
=
=
.
У
динамици
флуида
статички
притисак
одговара
горњој
релацији
.
22
Потенцијално
струјање
M
одел
потенцијалног
струјања
са
граничним
слојем
Теорија
потенцијалног
струјања
најстарија
је
област
аналитичке
механике
флуида
и
датира
из
половине
XVIII
века
.
Све
до
Прантлове
хипотезе
математичка
решења
бројних
случајева
потенцијалног
струјања
нису
имала
(
или
су
имала
врло
ограничену
)
инжењерску
примену
.
Стварно
струјно
поље
аналитички
се
представља
комбинацијом
потенцијалног
поља
и
граничног
слоја
.
Коришћењем
модела
потенцијалног
струјања
са
граничним
слојем
добија
се
могућност
аналитичког
решавања
унутрашњих
и
спољашних
струјања
(
губици
у
цевима
,
силе
отпора
и
узгона
при
опструјавању
тела
).
Значај
рачунског
решавања
стварних
струјних
проблема
огледа
се
у
томе
што
нема
потребе
за
експерименталним
одређивањем
карактеристичних
,
често
најважнијих
елемената
струјања
као
што
су
коефицијенти
:
локалних
губитака
,
узгона
и
отпора
.
Постављање
модела
потенцијалног
струјања
са
граничним
слојем
уследило
је
због
знатних
одступања
која
су
се
јављала
при
аналитичком
решавању
струјања
са
високим
Рејнолдсовим
бројем
,
уз
претпоставку
о
занемареној
сили
трења
(
R
е
је
однос
инерцијских
и
вискозних
сила
).
Примери
занемаривања
вискозних
и
инерцијских
сила
су
Бернулијева
једначина
за
идеалан
флуид
и
потпуно
развијено
струјање
у
цевима
.
Прантл
је
1904.
године
први
поставио
хипотезу
граничног
слоја
по
којој
су
,
за
струјања
са
великим
Рејнолдсовим
бројевима
,
ефекти
вискозности
и
вртложења
ограничени
на
узан
слој
флуида
уз
граничну
површину
чврстог
тела
или
површину
дисконтинуитета
(
гранична
површина
млаза
).
Изван
ове
области
струјање
се
понаша
као
идеално
:
невискозно
и
невртложно
,
тј
.
потенцијално
.
На
слици
1
приказане
су
области
граничног
слоја
и
потенцијалног
струјања
за
пример
спољашњег
и
унутрашњег
струјања
.

24
струјна
слика
спољашњег
и
потенцијалног
струјања
потенцијалног
струјања
из
претходно
одређеног
распореда
притиска
.
Повезаност
решавања
потенцијалног
струјања
и
граничног
слоја
је
приказана
на
слици
2.
Повратна
спрега
(
корекције
)
прорачуна
граничног
слоја
и
потенцијалног
струјања
служи
за
добијање
тачнијег
распореда
притиска
што
је
посебно
значајно
при
решавању
проблема
унутрашњих
струјања
као
што
је
нпр
.
улазна
област
цеви
.
Сл
.2.
Шематски
приказ
прорачуна
при
великим
R
е
модела
потенцијалног
струјања
и
граничног
слоја
Иако
се
у
новије
време
мерењем
распореда
притиска
на
граници
и
кроз
гранични
слој
у
потпуности
избегава
решавање
потенцијалног
струјања
,
потреба
за
коришћењем
модела
потенцијалног
струјања
са
граничним
слојем
остаје
,
јер
представља
компоненту
теоријског
развоја
распоред
тангенцијалног
напона
на
чврстој
површини
,
τ
0
површински
распоред
притиска
притисак
корекције
(
према
избору
,
нпр
.
узимање
у
обзир
δ
за
смањењење
потенцијалног
струјног
тока
)
дебљина
граничног
слоја
прорачун
потенцијалног
струјања
прорачун
граничног
слоја
одвајање
у
граничном
слоју
?
важеће
решење
целокупног
струјног
проблема
почетак
25
механике
флуида
која
је
најдубље
продрла
у
реалне
струјне
проблеме
.
Иницирање
стварне
слике
струјања
и
могућност
комбиновања
и
описивања
нових
реалних
и
имагинарних
струјања
захтева
изучавање
потенцијалног
струјања
-
најзначајнијег
представника
кретања
идеалног
флуида
. (
Само
понекад
потенцијално
струјање
примењује
се
на
слабо
компресибилне
флуиде
.)
Гранични
слој
се
изучава
нешто
касније
јер
се
његов
практичан
значај
манифестује
кроз
отпоре
тела
која
се
крећу
кроз
флуид
-
најважнију
дисциплину
примењене
механике
флуида
у
техници
.
Функција
потенцијала
брзине
φφφφ
и
струјна
функција
ψ
ψ
ψ
ψ
Полазна
претпоставка
потенцијалног
струјања
да
невискозан
флуид
условљава
непостојање
вртлога
може
да
се
докаже
математичким
путем
или
објасни
на
следећи
начин
.
Да
би
се
у
невтрложном
струјању
појавили
вртлози
потребан
је
момент
,
међутим
,
у
одсуству
тангенцијалних
сила
(
невискозан
флуид
)
једине
силе
које
постоје
су
нормалне
(
притисак
који
делује
кроз
геометријско
тежиште
флуидног
делића
)
и
тежина
(
која
делује
кроз
центар
месе
флуидног
делића
).
Ниједна
од
ових
сила
у
флуиду
једнолике
густине
не
може
да
изазове
момент
који
ће
обртати
флуидни
делић
.
У
одсуству
вискозности
вртложење
не
може
да
се
појави
ни
уколико
кроз
миран
флуид
пролази
чврсто
тело
,
тако
да
важи
већ
раније
исказан
постулат
да
уколико
је
почетно
струјање
без
вртлога
оно
ће
и
остати
безвртложно
.
На
овом
месту
треба
поновити
да
стационарно
и
невртложно
струјање
Ојлерову
једначину
редукује
на
основни
облик
Бернулијеве
једначине
која
важи
за
цело
струјно
поље
.
2
d
.
2
v
p
gz
const
ρ
+
+
=
∫
Услов
невртложности
[
]
,
rot
0
v
v
∇
=
=
захтева
v
φ
= ∇

27
Слика
3.
Израчунавање
запреминског
протока
између
две
струјнице
За
одређену
струјницу
ψ
=const.
важи
d
d
d
0
x
y
x
y
ψ
ψ
ψ
∂
∂
=
+
=
∂
∂
тј
.
.
d
d
y
const
x
v
y
x
x
v
y
ψ
ψ
ψ
=
∂
∂
=
=
∂
∂
представља
нагиб
струјнице
.
Линије
φ
=const.
које
се
називају
еквипотенцијалним
линијама
нормалне
су
на
струјнице
пошто
се
из
d
d
d
0
x
y
x
y
φ
φ
φ
∂
∂
=
+
=
∂
∂
и
.
d
d
x
const
y
v
y
x
x
v
y
φ
φ
φ
=
∂
∂
=
= −
∂
∂
види
да
је
њихов
нагиб
реципрочно
негативан
према
струјницама
.
Ова
законитост
губи
се
само
у
тачкама
сингуларитета
(
тачке
у
којима
су
v
x
=0
и
v
y
=0).
На
слици
4
приказане
су
струјнице
и
еквипотенцијалне
линије
за
спољашње
и
унутрашње
струјање
.
За
невртложно
струјање
важи
28
2
0
y
x
z
v
v
x
y
x
x
y
y
ψ
ψ
ω
ψ
∂
∂
∂
∂
∂
∂
=
−
=
−
−
−
= −∆ =
∂
∂
∂
∂
∂
∂
односно
струјна
функција
мора
да
задовољи
Лапласову
једначину
као
и
функција
φ
2
2
0
x
y
ψ
ψ
∂
∂
+
=
∂
∂
Функције
φ
и
ψ
различите
су
,
иако
задовољавају
исту
једначину
због
различитих
граничних
услова
.
Слика
4.
Еквипотенцијалне
линије
и
струјнице
за
пример
спољашњег
и
унутрашњег
струјања
Струјна
функција
се
подудара
са
обликом
чврсте
површине
.
Тамо
је
ψ
=const.
Често
је
корисно
неку
интересантну
струјницу
или
више
њих
прогласити
чврстом
површином
.
Тада
остатак
струјница
представља
потенцијално
поље
које
се
анализира
.
Уобичајена
су
три
начина
за
решавање
потенцијалних
струјних
проблема
.
Сва
три
су
математичка
,
без
потребе
за
експерименталним
информацијама
.
То
су
:
o
метод
сингуларитета
o
метод
комплексних
променљивих
,
o
директни
апроксимативни
методи
.
Метод
сингуларитета
формира
потенцијално
струјање
комбинацијом
једноликог
струјања
и
изолованих
сингуларитета
.
Сингуларитети
(
тачке
у
којима
једначине
струјања
престају
да
важе
)
су
:
извор
-
понор
,
вртлог
и
двопол
.
Једнолико
струјање
није
сингуларитет
али
може
да
се

30
( )
(
)
f z
f x
iy
i
ξ
η
=
+
= +
Најважнији
тип
комплексних
функција
су
аналитичке
функције
које
имају
једнозначан
извод
у
свакој
тачки
.
За
аналитичке
функције
важе
следеће
две
теореме
:
o
Да
би
аналитичка
функција
( )
(
)
f z
f x
iy
i
ξ
η
=
+
= +
била
аналитичка
,
неопходно
је
да
x
y
ξ
η
∂
∂
=
∂
∂
и
y
x
ξ
η
∂
∂
= −
∂
∂
.
Ово
су
Коши
-
Риманове
једначине
.
o
Обе
компоненте
аналитичких
функција
задовољавају
Лапласову
једначину
,
тако
да
је
2
2
2
2
0
x
y
ξ
ξ
∂
∂
+
=
∂
∂
и
2
2
2
2
0
x
y
η
η
∂
∂
+
=
∂
∂
Пошто
потенцијал
брзине
и
струјна
функција
задовољавају
Коши
-
Риманове
једначине
и
Лапласову
једначину
;
x
y
v
v
x
y
y
x
φ
ψ
φ
ψ
∂
∂
∂
∂
=
=
=
= −
∂
∂
∂
∂
2
2
2
2
2
2
2
2
0;
0
x
y
x
y
φ
φ
ψ
ψ
∂
∂
∂
∂
+
=
+
=
∂
∂
∂
∂
оне
комбинацијом
дају
комплексну
функцију
(
комплексни
потенцијал
)
( )
( , )
( , )
W z
x y
i
x y
φ
ψ
=
+
Одавде
следи
закључак
да
се
реалан
и
имагинаран
део
једне
произвољне
аналитичке
функције
W
може
увек
сматрати
потенцијалом
брзине
φ
и
струјном
функцијом
ψ
раванског
потенцијалног
струјања
.
Из
комплексног
потенцијала
непосредно
се
налази
брзина
струјања
не
одређујући
претходно
ни
потенцијал
φ
ни
струјну
функцију
ψ
.
d
d
d
d
d
d
d
d
d
d
1
1
(1
)
(1
)
1
1
x
y
y
x
x
y
x
y
y
i
i
y
i
W
i
x
y
x
y
x
x
z
x
i y
iy
iy
v
v y
iv
iv y
v
iy
iv
iy
v
iv
v
iy
iy
φ
φ
ψ
ψ
φ
ψ
φ
ψ
∂
∂
∂
∂
′
′
+
+
+
+
+
∂
∂
∂
∂
=
=
=
=
′
′
+
+
+
′
′
′
′
+
−
+
+
−
+
=
=
=
−
=
′
′
+
+
која
се
назива
комплексном
брзином
.
Апсолутне
вредности
комплексне
брзине
31
2
2
d
d
x
y
W
v
v
v
v
z
=
=
+
=
представља
интензитет
вектора
брзине
v
.
За
позитиван
угао
θ
између
вектора
брзине
v
и
x
-
осе
cos ,
sin
x
y
v
v
v
v
θ
θ
=
=
следи
d
(cos
sin )
d
i
W
v
v
i
ve
z
θ
θ
θ
−
= =
−
=
.
Зауставне
тачке
у
струјању
налазе
се
на
местима
где
је
d
0
d
W
z
=
.
Комплексни
потенцијали
основних
струјања
Комплексни
потенцијал
униформне
струје
под
углом
α
у
односу
на
x
-
осу
(
слика
5
)
0
0
0
0
cos
sin
(
cos
sin )
W
i
v x
v y
i v y
v x
φ ψ
α
α
α
α
= +
=
+
+
−
тј
.
(
)
0
0
0
(
)(cos
sin )
i
i
W
v x
iy
i
v ze
v re
α
θ α
α
α
−
−
=
+
−
=
=
,
0
0
cos(
);
sin(
)
v r
v r
φ
θ α
ψ
θ α
=
−
=
−
.
Слика
5.
Униформна
струја
Комплексни
потенцијал
осамљеног
извора
(
понора
) (
слика
6)

33
Слика
7.
Струјање
у
пољу
вртлога
( )
ln
(
ln )
2
2
,
ln
2
2
2
i
W z
z
i
r
r
v
r
θ
π
π
φ
θ
ψ
π
π
π
Γ
Γ
=
=
− +
Γ
Γ
= −
=
Γ
=
За
вртложење
у
правцу
казаљке
на
сату
(
негативно
)
циркулација
је
мања
од
нуле
,
Γ
<0.
Када
произвољна
затворена
крива
не
обухвата
координатни
почетак
(
сингуларитет
),
циркулација
је
нула
Γ
=0,
а
када
га
обухвата
Γ
=
const.
,
јер
је
потенцијал
вишезначна
функција
d
,
0,
d
2
2
r
W
i
v
v
v
z
z
r
θ
π
π
Γ
Γ
=
=
=
= −
Комплексни
потенцијал
двопола
Варирањем
раздаљине
извора
и
понора
у
јадноликој
струји
добијају
се
екстремна
струјања
.
Уколико
се
извор
и
понор
удаљавају
у
бесконачност
добија
се
струјање
преко
танке
равне
плоче
,
што
одговара
једноликом
струјању
.
Уколико
се
извор
и
понор
споје
,
добија
се
34
струјање
око
двопола
које
одговара
опструјавању
малог
кружног
цилиндра
.
Опструјавање
кружног
цилиндра
коначне
димензије
је
врло
важно
,
а
добија
се
повећањем
издашности
извора
и
понора
уз
задржавање
а
→
0,
односно
одржавањем
константног
производа
ε
а
док
а
→
0.
Овако
се
добија
нов
тип
сингуларитета
који
се
назива
двопол
(
слика
8
).
Слика
8.
Струјање
у
пољу
двопола
Извор
и
понор
истих
издашности
±
ε
налазе
се
на
реалној
оси
на
растојању
±
а
од
координатног
почетка
.
Комплексни
потенцијал
збирног
струјања
износи
ln(
)
ln(
)
ln
2
2
2
z
a
W
z
a
z
a
z
a
ε
ε
ε
π
π
π
+
=
+
−
−
=
−
.
Ако
се
уведе
момент
двопола
0
lim 2
.
a
M
a
const
ε
→
=
=
и
потражи
гранична
вредност
када
а
→
0,
добија
се
0
lim
ln
4
2
a
M
z
a
M
W
a
z
a
z
π
π
→
+
=
=
−
комплексни
потенцијал
двопола
смештеног
у
координатном
почетку
.
Из
комплексног
потенцијала
следе
2
2
2
2
cos ,
sin
2
2
2
2
M
x
M
M
y
M
x
y
r
x
y
r
φ
θ
ψ
θ
π
π
π
π
=
=
= −
= −
+
+
Ако
је
двопол
оријентисан
под
неким
углом
α
према
x
-
оси
,
а
налази
се
у
тачки
z
1
,
онда
се
W
,
φ
и
ψ
дати
са

36
Слика
10.
Струјање
у
пољу
двопола
Комплексни
потенцијал
облика
степене
функције
Аналитичка
функција
облика
(cos
sin
)
n
n
W
az
ar
n
i
n
θ
θ
=
=
+
са
cos
n
ar
n
φ
θ
=
и
sin
n
ar
n
ψ
θ
=
,
представља
струјање
унутар
и
око
две
чврсте
равни
које
међусобно
заклапају
различите
углове
.
Нулте
струјнице
(
чврсте
равни
)
добијају
се
из
услова
ψ
=0,
тј
. sin
0
n
θ =
односно
,
0,1,...
n
k
k
θ
π
=
=
,
па
је
θ
=0
и
θ
=
π
/n
.
На
слици
11
приказани
су
карактеристични
случајеви
струјања
.
Комплексна
брзина
је
:
1
1
(
1)
d
d
n
n
i n
W
v
naz
nar
e
z
θ
−
−
−
=
=
=
а
интензитет
брзине
1
n
v
v
na r
−
=
=
,
тако
да
у
зауставној
тачки
(
v
=0
или
r
=0)
постоје
следећи
случајеви
:
за
n
> 1,
v
= 0
за
n
= 1,
v
=
a
за
n
< 1,
v
=
∞
.
37
Слика
11.
Изглед
струјница
степене
функције
за
различите
вредности
n
Одавде
следи
закључак
да
у
тачкама
дисконтинуитета
чврсте
површине
(
различите
тангенцијале
),
брзина
идеалног
флуида
једнака
је
нули
када
је
угао
са
флуидне
стране
мањи
од
π
,
а
брзина
је
бесконачна
када
је
угао
већи
од
π
.
Код
реалног
флуида
за
θ
<
π
наступиће
одвајање
флуида
од
чврсте
површине
због
успоравања
струјања
уз
појаву
струјећег
вртлога
у
углу
,
а
за
θ
>
π
ако
је
промена
тангенцијале
нагла
на
месту
дисконтинуитета
долази
до
одвајања
струје
.
Комплексни
потенцијали
сложених
струјања
Функције
комплексних
потенцијала
могу
да
се
сабирају
.
На
тај
начин
долази
се
до
великог
броја
разноликох
струјања
која
се
лакше
проучавају
.
Највећи
значај
имају
онда
када
нулта
струјница
образује
затворену
криву
јер
тада
може
да
се
замисли
да
она
представља
чврсту
површину
и
да
флуид
струји
око
ње
.
Струјање
око
полутела
Струјање
око
полутела
представља
збир
једнолике
струје
и
струјања
у

39
2
2
2
0
0
0
1
1
4
p
v
c
x
v
v
v
x
y
ε
ε
π
π
= −
= −
+
+
.
Слика
13.
Промена
притиска
и
брзине
при
струјању
у
пољу
једнолике
струје
и
понора
Струјање
око
кружног
цилиндра
Комбинација
једноликог
струјања
са
двополом
представља
струјање
идеалног
флуида
око
кружног
цилиндра
.
Опструјавање
кружног
цилиндра
је
значајно
за
анализу
отпора
тела
потопљених
у
реалном
флуиду
.
Струјање
без
циркулације
Комплексни
потенцијал
,
струјна
функција
и
функција
потенцијала
брзине
,
за
усвојену
брзина
долазеће
струје
v
0
и
момент
двопола
М
дати
су
следећим
изразима
40
0
2
M
W
v z
z
π
=
+
;
0
0
2
2
1
cos
2
2
M
M
x v
v r
x
y
r
φ
θ
π
π
=
+
=
+
+
;
0
0
2
2
1
sin
2
2
M
M
y v
v r
x
y
r
ψ
θ
π
π
=
−
=
−
+
.
За
нулту
струјницу
(
ψ
=0)
добија
се
2
2
2
0
0;
;
2
M
y
x
y
R
R
v
π
=
+
=
=
.
Струјнице
су
симетричне
у
односу
на
x
и
y
-
осу
(
видети
слику
14
).
Могу
се
конструисати
графички
,
слагањем
два
компонентна
струјања
(
Ренкинова
метода
).
Слика
14.
Ациклично
струјање
око
кружног
цилиндра
Сменом
0
2
M
v R
π
=
струјна
функција
може
се
изразити
зависно
од
полупречника
кружног
цилиндра
2
0
sin
R
v
r
r
ψ
θ
=
−
и
2
0
R
W
v
z
z
=
−
одакле
је
2
2
0
0
2
2
1
1
cos ;
1
sin
r
R
R
v
v
v
v
r
r
r
r
θ
ψ
ψ
θ
θ
θ
∂
∂
=
=
−
= −
= −
+
∂
∂
.
За
контуру
цилиндра
је
0
r
v
=
и
0
2
sin
v
v
v
θ
θ
= = −
;
а
зауставне
тачке
су

42
Слика
15.
Распоред
притиска
при
опструјавању
кружног
цилиндра
без
циркулације
идеалним
флуидом
Струјање
са
циркулацијом
Цилиндар
који
се
обрће
и
истовремено
налази
у
једноликој
струји
генерише
узгонску
силу
,
без
обзира
да
ли
је
у
питању
идеалан
или
реалан
(
вискозан
)
флуид
.
Симулација
овакве
врсте
струјања
постиже
се
постављањем
усамљеног
линијског
вртлога
у
осу
цилиндра
.
Уколико
је
једнолико
струјање
усмерено
у
правцу
позитивне
x
-
осе
,
за
добијање
позитивне
узгонске
силе
(
позитивна
y
-
оса
)
мора
се
увести
негативна
циркулација
(
−Γ
).
Комплексни
потенцијал
,
струјна
функција
и
потенцијал
брзине
дати
су
са
:
2
0
ln
2
R
i
W
v
z
z
z
π
Γ
=
+
+
2
0
2
1
cos
2
R
v r
r
φ
θ
θ
π
Γ
=
+
−
2
0
2
1
sin
ln
2
R
v r
r
r
ψ
θ
π
Γ
=
−
+
.
Из
комплексне
брзине
2
0
2
d
1
d
2
W
R
i
v
v
z
z
z
π
Γ
=
=
−
+
43
координате
зауставних
тачака
одређују
се
решавањем
једначине
2
2
0
0
2
2
0
v z
i z
v R
π
π
+ Γ −
=
тј
.
дискриминанте
2
2
2
2
0
16
v R
π
∆ = −Γ +
за
коју
важе
три
случаја
(
видети
слику
14
):
а
)
0
∆ >
;
0
4
v R
π
Γ <
-
брзина
је
једнака
нули
у
тачкама
:
2
2
1
2
2
0
0
4
16
i
z
R
v
v
π
π
Γ
Γ
= −
−
−
и
2
2
2
2
0
0
4
16
i
z
R
v
v
π
π
Γ
Γ
= −
+
−
Ове
тачке
су
симетричне
према
y
-
оси
,
а
негативан
знак
ординате
показује
да
леже
на
кружници
испод
x
-
осе
,
како
је
показано
на
слици
16
.
б
)
0
∆ =
;
0
4
v R
π
Γ =
-
зауставне
тачке
се
поклапају
1
2
0
4
i
z
z
v
π
Γ
=
= −
в
)
0
∆ <
;
0
4
v R
π
Γ >
-
зауставне
тачке
леже
на
ординатној
оси
,
али
су
различито
удаљене
од
координатног
почетка
.
Тачка
која
пада
унутар
кружнице
без
значаја
је
за
опструјавање
.

45
2
2
0
0
0
1
1
2
sin
2
2
2
p
p
v
v
R
ρ
ρ
θ
π
Γ
=
+
−
−
−
.
Елементарна
сила
притиска
на
цилиндар
(
слика
17
),
где
је
b
дужина
цилиндра
је
:
d
d
F
pbR
θ
= −
а
хоризонтална
и
вертикална
компонента
d
cos d
x
F
pbR
θ θ
−
=
; d
sin d
y
F
pbR
θ θ
−
=
односно
:
2
0
cos d
0
x
F
pbR
π
θ θ
= −
=
∫
,
2
0
0
sin d
y
F
pbR
v
b
π
θ θ
ρ
= −
= −
Γ
∫
Слика
17.
Елементарна
сила
притиска
на
цилиндар
и
њене
компоненте
Сила
0
x
F
=
јер
је
распоред
притиска
симетричан
у
односу
на
y
-
осу
.
Из
0
y
F
v
ρ
= −
Γ
за
негативну
циркулацију
и
једнолико
струјање
смера
позитивне
x
-
осе
добија
се
узгонска
сила
(
вертикално
навише
усмерена
сила
нормална
на
правац
брзине
)
којом
флуид
,
чак
и
идеалан
,
дејствује
на
кружни
цилиндар
.
Постојање
ове
силе
указује
да
је
притисак
мањи
на
горњој
половини
кружнице
(
цилиндра
јединичне
висине
)
него
на
доњој
.
Код
реалног
флуида
циркулација
се
изазива
обртањем
цилиндра
.
Обимна
брзина
цилиндра
U
R
ω
=
у
вези
је
са
циркулацијом
2
U
R
π
Γ
=
тако
да
је
46
0
2
y
L
F
F
v U Rb
π
ρ
=
=
(
уобичајена
ознака
за
силу
узгона
је
F
L
)
а
коефицијент
узгона
је
2
0
0
2
1
(2
)
2
L
L
F
U
c
v
v
Rb
π
ρ
=
=
Овако
добијена
узгонска
сила
са
кружним
цилиндром
који
се
обрће
има
ограничену
примену
(
нпр
.
коришћењем
Флетнерових
ротора
на
палуби
брода
могућа
је
уштеда
горива
до
30 % (
слика
18
)).
Уобичајено
је
да
се
узгонска
сила
реализује
аеропрофилима
.
Слика
18.
Брод
Брукау
са
Флетнеровим
роторима
(1924.)

48
Општи
облик
једначине
о
промени
количине
кретања
Пошто
једначина
o
промени
количине
кретања
важи
само
за
инерцијски
координатни
систем
(
непокретан
или
се
креће
једноликом
брзином
),
промена
количине
кретања
изражаваће
се
у
односу
на
апсолутан
-
непокретан
координатни
систем
.
Количина
кретања
(
импулс
)
елементарне
масе
је
:
d
d
d
K
v m
v V
ρ
=
=
За
укупну
масу
флуида
,
запремине
V
(
t
),
важи
( )
d .
V t
K
v V
ρ
=
∫
Према
другом
Њутновом
закону
сума
спољашњих
сила
у
једном
правцу
једнака
је
промени
количине
кретања
у
истом
правцу
,
( )
d
d
d
d
d
V t
K
v V
F
t
t
ρ
=
=
∑
∫
Применом
трансформације
временског
извода
на
једначину
о
промени
количине
кретања
,
добија
се
:
( )
(
)
( )
d
d
d
d
,
d
d
d
V t
V
A
v
K
v V
V
v v n
A
F
t
t
t
ρ
ρ
ρ
∂
=
=
+
=
∂
∑
∫
∫
∫
Први
члан
у
трећем
делу
једнакости
описује
локалну
промену
количине
кретања
у
запремини
V
,
за
чега
је
потребно
познавање
струјних
величина
у
унутрашњости
запремине
.
Други
члан
даје
резултујуће
струјање
кроз
граничне
површине
,
за
чега
је
потребно
познавање
свих
променљивих
само
на
граничним
површинама
запремине
V
.
За
стационарна
струјања
запремински
интеграл
једнак
је
нули
,
тако
да
су
потребни
само
струјни
подаци
на
границама
контролне
запремине
,
а
израз
постаје
:
(
)
,
d
A
v v n
A
F
ρ
=
∑
∫
Ако
се
сила
количине
кретања
дефинише
као
:
(
)
,
d
K
A
F
v v n
A
ρ
= −
∫
49
онда
(
)
,
d
A
v v n
A
F
ρ
=
∑
∫
може
да
се
напише
као
:
0
K
F
F
+
=
∑
.
Ово
је
основна
једначина
за
решавање
проблема
.
Слика
1.
Дефинисање
струјних
величина
на
границама
контролне
запремине
За
силу
количине
кретања
важи
да
је
локално
паралелна
са
вектором
брзине
и
увек
је
усмерена
ка
унутрашњости
контролне
запремине
(
слика
1
).
(
)
d
,
d
K
F
v v n
A
ρ
= −
У
једначинама
сума
сила
F
∑
представља
збир
свих
сила
које
су
узрок
промени
количине
кретања
.
Ове
силе
су
:
o
активне
-
запреминске
силе
(
сила
Земљине
теже
,
Њутнова
привлачна
сила
);
o
силе
притиска
-
деле
се
на
спољашње
и
унутрашње
.
Спољашње
силе
потичу
од
маса
које
су
ван
контролне
запремине
.
Унутрашње
силе
налазе
се
у
маси
посматране
запремине
;
њихово
узајамно
дејство
се
поништава
,
па
је
њихов
збир
једнак
нули
за
целу
запремину
;
o
силе
реакције
које
се
такође
деле
на
спољашње
и
унутрашње
,
али
и
овде
остају
реакције
од
спољашњих
маса
,
и
реакције
граничних
површина
-
силе
којима
се
граничне
површине
супротстављају
кретању
флуида
.
Једначина
о
промени
количине
кретања
користи
се
при
стационарном
струјању
кроз
(
преко
)
непокретних
елемената
или
елемената
који
се
крећу
константном
брзином
,
а
при
томе
се
маса
и
брзина
флуида
у

51
Примери
примене
једначине
о
промени
количине
кретања
на
непокретне
елементе
За
примену
једначине
неопходно
је
уочити
део
флуидног
простора
-
контролну
запремину
у
којој
долази
до
промене
неких
од
наведених
величина
;
па
је
затим
ограничити
од
остале
средине
граничним
или
замишљеним
површинама
–
пресецима
.
Уз
претпоставку
да
је
ток
струјница
познат
,
пресеци
се
постављају
нормално
на
њих
,
а
утицаји
средине
испред
и
иза
замењују
се
силама
притиска
p
i
A
i
које
су
уперене
ка
контролној
запремини
.
Непокретни
елементи
са
једним
улазом
и
излазом
флуида
Поступак
одређивања
силе
којом
флуид
дејствује
на
зидове
цеви
биће
представљен
на
примеру
обичног
цевног
колена
,
од
90
°
константне
површине
пресека
А
кроз
које
протиче
количина
флуида
Q
густине
ρ
(
слика
2
):
Слика
2.
Примена
једначине
о
промени
количине
кретања
на
колено
1.
Издвоји
се
контролна
запремина
(
пресеци
1-1, 2-2
и
зидови
колена
;
индекс
1-1
везује
се
за
пресек
кроз
који
улази
флуид
,
а
2-2
за
пресек
кроз
који
флуид
напушта
контролну
запремину
).
52
2.
Означе
се
смерови
брзина
у
улазном
и
излазном
пресеку
,
уцртају
се
вектори
сила
количине
кретања
(
F
k
1
и
F
k
2
)
и
силе
притисака
које
замењују
утицај
флуидне
струје
испред
и
иза
колена
(
F
p
1
и
F
p
2
).
3.
Произвољно
се
претпоставе
смерови
координатних
оса
x
и
y
.
4.
Уноси
се
сила
реакција
веза
F
.
То
је
сила
којом
веза
(
лук
,
млазник
,
плоча
и
др
.)
дејствује
на
флуид
.
Она
се
најчешће
разлаже
на
компоненте
у
правцу
координатних
оса
(
F
x
и
F
y
).
Општи
облик
једначине
о
промени
количине
кретања
у
векторском
облику
гласи
1
2
1
2
0
k
k
p
p
F
F
F
F
F
+
+
+
+
=
.
Пројектована
на
x
-
осу
у
случају
непокретне
контролне
запремине
испуњене
стационарном
струјом
нестишљивог
флуида
уз
занемаривање
локалног
губитка
(
p
2
=
p
1
=
p
),
своди
се
на
2
0
x
p A
F
Qv
ρ
−
+
−
=
где
је
производ
ρ
Qv
интензитет
силе
количине
кретања
у
случају
стационарне
униформне
струје
,
тј
.
за
p
2
=
p
Qv
pA
F
x
ρ
+
=
.
Слично
је
за
y
-
осу
1
0
y
p A
F
Qv
ρ
−
+
−
=
тј
.
y
F
pA
Qv
ρ
=
+
.
Пошто
је
у
резултату
за
F
x
и
F
y
добијена
позитивна
вредност
,
реакције
везе
су
добро
претпостављеног
смера
.
Силе
којима
флуид
дејствује
на
зидове
колена
(
x
F
′
и
y
F
′
)
истог
су
интензитета
као
и
одређене
силе
реакције
колена
,
само
су
супротног
смера
.
Дакле
,
Qv
pA
F
x
ρ
+
=
′
←
и
Qv
pA
F
y
ρ
+
=
′
↓
а
резултујућа
сила
(
)
2
Qv
pA
F
ρ
+
=
′
.

54
Слика
4.
Дефинисање
контролне
запремине
при
струјању
кроз
нагло
проширење
У
пресеку
1’-1’ (
слика
4
)
влада
исти
притисак
као
и
у
пресеку
1-1
јер
су
то
два
пресека
који
су
бесконачно
близу
један
другом
.
Једначину
о
промени
количине
кретања
треба
поставити
за
пресеке
1’-1’
и
2-2,
јер
се
на
тај
начин
елиминише
величина
А
1
,
а
у
физичкој
интерпретацији
проблема
не
прави
се
никаква
апроксимација
због
тога
што
је
носилац
притиска
струја
флуида
,
а
не
геометрија
простора
(
пресек
1’-1’)
кога
струјање
делимични
испуњава
.
1
2
2
2
1
2
0
p A
p A
Qv
Qv
ρ
ρ
−
+
−
=
.
Једначина
о
промени
количине
кретања
допуњава
се
Бернулијевом
једначином
за
пресеке
1-1
и
2-2. (
При
протицању
реалног
флуида
кроз
контролну
запремину
,
губици
енергије
не
улазе
експлицитно
у
једначину
о
промени
количине
кретања
,
јер
она
описује
промене
на
границама
контролне
запремине
.
Експлицитно
ови
губици
се
јављају
у
Б
.
ј
.,
јер
она
важи
за
струјнице
унутар
контролне
запремине
.
У
једначини
о
промени
количине
кретања
губици
енергије
су
имплицитно
представљени
вредностима
притисака
и
брзина
из
Б
.
ј
.).
2
2
1
1
2
2
2
2
i
p
v
p
v
gh
ρ
ρ
+
=
+
+
Из
једначине
о
промени
количине
кретања
замењује
се
израз
1
2
(
) /
p
p
ρ
−
у
Б
.
ј
.
2
2
2
1
2
1
2
(
)
2
i
Q v
v
v
v
gh
A
−
−
=
+
2
2
2
1
2
2
1
(
)
2
2
i
v
v
v v
v
gh
−
−
+
=
55
2
2
2
1
1 2
2
2
2
2
i
v
v
v v
gh
+
−
=
2
1
2
(
)
2
i
v
v
gh
−
=
.
Млазник
и
дифузор
са
слободним
истицањем
У
овом
примеру
потребно
је
одредити
силу
којом
флуид
делује
на
везу
млазника
са
цевоводом
при
слободном
истицању
флуида
из
њега
.
Брзина
и
притисак
су
познати
у
пресецима
А
1
и
А
2
,
као
и
ρ
=const.
и
p
2
=
p
a
.
Погодан
избор
контролне
запремине
дат
је
на
слици
5
,
док
су
на
слици
6
приказане
силе
које
делују
на
њу
.
Слика
5.
Контролна
запремина
при
струјању
кроз
млазник
Слика
6.
Силе
на
границама
контролне
запремине
2
2
1
1
1
2
2
2
1
1
1
2
1
,
,
,
k
k
p
p
a
B
x
F
v A F
v A
F
p A F
p A
F
F
ρ
ρ
=
=
=
=
=
Смер
реакције
ослонца
произвољно
је
усвојен
,
а
одређује
се
решавањем
једначине
о
промени
количине
кретања
у
правцу
х
осе
.
2
2
1
1
1
1
2
2
1
0
a
x
v A
p A
v A
p A
F
ρ
ρ
+
−
−
+
=
.
Из
једначине
континуитета
1
1
2
2
v A
v A
=
следи
(
)
2
2
2
2
2
1
1
1
x
a
A
F
v
A
p
p A
A
ρ
=
−
+
−
Ако
се
додатно
претпостави
да
је
струјање
идеално
,
из
Бернулијеве
једначине
следи
:

57
2
min
A
gh
A
v
=
.
Ако
се
искористи
Торичелијева
формула
за
брзину
истицања
2
v
gh
=
,
следи
:
0,5
min
A
A
=
.
Слика
8.
Струјање
кроз
Бордин
наглавак
Слика
9.
Контролна
запремина
при
истицању
из
Бординог
наглавка
Експериментима
су
добијене
веће
вредности
: 0,61
до
0,64.
Ово
потиче
отуда
,
што
је
,
због
губитака
,
брзина
истицања
мања
од
Торичелијеве
идеалне
брзине
.
Отпор
полутела
у
каналу
Разматра
се
нестишљиво
,
невискозно
опструјавање
непокретног
полутела
у
каналу
(
слика
10
).
Сила
отпора
којом
флуид
делује
на
тело
може
да
се
одреди
из
једначине
о
промени
количине
кретања
.
За
одређивање
ове
силе
потребно
је
познавање
услова
на
задњој
страни
тела
.
Код
полутела
није
могуће
да
се
дефинишу
ти
услови
,
зато
се
претпоставља
да
у
попречном
пресеку
(2),
довољно
далеком
од
врха
тела
,
влада
притисак
p
=
p
2
.
58
x
→
Слика
10.
Опструјавање
полутела
у
каналу
Једначина
о
промени
количине
кретања
гласи
:
2
2
1
1
2
1
1
1
2
2
0
x
p A
p A
F
v A
v A
ρ
ρ
−
−
+
−
=
где
је
2
1
A
A
A
=
−
,
а
једначина
континуитета
је
1
1
2
2
v A
v A
=
.
Сила
x
F
је
реакција
ослонца
,
а
сила
отпора
којом
вода
делује
на
тело
је
x
x
F
F
′ =
:
2
2
1
1
2
2
1
2
1
(
)
x
F
v A
v A
p
p A
ρ
ρ
=
−
+
−
.
Бернулијева
једначина
и
једначина
континуитета
за
разматрано
невискозно
струјање
дају
:
(
)
2
2
2
2
1
1
2
2
1
1
2
2
1 ,
2
2
A
p
p
v
v
v
A
ρ
ρ
−
=
−
=
−
односно
2
2
2
1
1
1
1
1
2
2
1
1
2
2
1
x
A
A
A
F
v A
v A
A
A
A
ρ
ρ
=
−
=
−
.
Сила
којом
вода
делује
на
полутело
(
x
F
′
)
представља
се
познатим
изразом
2
1
2
x
D
F
v c A
ρ
′ =
,
где
је
брзина
v
=
v
1
.
Бездимензијски
коефицијент
c
D
представља
коефицијент
отпора
тела
.
Изједначавањем
последње
две
једначине
следи
да
је

60
а
из
Бернулијевих
једначина
(
водећи
рачуна
да
је
p
1
=
p
4
=
p
0
притисак
околне
средине
)
за
пресеке
1-2
и
3-4,
добија
се
редом
2
2
2
2
1
1
2
3
4
4
1
1
1
1
,
2
2
2
2
p
v
p
v
p
v
p
v
ρ
ρ
ρ
ρ
+
=
+
+
=
+
.
Слика
11.
Примена
једначине
о
промени
количине
кретања
на
пропелер
Сабирањем
претходне
две
једначине
,
водећи
рачуна
да
су
p
1
=
p
4
=
p
0
,
следи
2
2
3
2
4
1
1
(
)
2
p
p
v
v
ρ
−
=
−
што
у
комбинацији
са
једначином
4
1
3
2
(
)
v v
v
p
p
ρ
−
=
−
даје
1
4
2
v
v
v
+
=
Брзина
кроз
роторски
диск
,
значи
,
има
средњу
вредност
брзине
узводно
и
низводно
од
пропелера
.
61
Коефицијент
корисног
дејства
пропелера
добија
се
из
односа
корисне
снаге
(
добијене
из
једначине
о
промени
количине
кретања
)
и
утрошене
снаге
.
Корисна
снага
је
(
)
1
4
1
1
F v
Q v
v v
ρ
⋅ =
−
као
да
се
пропелер
креће
брзином
v
1
кроз
миран
флуид
.
Уложена
снага
је
(
)
2
2
4
1
2
Q
v
v
ρ
−
.
Степен
корисног
дејства
пропелера
је
однос
корисне
и
уложене
енергије
(
)
(
)
4
1
1
1
1
2
2
4
1
4
1
2
2
Q v
v v
v
v
Q
v
v
v
v
v
ρ
η
ρ
−
=
=
=
+
−
.
Ако
се
са
∆
v
означи
4
1
v
v
v
∆ =
−
,
тј
.
4
1
v
v
v
= ∆ +
,
степен
корисности
може
да
се
напише
у
другом
облику
1
1
2
v
v
v
η =
∆
+
одакле
се
види
да
се
највећи
степен
корисности
постиже
за
најмање
повећање
брзине
флуидне
струје
.
Непокретни
елементи
са
различитим
бројем
улаза
и
излаза
флуида
Једначина
о
промени
количине
кретања
,
без
икаквих
ограничења
,
може
да
се
примени
за
различит
број
улаза
и
излаза
флуида
.
Под
претпоставком
да
су
притисци
у
свим
означеним
улазним
и
излазним
пресецима
представљене
контролне
запремине
исти
,
ЈПКК
за
претпостављене
силе
реакције
гласи
:
2
1
1
1
n
m
i
xi
j
xj
i
j
X
Q v
Q v
ρ
=
=
=
−
∑
∑
∑
2
1
1
1
n
m
i
yi
j
yj
i
j
Y
Q v
Q v
ρ
=
=
=
−
∑
∑
∑
тј
.

63
Слика
13.
Примена
једначине
о
промени
количине
кретања
млазни
мотор
Због
тога
настаје
сила
количине
кретања
,
чија
компонента
у
х
правцу
износи
:
(
)
(
)
2
,
M
,
d
K x
x
m
F
v v n
A
v m
v
A
A
ρ
ρ
∞
∞ ∞
∞
= −
=
=
−
∫
ЈПКК
гласи
:
(
)
(
)
2
2
2
2
0
m
m m
m
m
R
v A
v
A
A
v A
v
A
A
F
ρ
ρ
ρ
ρ
∞ ∞
∞ ∞
∞
∞ ∞
+
−
−
−
−
+
=
где
је
F
R
сила
реакције
коју
трпи
мотор
.
Она
је
једнака
(
)
2
2
R
m m
m
m
m
F
v A
v A
m
v
v
ρ
ρ
∞ ∞
∞
∞
=
−
=
−
.
Погонска
сила
R
F
F
= −
.
Директно
је
пропорционална
масеном
протоку
и
порасту
брзине
у
млазу
у
односу
на
околину
.
Елементи
који
се
крећу
константном
брзином
Сила
којом
флуид
дејствује
на
покретан
елемент
,
и
под
чијим
се
дејством
тај
елемент
креће
,
може
да
буде
:
o
активна
сила
флуидног
млаза
(
нпр
.
турбинске
лопатице
)
o
реактивна
сила
флуидног
млаза
(
нпр
.
ракетни
мотор
).
Ова
покретачка
сила
иста
је
непокретном
(
апсолутном
)
и
покретном
(
релативном
)
систему
ако
се
овај
креће
константном
брзином
,
под
условом
да
се
изражава
преко
одговарајућих
апсолутних
и
релативних
величина
на
улази
и
излазу
тог
елемента
.
Када
је
покретачка
сила
активна
сила
млаза
,
обично
се
употребљава
релативни
координатни
64
систем
.
Када
је
покретачка
сила
реактивна
сила
млаза
,
употребљава
се
апсолутни
координатни
систем
.
У
случајевима
променљивог
кретања
(
убрзаног
или
успореног
)
употребљаваће
се
општи
облик
једначине
о
промени
количине
кретања
.
Активна
покретачка
сила
Брзина
и
проток
који
се
налазе
у
једначини
о
промени
количине
кретања
у
релативном
систему
имају
вредности
које
се
разликују
од
апсолутних
.
Брзина
у
релативном
координатном
систему
,
тј
.
брзина
преко
или
кроз
покретни
елемент
,
векторска
је
разлика
апсолутне
и
преносне
брзине
,
и
представља
релативну
брзину
.
Проток
преко
покретног
елемента
назива
се
редукованим
.
У
случајевима
великог
броја
лопатица
које
се
крећу
под
дејством
млаза
(
нпр
.
турбинско
коло
које
се
обрће
великом
брзином
,
или
има
велик
број
лопатица
),
редукован
проток
треба
заменити
стварним
,
апсолутним
протоком
,
најчешће
обележеним
са
Q
0
,
јер
у
сваком
и
најмањем
тренутку
времена
млаз
флуида
удара
у
нову
лопатицу
,
а
да
је
при
томе
не
сустиже
.
Међутим
,
без
обзира
на
проток
,
обавезно
се
у
једначину
о
промени
количине
кретања
уноси
релативна
брзина
,
уколико
се
ради
о
релативном
систему
,
v
=
v
0
-
u
,
јер
лопатица
у
сваком
кратком
временско
интервалу
прима
количину
флуида
брзином
v
0
,
а
даље
је
преноси
брзином
u
.
Реактивна
покретачка
сила
Реактивна
покретачка
сила
одређује
се
у
апсолутном
непокретном
координатном
систему
преко
вредности
апсолутног
протока
и
апсолутне
брзине
млаза
,
јер
покретан
систем
није
или
не
мора
да
буде
инерцијски
.
Проток
је
изражен
као
Q=A
⋅
v
rel
,
где
је
А
део
површине
омотача
контолне
запремине
кроз
коју
истиче
(
утиче
)
флуид
релативном
брзином
v
rel
.

66
Момент
количине
кретања
локално
је
паралелан
вектору
[ ]
,
r v
.
Струјање
у
радијалним
колима
Промена
момента
количине
кретања
користи
се
при
анализи
рада
турбомашина
.
Турбина
одузима
енергију
,
а
пумпа
предаје
енергију
флуидној
струји
уз
помоћ
роторског
кола
чије
се
лопатице
померају
само
у
тангенцијалном
правцу
.
Због
тога
је
користан
рад
последица
деловања
тангенцијалних
компоненти
сила
на
ротор
.
Радијалне
силе
не
могу
да
померају
ротор
у
радијалном
правцу
,
пошто
је
чврсто
причвршћен
на
осовину
,
па
не
производе
рад
.
Снага
која
се
преноси
са
или
на
осовину
,
добија
се
када
се
резултујући
момент
тангенцијалних
компонената
сила
помножи
са
угаоном
брзином
обртања
осовине
P=M
ω
.
Резултујући
момент
тангенцијалних
сила
,
пошто
је
проток
флуида
кроз
роторско
коло
константан
,
одређен
је
изразом
(
)
2
1
2 2
2
1 1
1
cos
cos
M
M
M
Q R v
R v
ρ
α
α
=
−
=
−
где
се
индекси
1
и
2
односе
на
излаз
и
улаз
флуида
.
Карактеристични
троуглови
брзине
за
улазне
и
излазне
пресеке
центрифугалне
пумпе
и
турбине
дати
су
на
слици
14
.
Из
претходне
једначине
види
се
да
пумпу
треба
пројектовати
тако
да
је
М
1
=0,
а
турбину
да
је
М
2
=0.
Тада
је
степен
ефикасности
η
највећи
,
а
пошто
је
M
P
QY
ω
ρ
=
=
напор
пумпе
је
2 2
2
cos
P
Y
u v
α
=
а
напор
турбине
1 1
1
cos
T
Y
u v
α
=
.
67
Слика
14.
Карактеристични
троуглови
брзина
за
роторско
коло
пумпе
(
лево
)
и
турбине
(
десно
)
v –
апсолутна
брзина
, u –
обимна
брзина
, v
rel
–
релативна
брзина
, v
2u
–
тангенцијална
компонента
апсолутне
брзине
, v
2r
–
радијална
компонента
апсолутне
брзине
Турбомашине
класификује
се
према
врсти
струјања
на
:
радијалне
,
аксијалне
и
радијално
-
аксијалне
.
На
слици
15
приказане
су
скице
типичних
турбомашина
.
Снага
турбине
је
:
(
)
(
)
(
)
1 1
1
2 2
2
1
1
2
2
1
1
2
2
cos
cos
T
T
u
u
u
u
P
M
m u v
u v
m R v
R
v
m u v
u v
ω
α
α
ω
ω
=
=
−
=
−
=
−
.
Специфичан
рад
кола
струјне
машине
дефинише
се
као
:
1
1
2
2
T
u
u
P
u v
u v
m
=
−
.
Ово
је
Ојлерова
једначина
за
турбине
.
Важи
и
за
пумпе
,
када
се
промене
предзнаци
.
Специфичан
рад
(
напор
)
је
рад
који
флуид
предаје
колу
.

69
Струјање
компресибилног
флуида
Једначина
континуитета
,
брзина
слабог
еластичног
поремећаја
,
Махов
конус
.
Уколико
у
флуидном
систему
долази
до
значајне
промене
густине
или
запремине
,
као
најважнији
струјни
ефект
,
јавља
се
компресибилност
.
Ово
се
дешава
при
наглим
убрзањима
,
појави
еластичних
таласа
,
брзинама
које
достижу
и
премашују
брзину
звука
,
струјање
у
разређеним
срединама
(
велике
висине
).
Овде
се
разматрају
само
реверзибилне
промене
–
компресибилно
струјање
без
губитака
представљено
изентропским
законом
.
Једначина
континуитета
за
компресибилан
флуид
је
:
1 1
1
2 2
2
v A
v A
m
ρ
ρ
=
=
где
се
индекси
1
и
2
односе
се
на
густину
,
брзину
и
површину
проточног
пресека
нормалну
на
правац
брзине
у
пресеку
1
и
2,
или
(
)
d
d
m
vA
ρ
=
.
Брзина
слабог
еластичног
поремећаја
–
брзина
звука
Брзина
звука
дефинише
се
као
брзина
простирања
малих
поремећаја
величина
стања
(=
звука
)
у
мирном
флуиду
.
То
је
брзина
сигнала
коју
треба
разликовати
од
брзине
струјања
самог
флуида
у
којем
се
сигнал
преноси
.
Слика
1.
Ударно
-
таласна
цев
Разматра
се
напредовање
таласног
фронта
у
каналу
константног
струјног
пресека
,
тзв
.
ударно
-
таласну
цев
(
слика
1
).
Према
излазу
цев
је
мембраном
подељена
на
два
дела
.
У
десном
делу
влада
ниски
притисак
,
70
а
у
левом
висок
.
Када
се
мембрана
уклони
,
долази
до
згушњавања
у
области
ниског
притиска
и
разређивања
у
области
високог
притиска
.
Како
су
у
питању
мали
поремећаји
,
сигнали
се
преносе
брзином
звука
c
(
слика
2
).
Разматра
се
околина
таласног
фронта
који
се
креће
удесно
.
Слика
2.
Деснокретни
талас
сабијања
То
је
нестационарни
процес
који
се
променом
брзине
–
с
може
разматрати
као
стационаран
. (
Стационаран
приказ
добија
се
када
се
посматрач
налази
на
фронту
таласа
.
У
сусрет
,
с
десна
на
лево
,
стиже
флуид
брзином
–
с
,
а
ниструјно
је
–
с
+d
v
).
Користе
се
основне
једначине
струјне
теорије
,
али
линеаризоване
.
Једначина
континуитета
је
:
(
)(
)
.
d
d
m
const
cA
c
v A
ρ
ρ
ρ
=
−
=
+
− +
одакле
,
уз
линеаризацију
,
следи
:
d
d
v
c
ρ
ρ
=
Применом
Бернулијеве
једначине
за
компресибилно
струјање
2
d
.
2
v
p
const
ρ
+
=
∫
добија
се
( )
(
)
2
2
d
d
d
d
2
2
p
p
p
c
c
v
p
p
ρ
ρ
+
−
− +
+
=
+
∫
∫
d
2
2
d
d
d
2
2
p
p
p
c
p
c
p
c v
ρ
ρ
+
+
=
−
+
∫
∫
.
Последњи
члан
једнакости
може
да
се
трансформише
на
следећи
начин
:
d
d
d(
d )
d
d(d )
d
d
p
p
p
p
p
p
p
p
p
p
p
p
p
ρ
ρ
ρ
ρ
ρ
ρ
+
+
=
=
+
=
+
∫
∫
∫
∫
∫
.

72
пораста
брзине
кретања
извора
успорава
се
ширење
поремећајног
таласа
уструјно
од
извора
поремећаја
,
а
убрзава
ниструјно
.
При
брзинама
кретања
извора
v
>
c
талас
поремећаја
простире
се
у
виду
конуса
(
Махов
конус
)
чија
је
карактеристика
угао
α
,
1
sin
c
v
Ma
α
= =
.
Слика
3.
Формирање
Маховог
конуса
Слика
формирања
Маховог
конуса
идентична
је
када
је
извор
поремећаја
непокретан
,
а
брзина
флуида
се
повећава
од
нуле
до
v
>
c
.
Основни
закони
компресибилног
струјања
Основни
закони
компресибилног
флуида
могу
се
добити
из
Бернулијевог
интеграла
Ојлерове
једначине
за
баротропан
флуид
при
занемареним
спољашњим
силама
:
2
d
.
2
v
p
const
ρ
+
=
∫
73
За
изотермско
струјање
,
p/
ρ
=const.
следи
2
2
2
1
1
2
ln
2
v
v
RT
ρ
ρ
−
=
или
2
2
2
1
1
1
1
2
ln
2
v
v
p
p
p
ρ
−
=
где
се
индекси
односе
на
две
тачке
исте
струјнице
.
За
изентропску
струју
,
p/
ρ
κ
=
const.
интеграцијом
једначине
2
d
.
2
v
p
const
ρ
+
=
∫
може
се
добити
један
од
следећих
израза
:
2
2
1
1
2
2
1
2
2
2
2
2
1
1
2
2
2
2
1
2
1
2
1
2
2
2
1
1
2
1
1
1
2
2
2
1
1
2
1
1
2
1
2
1
2
1
2
1
2
1
2
1
1
2
2
1
1
.
2
2
1
v
p
v
p
v
c
v
c
v
v
RT
RT
v
v
p
p
p
v
v
p
κ
κ
κ
κ
κ
κ
ρ
κ
ρ
κ
κ
κ
κ
κ
κ
κ
κ
ρ
ρ
κ
κ
ρ
ρ
−
−
+
=
+
−
−
+
=
+
−
−
+
=
+
−
−
−
=
−
−
−
=
−
−
За
случај
да
је
тачка
2 (
слика
4
)
зауставна
(
v
2
=0;
p
2
=
p
t
)
добијају
се
следеће
релације
у
којима
фигуришу
тоталне
величине
стања
:
Слика
4.
Зауставна
струјница

75
1
1
2
1
1
1
1
2
t
Ma
κ
ρ
κ
ρ
−
−
=
+
2
1
1
1
1
2
t
T
Ma
T
κ
−
= +
.
Бернулијева
једначина
показује
да
је
брзина
струје
мања
у
тачкама
где
влада
већа
брзина
звука
,
и
обрнуто
.
2
2
.
2
1
v
c
const
κ
+
=
−
Максимална
брзина
звука
c
max
припада
тачкама
где
флуид
мирује
,
тј
.
тачкама
v
=0.
t
max
t
t
p
c
RT
κ
κ
ρ
=
=
За
миран
стандардни
ваздух
на
нивоу
мора
, 15
0
С
,
c
max
=341 m/s.
Значајн
e
тачке
су
и
оне
где
се
струјна
брзина
изједначује
са
брзином
звука
. To
су
критичне
тачке
,
а
одговарају
им
критичне
брзине
v
kr
и
c
kr
.
Из
једначине
2
2
2
2
1
1
max
c
v
c
κ
κ
+
=
−
−
за
kr
kr
v
c
v
c
= =
=
добија
се
2
1
kr
kr
max
v
c
c
κ
=
=
+
.
За
ваздух
,
при
стандардним
условима
2
m
341 311
1, 408 1
s
kr
kr
v
c
=
=
=
+
Одређивање
критеријума
нестишљивог
струјања
гаса
Поставља
се
питање
до
ког
Маховог
броја
,
односно
брзине
,
може
неко
струјање
са
довољном
тачношћу
да
се
разматра
као
нестишљиво
.
Нека
се
претпостави
да
је
у
таквом
струјању
релативна
промена
густине
увек
мања
од
1 %.
76
2
1
2
1
2
1
1
1
...
2
1
1
...
1
2
2
t
Ma
Ma
Ma
κ
ρ
ρ
κ
−
=
=
= −
+
−
+
+
+
одакле
се
налази
2
... 0, 01
0,14
2
t
t
Ma
Ma
ρ ρ
ρ
−
=
+ ≤
⇒
≤
.
За
ваздух
на
собној
температури
,
гранична
брзина
струјања
нестишњивог
флуида
је
v
<50 m/s=180 km/h.
Истицање
кроз
конвергентан
млазник
Уобичајено
је
да
се
разматра
адијабатска
експанзија
(
истицање
)
из
великог
резервоара
у
коме
се
брзина
флуида
може
занемарити
(
v
1
=0)
(
слика
6
).
Слика
6.
Иситацање
гаса
из
резервоара
Карактеристичан
случај
је
када
на
изласку
из
млазника
максимална
брзина
струјања
достиже
брзину
звука
.
За
ову
брзину
потребан
однос
притиска
у
резервоару
и
на
излазном
пресеку
млазника
одређује
се
из
:
2
2
1
1
2
2
1
2
2
1
2
1
v
p
v
p
κ
κ
κ
ρ
κ
ρ
+
=
+
−
−
,
2
2
2
2
2
2
p
v
c
κ
ρ
=
=
као
1
2
2
1
2
2
1
1
1
1
2
1
2(
1)
p
p
p
κ
κ
κ
κ
κ
ρ
ρ
κ
ρ
κ
+
=
+
=
−
−
−
тј
.

78
2
1
2
1
kr
T
T
κ
=
+
.
У
овом
случају
постоји
скок
притиска
и
температуре
на
излазном
пресеку
млазника
са
2
kr
p
на
2
p
′
и
са
2
kr
T
на
2
T
′
.
Масени
проток
је
одређује
се
према
познатој
једначини
2 2
2
m
v A
ρ
=
.
За
случај
2
2
1
1
kr
p
p
p
p
′
<
када
је
2
2
max
v
c
v
=
=
,
масени
проток
износи
1
1
1
1
2
2
2
1
1
m
RT
A
κ
ρ
κ
κ
κ
−
=
+
+
.
За
случај
2
2
1
1
kr
p
p
p
p
′
>
тј
.
2
2
v
c
<
,
масени
проток
је
1
1
2
2
2
1
1
2
1
1
2
1
1
p
p
m
c
A
p
p
κ
κ
κ
ρ
κ
−
′
′
=
−
−
.
Струјање
кроз
сужени
пресек
Ако
компресибилни
флуид
струји
кроз
цевовод
променљивог
пресека
,
када
се
не
може
занемарити
квадрат
брзине
у
пресеку
1 (
слика
7
),
проток
се
може
одредити
решавањем
једначина
:
1
2
2
2
1
1
2
1
1
1
2
2
1
v
v
p
p
p
κ
κ
κ
κ
ρ
−
−
=
−
−
1 1
1
2 2
2
m
v A
v A
ρ
ρ
=
=
1
2
2
1
1
p
p
κ
ρ
ρ
=
.
79
Слика
7.
Компресибилно
струјање
кроз
сужени
пресек
За
резултат
се
добија
:
2
1
2
2
2
1
1
2
2
2
1
1
2
2
1
1
2
1
1
A
p
p
m
p
A
p
p
p
A
p
A
κ
κ
κ
κ
κ
ρ
κ
+
=
−
−
−
.
Струјање
кроз
конвергентно
–
дивергентан
млазник
Слика
8.
Струјање
компресибилног
флуида
кроз
конвергентно
-
дивергентан
млазник

81
Закон
промене
брзине
зависно
од
променљиве
површине
проточног
пресека
Из
једначине
континуитета
за
стационарно
компресибилно
струјање
следи
.
vA
const
ρ
=
log
log
log
log
v
A
C
ρ
+
+
=
d
d
d
A
v
A
v
ρ
ρ
= −
−
(*)
Користећи
се
једном
од
основних
Бернулијевих
једначина
,
2
.
2
1
v
p
const
κ
κ
ρ
+
=
−
диференцирањем
се
успоставља
зависност
:
2
d
d
d
0
1
p
p
v v
κ
ρ
κ
ρ
ρ
+
−
=
−
која
у
комбинацији
са
:
2
2
d
d
p
p
c
c
p
ρ
κ
κ
ρ
ρ
=
=
⇒
=
доводи
до
:
2
2
d
1
d
d
0
1
v
c
p
v
v
p
ρ
κ
ρ
+
−
=
−
.
**)
Из
једначине
изентропске
промене
.
p
const
κ
ρ
=
,
логаритмовањем
и
диференцирањем
следи
d
d
p
p
ρ
κ
ρ
=
.
Уврштавањем
у
једначину
(**)
добија
се
:
2
d
d
v
Ma
v
ρ
ρ
= −
Заменом
ове
једначине
у
(*),
добија
се
:
2
d
d
(1
)
A
v
Ma
A
v
= − −
.
За
струјање
некомпресибилног
флуида
,
једначина
континуитета
,
уз
трансформације
постаје
:
.
vA
const
=
82
log
log
log
v
A
C
+
=
d
d
A
v
A
v
= −
(***)
Поређење
једначина
за
d
A/A
за
компресибилно
(*)
и
некомпресибилно
(***)
струјање
,
указује
да
се
знатне
промене
проточног
пресека
компресибилног
струјања
,
у
односу
на
некомпресибилно
,
јављају
са
порастом
Маховог
броја
.
Промена
притиска
,
густине
и
брзине
приказана
је
на
слици
9
.
За
дозвучно
струјање
Ма
<1,
следи
:
o
притисак
и
густина
повећавају
се
са
повећањем
пресечне
површине
o
брзина
опада
са
повећањем
пресечне
површине
.
За
надзвучно
струјање
Ма
>1,
следи
:
o
притисак
и
густина
опадају
са
повећањем
пресечне
површине
o
брзина
расте
са
повећањем
пресечне
површине
.
Када
се
достигне
брзина
звука
,
промена
притиска
и
густине
веома
је
велика
чак
и
за
врло
мале
промене
пресека
.
У
близини
Ма
=1,
промене
у
брзини
и
густини
се
компензују
,
што
се
види
из
једначине
2
d
d
v
Ma
v
ρ
ρ
= −
.
Слика
9.
Промена
притиска
,
густине
и
брзине
у
изентропској
струји

84
Динамика
вискозног
флуида
Слика
1.
Струјни
режими
(
лево
–
ламинаран
;
десно
–
турбулентан
)
Динамика
вискозних
флуида
описује
кретање
реалних
течности
и
ређе
,
гасова
при
малим
брзинама
.
Ојлерова
једначина
за
кретање
идеалног
флуида
проширује
се
члановима
који
обухватају
утицај
унутрашњег
трења
флуидних
слојева
који
се
крећу
различитим
брзинама
.
Разликују
се
два
струјна
режима
:
ламинаран
и
турбулентан
(
слика
1
).
Ламинарно
струјање
карактерише
слојевито
кретање
флуидних
делића
без
њиховог
прелажења
из
слоја
у
слој
,
а
турбулентно
–
интензивно
,
хаотично
и
пулзацијско
мешање
флуидних
делића
различитих
слојева
.
Оба
струјна
режима
за
одржавање
својих
карактеристичних
профила
брзине
одузимају
од
флуидне
струје
енергију
(
енергија
притиска
),
којом
се
савлађује
унутрашње
трење
између
флуидних
слојева
,
односно
флуида
и
чврстих
граничних
површина
.
Ламинарно
струјање
може
да
буде
потенцијално
и
вртложно
,
док
је
турбулентно
увек
вртложно
.
Потенцијално
струјање
подразумева
постојање
профила
брзине
у
струјном
пољу
без
вртлога
(
елементарних
или
макроскопских
).
У
вртложном
ламинарном
струјању
најчешће
постоји
правилан
стабилан
низ
елементарних
вртлога
између
два
слоја
флуидних
делића
који
се
крећу
различитим
брзинама
или
макроскопски
вртлози
.
Постојање
оба
струјна
режима
Рејнолдс
је
доказао
једноставним
огледом
1883.
године
(
слика
2
).
У
хоризонталну
стаклену
цев
која
је
једним
крајем
постављена
у
суд
са
водом
,
уноси
се
танак
млаз
обојене
течности
.
Оба
тока
могу
да
се
регулишу
одговарајућим
славинама
.
За
мале
брзине
струјања
нит
обојене
течности
не
меша
се
са
водом
по
целој
дужини
стаклене
цеви
,
док
се
за
веће
брзине
интензивно
мешање
дешава
већ
на
краћем
растојању
од
улаза
у
стаклену
цев
.
85
Слика
2.
Рејнолдсов
оглед
Рејнолдс
је
пронашао
да
су
утицајне
величине
од
којих
зависи
врста
струјног
режима
:
v
,
d
,
ρ
,
η
,
односно
бездимензијски
критеријум
vd/
ν
,
који
је
по
њему
добио
име
Рејнолдсов
број
(
Re
).
У
уобичајеним
условима
,
када
се
не
поклања
изузетна
пажња
експерименту
,
прелазак
из
ламинарног
у
турбулентно
струјање
дешава
се
између
Re
=2000 -
4000.
Екстремно
добијене
границе
припротицању
кроз
цев
су
: 2000
за
прелазак
из
турбулентног
у
ламинарно
и
40.000
из
ламинарног
у
турбулентно
струјање
.
Вредност
Re
=2300
често
се
употребљава
као
граница
између
ламинарног
и
турбулентног
струјања
у
цевима
кружног
пресека
.
Гранични
Re
број
зависи
од
врсте
струјања
;
нпр
.
o
за
струјање
око
лопте
Re=vd/
ν
=
0,1 (
d
је
пречник
лопте
),
o
за
равну
плочу
Re=vl/
ν
=
500.000 (
l
је
дужина
у
правцу
струјања
).
При
ламинарном
струјању
губитак
енергије
пропорционалан
је
средњој
брзини
,
а
за
турбулентно
струјање
пропорционалан
је
v
1,7÷2,0
.
Носиоци
утрошка
флуидне
енергије
су
елементарни
вртлози
када
је
у
питању
трење
при
струјању
кроз
цевоводе
и
преко
равних
танких
површина
.
Међутим
,
при
опструјавању
тела
,
када
долази
до
отцепљења
флуидне
струје
,
јављају
се
вртлози
,
многоструко
већи
и
невезани
са
елементарним
вртлозима
,
на
чије
се
одржавање
троше
енергија
флуидне
струје
или
енергија
тела
које
се
креће
кроз
миран
флуид
.
Овим
питањем
бави
се
поглавље
о
отпорима
.

87
2
2
2
2
2
2
2
3
y
y
x
x
x
x
x
z
z
v
v
v
v
v
v
v
v
v
x
y
z
x
x
y
z
x
x
y
z
ν
ν
ν
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
+
+
+
+
+
−
+
+
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
што
доводи
до
:
1
div
3
x
v
v
x
ν
ν
∂
∆ +
∂
.
Слика
3.
Деловање
вискозних
сила
на
издвојени
флуидни
делић
За
стационарно
струјање
нестишљивог
флуида
важи
:
div
0
v
=
,
па
Навије
-
Стоксова
једначина
изгледа
овако
у
:
векторском
d
1
grad
d
v
f
p
v
t
ν
ρ
= −
+ ∆
,
односно
у
скаларном
облику
d
1
d
d
1
d
d
1
,
d
x
x
x
y
y
y
z
z
z
v
p
f
v
t
x
v
p
f
v
t
y
v
p
f
v
t
z
ν
ρ
ν
ρ
ν
ρ
∂
=
−
+ ∆
∂
∂
=
−
+ ∆
∂
∂
=
−
+ ∆
∂
где
су
на
левој
страни
једначине
инерцијске
силе
,
а
редом
на
десној
страни
:
спољашње
силе
,
силе
притиска
и
вискозне
силе
.
Развијени
облик
Навије
-
Стоксове
једначине
је
:
88
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
1
1
1
x
x
x
x
x
x
x
y
z
x
y
y
y
y
y
y
x
y
z
y
z
z
z
z
z
z
x
y
z
z
v
v
v
v
v
v
p
v
v
v
f
x
y
z
x
x
y
z
v
v
v
v
v
v
p
v
v
v
f
x
y
z
y
x
y
z
v
v
v
v
v
v
p
v
v
v
f
x
y
z
z
x
y
z
ν
ρ
ν
ρ
ν
ρ
∂
∂
∂
∂
∂
∂
∂
+
+
=
−
+
+
+
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
+
+
=
−
+
+
+
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
+
+
=
−
+
+
+
∂
∂
∂
∂
∂
∂
∂
са
пратећом
једначином
континуитета
:
0
y
x
z
v
v
v
x
y
z
∂
∂
∂
+
+
=
∂
∂
∂
.
При
ламинарном
струјању
нестишљивог
флуида
Навије
-
Стоксове
једначине
разликују
се
од
Ојлерове
за
члан
v
∆
.
Како
из
rot rot
grad div
v
v
v
=
− ∆
а
при
потенцијалном
струјању
важи
grad , rot
0,
div
0
v
v
v
φ
=
=
=
;
следи
да
је
:
0
v
∆ =
,
тј
.
да
велика
група
проблема
из
струјања
,
како
вискозног
тако
и
идеалног
флуида
има
заједничко
решење
,
јер
се
описује
истим
једначинама
.
Међутим
,
због
различитих
граничних
услова
и
решења
су
сасвим
различита
.
До
израза
за
распоред
притисака
и
брзине
може
се
доћи
постављањем
услова
динамичке
равнотеже
за
издвојену
елементарну
флуидну
запремину
.
Ову
анализу
погодно
је
применити
када
инерцијске
силе
могу
да
се
занемаре
.
Тангенцијални
напон
који
одређује
силу
трења
налази
се
посредством
Њутновог
закона
.
Ламинарно
струјање
–
Хаген
-
Поазејево
струјање
кроз
кружну
цев
Посматра
се
део
хоризонталне
цеви
.
Нека
је
струјање
развијено
(
нема
промене
профила
брзине
у
правцу
струјања
-
х
осе
).
Тиме
се
утврђује
да
је
разматрани
део
цеви
довољно
далеко
од
њеног
улаза
.
При
слојевитом
струјању
флуида
кроз
цев
,
притисак
је
константан
по
целом
струјном
пресеку
.
За
одређивање
профила
брзине
примениће
се
једначина
о
промени
количине
кретања
(
ЈПКК
)
на
коаксијални

90
2
4
max
p R
v
l
η
∆
=
.
Слика
5.
Профил
брзине
и
апсолутне
вредности
тангенцијалног
напона
при
ламинарном
струјању
кроз
цеви
Интеграцијом
профила
брзине
по
струјном
пресеку
добија
се
проток
који
је
једнак
производу
средње
брзине
и
површине
струјног
пресека
:
2
2
2
0
d
1
2
d
2
2
R
max
max
sr
max
r
v
v
r
Q
v A
v A
v
r r
R
A
R
π
π
=
=
=
=
−
=
=
∫
∫
.
Одавде
се
лако
закључује
да
важи
2
max
sr
v
v
=
.
Сада
за
запремински
проток
следи
једнакост
:
4
2
8
max
sr
v
pR
Q
v A
A
l
π
η
∆
=
=
=
одакле
се
види
да
је
проток
пропорционалан
паду
притиска
и
четвртом
степену
полупречника
цеви
:
Q
∼∆
p
,
Q
∼
R
4
.
Ови
закључци
позанти
су
као
Хаген
-
Поазејев
закон
.
Он
показује
карактеристичну
зависност
запреминског
протока
.
Нарочито
је
важно
да
се
проток
мења
као
четврти
степен
полупречника
цеви
.
Ово
указује
да
умањење
полупречника
води
драстичном
смањењу
протока
(
примена
у
медицини
).
Поставља
се
питање
колики
је
пад
притиска
(
губитак
притиска
)
у
цеви
кружног
попречног
пресека
(
кружној
цеви
)
при
задатом
запреминском
протоку
.
91
У
паду
притиска
испољава
се
,
у
правцу
струјања
,
утицај
трења
.
При
томе
профил
брзине
остаје
непромењен
.
Из
2
4
max
p R
v
l
η
∆
=
следи
2
2
4
8
max
sr
lv
vlv
p
R
R
η
ρ
∆ =
=
Имајући
на
уму
следеће
:
64
lam
Re
λ
=
и
sr
v D
Re
ν
=
добија
се
:
2
2
sr
lam
v
l
p
D
ρ
λ
∆ =
.
Први
чинилац
на
десној
страни
је
коефицијент
губитака
који
узима
у
обзир
трење
у
цевима
(
коефицијент
трења
),
други
чинилац
карактерише
геометрију
,
а
трећи
динамички
притисак
.
Пад
притиска
је
линеарна
функција
дужине
цеви
.
Турбул
e
нтно
струјање
Снимања
покретним
камерама
и
великим
бројем
снимака
у
секунди
потврдила
су
још
у
19.
веку
постављене
законитости
о
пулсирајућем
карактеру
струјних
карактеристика
турбулентног
струјања
брзине
,
притиска
и
напона
.
Тренутна
вредност
брзине
v
(
x,y,z
)
може
да
се
растави
на
два
сабирка
,
а
то
су
:
o
временски
o
средњена
вредност
–
просечна
вредност
( , , )
v x y z
o
брзина
тренутног
одступања
(
одступање
од
просечне
-
осредњене
вредности
,
брзина
пулсације
,
флуктуације
,
осциловања
)
( , , , )
v x y z t
′
.
Тако
следи
:

93
0
1
d
T
p
p t
T
=
∫
и
0
1
d
T
t
T
τ
τ
=
∫
.
Тренутна
одступања
притиска
и
тангенцијалних
напона
су
:
p
p
p
τ
τ τ
′
′
= −
= −
а
њихове
просечне
вредности
,
као
што
је
већ
речено
:
0
0
1
1
d
0
d
0
T
T
p
p t
t
T
T
τ
τ
′
′
′
′
=
=
=
=
∫
∫
.
Увођењем
просечних
(
по
времену
осредњених
)
параметара
турбулентног
струјања
могу
се
решити
многи
практични
задаци
,
међутим
,
није
добијена
довољна
представа
о
унутрашњој
структури
кретања
течности
,
па
се
уводе
и
други
параметри
од
којих
су
три
основна
:
1.
Интензитет
турбуленције
–
бездимензијска
величина
У
бројиоцу
се
налази
корен
средњег
квадратног
одступања
као
карактеристична
мера
,
у
општем
случају
,
ма
које
осцилаторне
величине
.
Она
је
овде
дефинисана
у
односу
на
средњу
брзину
струјања
на
посматраном
месту
.
Интензитет
турбуленције
заједно
са
фреквенцијом
промене
знака
представља
кванититативну
меру
турбуленције
.
2
v
I
v
′
=
.
Слика
7.
Корен
средњег
квадратног
одступања
као
показатељ
јачине
турбуленције
На
слици
7
приказана
су
два
записа
брзине
истих
просечних
вредности
,
али
запис
са
леве
стране
има
већи
интензитет
турбуленције
,
јер
је
корен
квадратног
одступања
2
v
′
већи
.
2
v
′
2
v
′
94
2.
Коефицијент
корелације
(
параметар
турбулентности
)
Одређује
степен
повезаности
брзина
тренутних
одступања
.
На
пример
1
x
v
′
и
2
x
v
′
су
пројекције
осцилација
на
осу
х
у
тачкама
М
1
и
М
2
на
слици
8
.
( ) ( )
_______
1
2
1 2
_______ _______
2
2
1
2
x
x
x
x
v v
R
v
v
−
′ ′
=
′
′
.
Слика
8.
Две
трајекторије
непредвидивог
кретања
делића
који
стижу
у
различитим
временским
тренуцима
у
исту
тачку
пресека
II
3.
Учестаност
осцилација
(
брзине
)
Експерименти
су
показали
да
је
у
турбулентним
струјањима
присутан
широк
спектар
фреквенција
осцилација
.
Анализа
турбулентних
појава
показала
је
да
треба
разликовати
две
врсте
осцилација
:
макроскопске
и
микроскопске
.
Макроскопске
осцилације
су
са
знатним
амплитудама
и
релативно
малим
фреквенцијама
.
Микроскопске
осцилације
су
са
малим
амплитудама
и
великим
фреквенцијама
.
Макроскопске
осцилације
генерише
врста
струјања
,
њени
почетни
и
гранични
услови
и
карактер
локалних
поремећаја
.
Микроскопске
осцилације
су
универзалне
–
опште
присутне
и
играју
важну
улогу
при
дисипацији
енергије
у
струји
вискозног
флуида
.
Ова
микроскопска
турбуленција
,
компонована
од
великог
броја
малих
вртлога
,
брзо
претвара
механичку
енергију
у
топлотну
.
Примери
макроскопских
вртлога
или
вртлога
великих
размера
су
речни
вртлози
и
“
ваздушне
пијавице
”.
Велики
вртлози
генеришу
вртложење
малих
размера
(
елементарни
вртлози
,
микроскопске
пулсације
).

96
И
даље
,
(
)
(
)
(
)
_____
_____
_____
1
1
x
x
x
y
x z
x
x
x
x
x
y
x
z
v v
v v
v v
x
y
z
p
f
v
v v
v v
v v
x
x
y
z
ν
ρ
ρ
ρ
ρ
ρ
∂
∂
∂
+
+
=
∂
∂
∂
∂
∂
∂
∂
′ ′
′ ′
′ ′
−
+ ∆ −
+
+
∂
∂
∂
∂
пошто
је
0
0
0
0
0
0
0
1
0, a
због
d
0
и
d
0
1
d
d
1
d
d
0,
због
d
0
T
T
x
x
x
T
T
x
y
x
y
x
y
T
T
T
x
x
y
y
y
v
v
v t
t
t
T
t
v v
v v t
t v v
T
T
v
v v t
v t
v t
T
T
′
∂
∂
′
=
=
=
∂
∂
=
=
′
′
′
=
=
=
∫
∫
∫
∫
∫
∫
∫
_____
0
0
0
1
d
1
1
d
,
d
T
x
y
x
y
T
T
x
v v t
v v
T
f
t
f
p t
p
T
T
′ ′
′ ′
=
=
=
∫
∫
∫
Понављајући
поступак
за
преостале
две
пројекције
добијају
се
Рејнолдсове
једначине
за
стационарно
просечно
турбулентно
кретање
нестишљиве
течности
.
(
)
(
)
(
)
(
)
(
)
(
)
____
____
____
2
____
x x
x
y
x
z
x
x
x
x
x
x
y
x
z
y
x
y
y
y
z
y
y
y
y
x
v v
v v
v v
x
y
z
v
v
v
p
f
v
v v
v v
x
x
x
y
y
z
z
v v
v v
v v
x
y
z
v
v
p
f
v v
y
x
x
y
y
ρ
ρ
η
ρ
η
ρ
η
ρ
ρ
ρ
η
ρ
η
ρ
∂
∂
∂
+
+
=
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
′
′ ′
′ ′
−
+
−
+
−
+
−
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
+
+
=
∂
∂
∂
∂
∂
∂
∂
∂
′ ′
−
+
−
+
−
∂
∂
∂
∂
∂
(
)
(
)
(
)
____
____
2
____
____
____
2
y
y
y
z
z
x
z
y
z
z
z
z
z
z
z
x
z
y
z
v
v
v v
z
z
v v
v v
v v
x
y
z
v
v
v
p
f
v v
v v
v
z
x
x
y
y
z
z
η
ρ
ρ
ρ
η
ρ
η
ρ
η
ρ
∂
∂
′
′ ′
+
−
∂
∂
∂
∂
∂
+
+
=
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
′ ′
′ ′
′
−
+
−
+
−
+
−
∂
∂
∂
∂
∂
∂
∂
97
0
y
x
z
v
v
v
x
y
z
∂
∂
∂
+
+
=
∂
∂
∂
Чланови
:
____
____
____
2
2
2
,
,
y
x
z
x
y
z
v
v
v
v
v
v
x
y
z
η
ρ
η
ρ
η
ρ
∂
∂
∂
′
′
′
−
−
−
∂
∂
∂
представљају
нормалне
напоне
у
течности
,
проузроковане
вискозношћу
и
временски
осредњеним
стварним
брзинама
течности
.
А
чланови
____
____
____
____
____
____
,
,
,
,
y
x
x
x
y
x
z
y
x
y
z
z
y
z
z
x
z
y
v
v
v
v v
v v
v v
y
z
x
v
v
v
v v
v v
v v
z
x
y
η
ρ
η
ρ
η
ρ
η
ρ
η
ρ
η
ρ
∂
∂
∂
′ ′
′ ′
′ ′
−
−
−
∂
∂
∂
∂
∂
∂
′ ′
′ ′
′ ′
−
−
−
∂
∂
∂
су
тангенцијални
напони
услед
вискозности
и
осредњених
стварних
брзина
.
Чланови
____
____
____
2
2
2
,
,
x
y
z
v
v
v
ρ
ρ
ρ
′
′
′
називају
се
нормални
турбулентни
напони
,
а
____
____
____
,
,
x
y
x
z
y
z
v v
v v
v v
ρ
ρ
ρ
′ ′
′ ′
′ ′
тангенцијални
турбулентни
напони
.
За
нестационарно
,
временски
осредњено
турбулентно
кретање
нестишљиве
течности
могу
да
се
напишу
Навије
-
Стоксове
једначине
са
осредњеним
вредностима
брзине
и
притиска
уз
додатак
нормалних
и
тангенцијалних
турбулентних
напона
:
(
)
(
)
(
)
(
)
(
)
(
)
2
2
1
1
1
1
.
x
x
x
x
x
y
z
x
x
x
x
y
x
z
y
y
y
y
x
y
z
y
y
x
y
y
y
z
v
v
v
v
v
v
v
t
x
y
z
p
f
v
v
v v
v v
x
x
y
z
v
v
v
v
v
v
v
t
x
y
z
p
f
v
v v
v
v v
y
x
y
z
ν
ρ
ρ
ρ
ρ
ρ
ν
ρ
ρ
ρ
ρ
ρ
∂
∂
∂
∂
+
+
+
=
∂
∂
∂
∂
∂
∂
∂
∂
′
′ ′
′ ′
−
+ ∆ +
−
+
−
+
−
∂
∂
∂
∂
∂
∂
∂
∂
+
+
+
=
∂
∂
∂
∂
∂
∂
∂
∂
′ ′
′
′ ′
−
+ ∆ +
−
+
−
+
−
∂
∂
∂
∂

99
_____
j
i
i
j
v
v
v v
j
i
ρ
η
∂
∂
′ ′
′
= −
+
∂
∂
где
је
η
′
турбулентна
вискозност
и
разликује
се
од
динамичке
вискозности
η
.
Тангенцијални
напон
услед
вискозности
и
осредњавања
стварних
брзина
(
тј
.
турбулентности
)
добија
се
збиром
j
j
j
i
i
i
ij
p
v
v
v
v
v
v
j
i
j
i
j
i
τ
η
η
η
∂
∂
∂
∂
∂
∂
′
=
+
+
+
=
+
∂
∂
∂
∂
∂
∂
где
је
p
η
привидна
вискозност
дата
са
:
p
η
η η
′
= +
.
Структура
η
′
није
позната
.
Њена
вредност
мења
се
у
зависности
од
области
струјног
поља
.
Нпр
.
у
близини
зида
η
′
је
мала
у
односу
на
η
,
а
у
матици
струјања
,
обрнуто
,
η
′
је
неколико
пута
веће
од
η
,
тако
да
у
том
случају
важи
d
d
x
v
y
τ η
′
=
.
Увођењем
привидне
вискозности
тангенцијални
напон
задржава
свој
једноставан
и
логичан
облик
који
је
Њутн
дефинисао
.
Међутим
,
за
његово
коришћење
потребна
су
експериментална
мерења
,
па
му
је
значај
сужен
на
случајеве
за
које
постоје
мерни
подаци
.
Прикупљања
и
обрада
статистичких
података
основни
су
извор
информација
о
механизмима
турбуленције
,
усмерени
су
ка
проналажењу
уопштене
форме
турбулентних
напона
.
Од
многобројних
постављених
хипотеза
изложиће
се
најчешће
коришћена
Прантлова
хипотеза
.
Прантлова
хипотеза
(Prandtl, 1925)
Аналогно
кинетичкој
теорији
гасова
где
молекули
,
пре
него
што
се
сударе
са
другим
честицама
,
пређу
средње
растојање
,
уводи
се
Прантлова
путања
мешања
која
представља
растојање
коју
делић
флуида
пређе
у
турбулентној
средини
,
пре
него
што
се
измеша
са
околином
и
тако
јој
преда
своју
количину
кретања
(
слика
9
).
Према
Прантлу
,
један
делић
из
слоја
у
доспе
при
y
v
′
>0
у
ниво
y+l
1
.
Тако
он
поседује
мањак
брзине
у
односу
на
околину
,
чија
је
величина
1
1
d
( )
(
)
d
x
x
x
v
v y
v y
l
l
y
−
+
= −
.
100
Слика
9.
Прантлова
путања
мешања
Прантл
је
овај
мањак
брзине
схватио
као
осцилујућу
компоненту
брзине
(
брзину
одступања
)
у
нивоу
y+ l
1
,
тј
.
1
d
( )
d
x
x
v
v y
l
y
′
= −
.
По
аналогији
следи
и
:
2
d
( )
d
x
y
v
v
y
l
y
′
=
.
Сада
се
за
Рејнолдсов
привидни
тангенцијални
напон
добија
2
2
_____
___
2
1 2
d
d
d
d
x
x
x
y
v
v
v v
l l
l
y
y
τ
ρ
ρ
ρ
′ ′
= −
=
=
где
је
l
карактеристична
мера
дужине
мешања
у
турбулентном
струјању
,
једнака
Прантловој
путањи
мешања
.
Да
је
она
функција
од
у
мора
се
експериментално
показати
,
пошто
је
ова
теорија
изнесена
као
полуемпиријска
.
У
претходном
изразу
важна
је
зависност
од
квадрата
градијента
брзине
.
То
упућује
на
типичне
разлике
у
односу
на
ламинарно
струјање
.
Дакле
,
2
2
d
d
d
d
x
x
uk
v
v
l
y
y
τ
η
ρ
=
+
Прантлова
хипотеза
и
поред
бројних
оправданих
приговора
најчешће
се
користи
у
пракси
.
Један
од
приговора
је
да
би
у
близини
граничне
површине
и
на
њој
,
где
не
долази
до
промене
количине
кретања
,
путања
мешања
требало
да
тежи
нули
,
а
то
није
обухваћено
Прантловом
теоријом
.
Профил
брзине
и
пад
притиска
у
кружној
цеви

102
(
)
6
1
2 log
0,8
важи до
3 10
D
turb
D
turb
Re
Re
λ
λ
=
−
≈ ⋅
.
Никурадзе
је
направио
експеримент
,
користећи
се
вештачком
храпавошћу
(
просејан
песак
лаком
залепљен
за
зидове
цеви
).
Резултати
његових
испитивања
представљени
су
дијаграмом
(
слика
11
)
у
зависности
од
Рејнолдосвог
броја
и
R/e
(
R
је
полупречник
цеви
,
а
е
средња
висина
неравнина
).
Слика
11.
Никурадзеов
дијаграм
зависности
коефицијента
храпавости
од
Рејнолдсовог
броја
и
R/e
Из
Блазијусовог
закона
следи
интересантна
последица
за
просечну
брзину
турбулентног
струјања
.
Из
Блазијусове
формуле
узме
се
λ
turb
и
уврсти
у
израз
за
тангенцијални
напон
и
уоче
се
експоненти
:
7
1
1
7
1
1
2
4
4
4
4
4
4
2
4
turb
z
sr
sr
max
v
v
R
v
R
λ
ρ
τ
ρ ν
ρ
ν
−
−
=
≈
≈
.
Профил
брзине
,
за
практичне
потребе
,
може
да
се
представи
степеним
законом
са
независним
експонентом
m
:
( )
m
max
y
v y
v
R
=
где
је
у
нормално
растојање
од
зида
цеви
y=R-r
.
Ако
се
претходни
израз
реши
по
максималној
просечној
брзини
и
уврсти
у
израз
за
тангенцијални
напон
,
добија
се
103
7
7
1
7
1
4
4
4
4
4
2
m
m
z
v y
R
ρ
τ
ν
−
−
=
.
Прантл
и
Карман
изнели
су
хипотезу
да
тангенцијални
напон
на
зиду
цеви
не
би
требало
да
зависи
од
њеног
радијуса
,
тј
.
да
је
турбулентно
струјање
одређено
,
више
или
мање
,
локалним
подацима
струјног
поља
.
За
R
0
следи
m
=1/7.
1
1
7
7
( )
1
max
max
y
r
v y
v
v
R
R
=
=
−
Ово
је
важан
1/7
степени
закон
који
одговара
Блазијусовом
закону
.
Слика
12.
Ламинаран
и
турбулентан
профил
брзине
у
цеви
На
слици
12
приказани
су
ламинарни
и
турбулентни
профили
за
исти
запремински
проток
.
1.
За
m
=1/7
профил
има
два
карактеристична
места
.
У
близини
зида
постоји
бесконачан
успон
.
Несумњиво
је
да
је
у
непосредној
близини
зида
струјање
ламинарно
(
гранични
подслој
)
и
ту
не
може
да
се
користи
степени
закон
.
На
профилу
брзине
на
оси
цеви
настаје
превој
(
екстрем
).
2.
Са
порастом
Рејнолдсовог
броја
смањује
се
експонент
,
профил
брзине
постаје
све
више
правоугаоног
облика
.
За
ламинарно
и
турбулентно
струјање
кроз
храпаву
цев
важе
следеће
констатације
1.
За
ламинарно
струјање
је
λ
=
f
(
Re
),
тј
.
коефицијент
трења
не
зависи
од
храпавости
.
2.
За
турбулентно
струјање
важе
следеће
алтернативе

105
Слика
13.
Карактеристична
струјања
турбулентног
режима
106
Гранични
слој
Динамика
флуида
треба
да
одговори
на
питања
колики
је
минималан
утрошак
енергије
при
протицању
флуида
или
при
кретању
тела
кроз
флуид
и
на
који
начин
се
то
постиже
.
Врста
,
порекло
и
величина
отпорних
сила
,
од
којих
зависе
енергијска
улагања
,
анализирају
се
у
различитим
струјним
моделима
:
идеалном
флуиду
–
потенцијално
струјање
и
вискозном
флуиду
.
Ако
је
облик
тела
такав
да
се
при
кретању
кроз
флуид
не
изазива
цепање
струјница
,
теорија
потенцијаланог
струјања
и
Ојлерова
једначина
довољне
су
да
се
добије
стварна
слика
о
отпорима
.
За
кретање
кроз
вискозан
флуид
теорија
ламинарног
струјања
је
довољна
уколико
струјање
не
напушта
овај
струјни
режим
;
али
за
турбулентно
опструјавање
тела
,
које
увек
прати
већа
брзина
,
за
пунија
незаобљена
тела
,
отворен
систем
једначина
не
може
да
да
ни
приближан
одговор
о
квалитету
отпорних
сила
.
Пошто
су
у
природи
и
техници
присутна
углавном
оваква
струјања
,
њихова
математичка
анализа
са
пратећим
експерименталним
истраживањима
и
данас
представља
основну
делатност
бројних
истраживачких
центара
.
Прантл
је
1904.
године
,
постављајући
теорију
граничног
слоја
направио
највећи
корак
,
јер
је
успоставио
везу
између
кретања
идеалног
и
реалног
флуида
.
Та
веза
је
гранични
слој
,
област
флуида
у
непосредној
околини
тела
које
се
креће
кроз
флуид
,
у
којој
се
вискозни
утицаји
не
могу
занемарити
(
слике
1
и
2
).
Према
Прантлу
ефекат
унутрашњег
трења
за
флуиде
са
релативно
малом
вискозношћу
,
значајан
је
само
у
уској
области
која
окружује
границе
флуида
.
Према
овој
хипотези
струјање
изван
овог
,
уског
региона
,
може
се
сматрати
као
струјање
идеалног
флуида
,
тј
.
као
потенцијална
струја
.
У
поглављу
динамика
вискозног
флуида
речено
је
да
се
за
велике
Рејнолдсове
бројеве
вискозност
може
занемарити
у
односу
на
инерцијске
силе
и
обрнуто
.
Међутим
,
у
области
унутар
граничног
слоја
присутни
су
и
велики
Рејнолдсови
бројеви
(
инерцијске
силе
)
и
вискозни
ефекти
,
што
је
последица
граничних
услова
.
Гранични
услови
за
идеалну
течност
захтевају
само
одсуство
нормалне
компоненте
брзине
док
тангенцијална
компонента
на
површини
тела
,
око
кога
се
врши
оптицање
,
може
да
постоји
.
Међутим
,
у
вискозној
реалној

108
o
Распоред
брзине
и
притиска
практично
се
мало
разликује
од
њиховог
распореда
у
потенцијалном
струјању
невискозне
течности
;
вртложење
услед
вискозности
је
овде
слабо
.
За
практична
разматрања
у
области
изван
граничног
слоја
вискозне
силе
и
турбулентне
силе
трења
могу
се
занемарити
,
струјање
сматрати
потенцијалним
,
а
течност
невискозном
.
Под
дебљином
граничног
слоја
δ
подразумева
се
дебљина
флуида
у
близини
чврсте
површине
на
чијој
граници
је
брзина
скоро
једнака
брзини
основне
флуидне
струје
(
U
).
Уобичајено
одступање
је
1%.
Брзина
унутар
граничног
слоја
v
x
креће
се
од
нуле
до
99%
U
(
U
je
једнолика
брзина
изван
граничног
слој
a),
а
изван
граничног
слоја
v
x
=
U
,
како
је
представљено
на
слици
3
.
Слика
3.
Промена
брзине
унутар
и
изван
граничног
слоја
Слика
4.
Дебљина
изгубљеног
протока
За
анализу
граничног
слоја
потребно
је
познавати
још
две
геометријске
карактеристике
:
дебљину
изгубљеног
протока
δ
1
и
дебљину
изгубљене
количине
кретања
δ
2
.
Ове
дебљине
не
треба
мешати
са
дебљином
граничног
слоја
δ
.
Дебљина
изгубљеног
протока
δ
1
(
слика
4
)
је
замишљено
померање
непокретне
површине
према
основној
струји
,
толико
да
проток
флуида
кроз
ту
смањену
површину
попречног
пресека
,
са
једноликом
брзином
U
,
буде
једнак
стварном
протоку
.
Једнакост
осенчених
површина
на
слици
испуњава
овај
захтев
.
109
(
)
1
0
d
x
U
U
v
y
δ
δ =
−
∫
односно
1
0
1
d .
x
v
y
U
δ
δ
=
−
∫
Дебљина
изгубљене
количине
кретања
δ
2
(
слика
5
)
добија
се
изједначавањем
изгубљене
количине
кретања
кроз
гранични
слој
која
је
последица
опадања
брзине
v
x
<
U
(
у
односу
на
основну
струју
),
са
замишљеном
количином
кретања
у
струји
граничног
слоја
дебљине
δ
2
и
непроменљиве
брзине
U
.
Израз
за
δ
2
следи
из
(
)
2
2
0
d
x
x
U
U
v v y
δ
ρδ
ρ
=
−
∫
Односно
2
0
1
d
x
x
v
v
y
U
U
δ
δ
=
−
∫
.
Слика
5.
Дебљина
изгубљене
количине
кретања
На
слици
6
приказане
су
илустрације
различитих
струјних
режима
,
формирање
граничног
слоја
и
вртложног
трага
,
формирање
ламинарног
струјања
на
улазу
у
цев
.

111
Једначине
кретања
у
ламинарном
граничном
слоју
За
стационарно
,
раванско
струјање
,
уз
занемаривање
спољашњих
сила
,
Навије
-
Стоксове
једначине
са
пратећом
једначином
континуитета
имају
облик
:
2
2
2
2
2
2
2
2
1
1
x
x
x
x
x
y
y
y
y
y
x
y
v
v
v
v
p
v
v
x
y
x
x
y
v
v
v
v
p
v
v
x
y
y
x
y
ν
ρ
ν
ρ
∂
∂
∂
∂
∂
+
= −
+
+
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
+
= −
+
+
∂
∂
∂
∂
∂
0
y
x
v
v
x
y
∂
∂
+
=
∂
∂
.
Пошто
је
гранични
слој
танак
(
његова
дебљина
ће
се
касније
одредити
)
кретање
у
њему
ће
у
основи
бити
паралелно
са
равни
(
х
-
оса
)
око
које
се
врши
оптицање
,
па
је
брзина
v
y
мала
у
односу
на
v
x
.
У
правцу
y
-
осе
брзина
v
x
се
нагло
мења
,
док
се
у
правцу
х
-
осе
брзина
v
x
мења
споро
.
Из
наведеног
следи
да
се
може
занемарити
:
2
2
x
v
x
∂
∂
у
односу
на
2
2
x
v
y
∂
∂
и
2
2
,
,
y
y
y
x
y
v
v
v
v
v
x
y
y
∂
∂
∂
∂
∂
∂
у
односу
на
2
2
,
,
x
x
x
x
y
v
v
v
v
v
x
y
y
∂
∂
∂
∂
∂
∂
;
тако
да
Навије
-
Стоксове
једначине
постају
2
2
1
0
x
x
x
x
y
v
v
v
p
v
v
x
y
x
y
p
y
ν
ρ
∂
∂
∂
∂
+
= −
+
∂
∂
∂
∂
∂
=
∂
0
y
x
v
v
x
y
∂
∂
+
=
∂
∂
како
је
први
извео
Прантл
.
Једначина
0
p
y
∂
=
∂
показује
да
се
притисак
преноси
непромењен
кроз
гранични
слој
.
112
Притисак
p
сматра
се
познатим
јер
се
одређује
из
струјања
невискозног
флуида
у
близини
граничног
слоја
из
једначине
2
1
.
2
p
U
const
ρ
+
=
где
је
са
U
означена
брзина
у
тачкама
на
спољашњој
граници
слоја
.
Одатле
следи
1
1 d
d
d
d
p
p
U
U
x
x
x
ρ
ρ
∂
=
= −
∂
Полазне
једначине
своде
се
на
2
2
d
d
x
x
x
x
y
v
v
v
U
v
v
U
x
y
x
y
ν
∂
∂
∂
+
=
+
∂
∂
∂
0
y
x
v
v
x
y
∂
∂
+
=
∂
∂
.
Иако
је
карактер
последњих
једначина
остао
нелинеаран
,
знатна
,
спроведена
упрошћења
омогућавају
њихово
решавање
.
Једначине
важе
за
опструјавање
равне
плоче
,
а
уколико
је
радијус
кривине
велик
у
односу
на
дебљину
граничног
слоја
,
могу
се
употребити
и
за
закривљене
струјнице
.
Пошто
се
притисак
сматра
познатом
величином
,
две
једначине
су
довољне
за
решавање
непознатих
v
x
и
v
y
.
Међутим
,
уколико
је
притисак
у
граничном
слоју
непознат
(
тела
са
малим
радијусом
кривине
и
оштрим
прелазима
),
због
већег
броја
непознатих
од
расположивих
једначина
,
мора
се
прибећи
поступном
решавању
проблема
који
се
састоји
у
наизменичном
одређивању
струјних
карактеристика
спољашњег
струјног
поља
проширеног
на
гранични
слој
и
поправљањем
њихових
вредности
на
основу
једначина
струјања
у
граничном
слоју
.
Ламинаран
режим
задржава
се
на
целокупној
дужини
граничног
слоја
ако
су
следећи
карактеристични
Рејнолдсови
бројеви
мањи
од
критичних
,
тј
.
0
0
0
1
0
2
1
1
2
2
;
;
;
;
l
l
lk
k
l
l
k
k
U l
U
Re
Re
Re
Re
U
U
Re
Re
Re
Re
δ
δ
δ
ν
ν
δ
δ
ν
ν
=
≤
=
≤
=
≤
=
≤
где
су
:
l
-
дужина
тела
,
δ
l
-
дебљина
граничног
слоја
на
задњем
делу
тела
,
тј
.
највећа
дебљина

114
почетна
нестабилност
се
развија
и
распростире
по
граничном
слоју
.
На
крају
прелазне
области
(
пресек
II)
турбулентни
струјни
режим
развијен
је
по
целој
дебљини
граничног
слоја
,
изузев
ламинарног
подслоја
.
Дебљина
ламинарног
подслоја
је
врло
мала
(
вискозни
или
молекуларни
слој
),
граница
са
турбулентним
слојем
је
прелазног
карактера
.
Промена
параметра
подслоја
није
до
данас
позната
због
немогућности
мерења
струјних
карактеристика
у
самој
близини
граничне
површине
.
Према
неким
ауторима
дебљина
ламинарног
подслоја
се
не
мења
.
Пошто
је
δ
на
известан
начин
условна
величина
,
правилно
је
,
како
је
већ
приказано
,
да
се
карактеристични
Рејнолдсови
бројеви
дефинишу
преко
δ
1
и
δ
2
.
Губитак
стабилности
ламинарног
струјања
дефинише
се
са
:
1
2
1
2
1
2
1
2
;
;
;
;
I
I
I
I
I
I
II
II
II
II
II
II
U
U
Re
Re
U
U
Re
Re
δ
δ
ν
ν
δ
δ
ν
ν
=
=
=
=
Где
су
U
I
и
U
II
брзине
струјања
на
спољашњој
граници
слоја
у
пресецима
I
и
II.
Величине
свих
критичних
Рејнолдсових
бројева
зависе
од
:
параметра
облика
граничног
слоја
,
броја
почетне
турбулентности
и
локалне
релативне
храпавости
.
Одређивање
ових
функција
је
тешко
с
обзиром
да
још
није
јасна
природа
саме
турбуленције
,
па
се
у
практичним
задацима
оне
одређују
експериментално
.
За
турбулентан
гранични
слој
важе
следеће
једначине
осредњеног
кретања
(
величине
су
означене
без
горње
цртице
):
(
)
(
)
2
2
2
d
d
x
x
x
x
y
x
x
y
v
v
v
U
v
v
U
v
v v
x
y
x
y
x
y
ρ
ρ
η
ρ
ρ
∂
∂
∂
∂
∂
′
′ ′
+
=
+
−
−
∂
∂
∂
∂
∂
0
y
x
v
v
x
y
∂
∂
+
=
∂
∂
.
Експерименти
су
показали
да
је
члан
(
)
2
x
v
x
ρ
∂
′
∂
мали
у
поређењу
са
другим
,
па
се
може
занемарити
.
Тада
се
добија
:
115
d
1
d
x
x
x
x
y
x
y
v
v
v
U
v
v
U
v v
x
y
x
y
y
η
ρ
ρ
∂
∂
∂
∂
′ ′
+
=
+
−
∂
∂
∂
∂
0
y
x
v
v
x
y
∂
∂
+
=
∂
∂
са
граничним
условима
за
у
=0,
v
x
=v
y
=0
за
у
=
δ
,
v
x
=U
,
v
y
=0.
Тангенцијални
напон
на
омотачу
тела
τ
0
,
дебљина
граничног
слоја
δ
и
места
која
одвајају
области
,
главне
су
карактеристике
које
треба
да
се
одреде
.
Познавање
профила
брзине
кроз
гранични
слој
омогућава
решавање
овог
задатка
.
Логаритамски
распоред
брзина
добијен
за
турбулентна
струјања
у
цевима
(
које
важи
за
турбулентно
струјања
дуж
бесконачног
зида
),
одговара
стварности
и
за
турбулентан
гранични
слој
коначних
тела
,
тј
.
за
струјања
која
нису
на
великим
растојањима
од
површине
тела
.
Дебљина
граничног
слоја
расте
дуж
површине
тела
око
кога
се
врши
опструјавање
у
смеру
протицања
течности
.
Ово
објашњава
зашто
при
протицању
кроз
цев
логаритамски
профил
постоји
кроз
цео
попречни
пресек
.
Почев
од
места
утицања
у
цев
,
гранични
слој
се
развија
да
би
на
коначном
растојању
испунио
цео
пресек
цеви
.
Због
тога
се
у
довољно
дугој
цеви
не
разликује
гранични
слој
од
остале
зоне
струјања
.
Аналогно
стање
постоји
и
при
ламинарном
протицању
кроз
цев
када
се
улога
вискозности
манифестује
по
целом
пресеку
цеви
,
а
није
ограничена
на
танак
слој
течности
уз
зид
цеви
.
Потреба
за
познавањем
τ
0
,
δ
,
δ
1
,
δ
2
,
чини
недовољним
системе
датих
једначина
за
кретање
у
ламинарном
и
турбулентном
слоју
.
Прантл
је
уочио
погодност
примене
једначине
о
промени
количине
кретања
за
решавање
задатака
у
ламинарном
и
турбулентном
граничном
слоју
,
па
је
тај
метод
због
своје
једноставности
и
директног
увођења
τ
0
опште
прихваћен
.
Једначине
кретања
у
ламинарном
и
турбулентном
слоју
могу
да
се
сведу
на
такав
облик
,
али
се
најчешће
користи
директна
примена
једначине
о
промени
количине
кретања
.

117
Једначина
о
промени
количине
кретања
не
даје
информације
о
профилу
брзине
.
Он
се
претпоставља
на
основу
експерименталних
мерења
,
тако
да
има
исти
облик
за
сваку
вредност
x
.
Распоред
брзине
мора
да
задовољи
и
граничне
услове
:
v
x
=0
за
y
=0
и
v
x
=U
за
y=
δ
.
Слика
8.
Елемент
контролне
запремине
граничног
слоја
Ламинаран
гранични
слој
За
ламинаран
гранични
слој
преко
равне
површине
,
Прантл
је
предложио
општи
распоред
брзина
:
3
1
2
x
v
C y
C y
=
+
.
Коефицијенти
С
1
и
С
2
одређују
се
из
граничних
услова
:
за
,
,
0
x
x
v
y
v
U
y
δ
∂
=
=
=
∂
тако
да
следи
1
2
3
3
;
2
2
U
U
C
C
δ
δ
=
= −
.
Једначина
за
распоред
брзина
унутар
ламинарног
граничног
слоја
добија
облик
3
3
1
,
0
2
2
x
v
y
y
y
U
δ
δ
δ
=
−
≤ ≤
Једначина
за
тангенцијални
напон
,
када
се
у
њу
уврсти
добијени
профил
брзине
постаје
:
(
)
1
2
0
0
0
d
1
d
x
x
x
x
v
v
y
U
v v y
U
x
x
U U
δ
δ
τ
ρ
ρ
δ
∂
∂
=
−
=
−
∂
∂
∫
∫
односно
118
3
3
1
2
2
0
0
3
1
3
1
1
d
0,139
2
2
2
2
y
y
y
y
y
U
U
x
x
δ
δ
τ
ρ
ρ
δ
δ
δ
δ
δ
∂
∂
=
−
+
−
=
∂
∂
∫
На
равној
површини
такође
важи
:
0
0
3
2
x
y
v
U
y
τ
η
η
δ
=
∂
=
=
∂
,
те
изједначавањем
десних
страна
једначина
за
тангенцијални
напон
следи
2
3
0,139
2
U
U
x
δ
η
ρ
δ
∂
=
∂
тј
.
d
d
10, 78
x
U
η
δ δ
ρ
=
односно
после
инеграције
2
10,78
2
x
C
U
δ
ν
=
+
.
Ако
је
δ
=0
за
х
=0,
следи
да
је
и
С
=0,
па
је
4,65
4,65
x
x
Ux
Re
δ
ν
=
=
;
где
је
х
–
растојање
од
водеће
ивице
плоче
Re
x
–
Рејнолдсов
број
базиран
на
растојању
од
водеће
ивице
.
Претходна
једначина
показује
да
се
δ
повећава
са
квадратним
кореном
растојања
од
водеће
ивице
.
Заменом
δ
у
једначину
за
тангенцијални
напон
на
зиду
добија
се
3
0
0,322
U
x
ηρ
τ =
.
τ
0
мења
се
обрнуто
пропорционално
са
квадратним
кореном
растојања
од
водеће
ивице
,
а
директно
пропорционално
са
U
3/2
.
Отпорна
сила
са
једне
стране
плоче
износи
3
0
0
d
0,644
l
x
F
x
U l
τ
ηρ
=
=
∫
Сила
трења
обично
се
изражава
преко
коефицијента
отпора
(
поглавље
Отпори
)
динамичког
притиска
и
оквашене
површине
:

120
1
1
1
7
7
2
0
0
1
d
y
y
y
U
x
δ
τ
ρ
δ
δ
δ
∂
=
−
∂
∫
.
Увођењем
задатог
израза
за
τ
0
,
добија
се
диференцијална
једначина
за
одређивање
δ
1
1
4
4
d
0, 234
d
x
U
ν
δ δ
=
.
Интеграција
се
врши
под
претпоставком
да
је
турбулентан
слој
преко
целе
дужине
плоче
,
па
су
почетни
услови
х
=0,
δ
=0
1
5
4
4
0, 292
x
U
ν
δ
=
,
односно
1
1
5
5
5
0,37
0,37
x
x
x
U
Re
δ
ν
−
=
=
.
Дебљина
граничног
слоја
знатно
брже
расте
у
турбулентном
него
у
ламинарном
граничном
слоју
(
слика
9
),
х
4/5
:
х
1/2
.
Слика
9.
Развој
граничног
слоја
Отпор
трења
глатке
плоче
добија
се
коришћењем
израза
за
τ
0
1
4
2
0
0,0228
U
U
ν
τ
ρ
δ
=
и
заменом
δ
из
израза
5
0,37
x
x
Re
δ
=
па
се
добија
1
5
2
0
0, 0228
U
Ux
ν
τ
ρ
=
121
2
0
5
0
0, 036
d
l
l
U l
F
x
Re
τ
ρ
τ
=
=
∫
5
7
5 10
10
l
Re
⋅
<
<
а
коефицијент
отпора
је
5
0, 072
D
x
c
Re
=
.
Горње
једначине
важе
за
области
нижих
Re
.
За
Re
=400.000,
n
=1/8,
а
за
Re
=4.000.000,
n
=1/10.
Вредности
за
δ
1
и
δ
2
су
1
2
,
1
(1
)(1 2 )
n
n
n
n
n
δ
δ
δ
δ
=
=
+
+
+
.
Коришћењем
логаритамског
профила
где
је
v*
-
привидна
брзина
,
добија
се
*
*
0
*
*
2,5ln
5,5,
,
2,5ln
x
x
v
U
v
v y
y
v
v
v
τ
ν
ρ
δ
−
=
+
=
= −
.
Шлихтинг
је
добио
следећу
вредност
за
коефицијент
отпора
глатке
плоче
(
)
2,58
0, 455
log
D
l
c
Re
=
.
Вредности
δ
1
и
δ
2
за
логаритамски
профил
брзине
су
2
*
*
*
1
2
2,5
,
2,5
v
v
v
U
U
U
δ
δ
δ
δ
=
=
−
.
Прорачун
турбулентног
граничног
слоја
преко
храпаве
површине
спроводи
се
на
сличан
начин
,
користећи
се
мерењем
струјних
карактеристика
на
вештачки
храпавим
површинама
.
На
уструјном
делу
плоче
струјање
је
ламинарно
,
затим
у
турбулентном
граничном
слоју
,
где
је
слој
још
танак
,
а
однос
дебљине
граничног
слоја
и
храпавости
,
δ
/
е
,
још
мали
–
област
је
потпуно
развијене
храпавости
у
којој
је
отпор
пропорционалан
квадрату
брзине
.
За
дугачку
плочу
,
после
овог
дела
јавља
се
нова
област
у
којој
је
δ
/
е
знатно
веће
,
тако
да
плоча
може
да
постане
хидраулички
глатка
,
тј
.
губици
се
не
смањују
са
смањењем
храпавости
.
Вредности
τ
0
t
/
τ
0
l
дате
у
Табели
1
указују
на
овакву
тенденцију
.
Експериментима
је
установљено
да
утицај
храпавости
престаје
при
100
eU
ν
<
.

123
Слика
11.
Распореди
брзина
у
граничном
слоју
Поређење
карактеристика
граничног
слоја
при
ламинарном
и
турбулентном
струјању
дато
је
у
Табели
1
.
Поређење
је
извршено
при
истим
дебљинама
граничног
слоја
δ
и
истим
брзинама
ван
граничног
слоја
U
.
Високе
вредности
τ
0
t
/
τ
0
l
у
турбулентном
слоју
показују
колико
се
исплати
одржавање
ламинарног
слоја
ради
смањења
отпорних
сила
.
Мале
релативне
храпавости
(
δ
/
е
)
односе
се
на
струјање
преко
изразито
храпавих
површина
,
што
узрокује
високе
тангенцијелне
напоне
.
Табела
1.
Поређење
карактеристика
ламинарног
и
турулентног
граничног
слоја
ламинаран
гранични
слој
турбулентан
гранични
слој
логаритамска
промена
параболичан
закон
1/7
степени
закон
δ
/
е
=10
2
δ
/
е
=10
3
δ
/
е
=10
4
δ
1
/
δ
0,33
0,13
0,21
0,14
0,11
δ
2
/
δ
0,13
0,097
0,13
0,1
0,085
τ
0
t
/
τ
0
l
1,0
6,7
22
9,6
5,4
Одвајање
у
граничном
слоју
За
разлику
од
опструјавања
око
плоче
,
при
опструјавању
аеропрофила
или
обртних
тела
,
знатно
се
мења
grad
p
уздуж
турбулентног
граничног
слоја
,
а
зависно
од
њега
распоред
брзине
и
тангенцијални
напон
τ
0
.
На
124
слици
12
приказан
је
облик
граничног
слоја
и
распоред
брзина
у
области
одвајања
.
Види
се
да
уколико
је
∂
p
/
∂
x
<0 (
повољан
grad
p
)
профил
брзине
је
још
увек
пунији
(
у
односу
на
равну
плочу
са
∂
p
/
∂
x
=0);
а
при
∂
p
/
∂
x
>0 (
неповољан
grad
p
)
профил
брзине
је
оштрији
(
стрмији
).
При
позитивним
grad
p
профил
брзине
нагло
се
мења
и
одлепљује
када
је
τ
0
≈
0.
Слика
12.
а
)
Одвајање
у
граничном
слоју
;
б
)
Распоред
брзине
у
граничном
слоју
за
различите
вредности
grad p
До
одлепљивања
флуидне
струје
(
тачка
0)
долази
због
тога
што
је
знатна
количина
флуида
успорена
,
па
се
према
Бернулијевој
једначини
јавља
повећање
притиска
,
које
заједно
са
тангенцијалним
напоном
на
површини
тела
тежи
да
врати
флуид
уназад
.
Тачка
одвајања
налази
се
на
чврстој
површини
,
пошто
је
ту
,
под
дејством
повећаног
притиска
и
тангенцијалног
напона
,
заустављено
кретање
.
Струјање
у
супротном
правцу
које
је
зачетак
вртлога
(
таласа
),
формира
се
ниструјно
од
тачке
одвајања
.
Ефекат
одвајања
и
вртложни
траг
су
најзначајнији
процес
претварања
кинетичке
енергије
у
енергију
притиска
и
губитка
механичке
енергије
такође
.
Аеродинамична
тела
обликују
се
тако
да
се
тачка
одвајања
спусти
што
је
могуће
ниже
.
Уколико
се
одвајање
избегне
,
а
гранични
слој
остане
танак
,
притисак
се
скоро
у
потпуности
врати
флуидној
струји
.
У
том
случају
отпор
је
последица
само
трења
.
У
вртложном
трагу
притисак
се
не
враћа
флуидној
струји
,
тако
да
постоји
отпор
услед
притиска
,
одн
.
отпор
облика
.
Сужавањем
вртложног
трага
смањује
се
отпор
облика
који
је
најчешће
доминантан
(
поглавље
Отпори
).

126
Слика
14.
Одвајање
а
)
ламинарног
и
б
)
турбулентног
граничног
слоја
на
лопти
U=7,5 m/s,
φ
20 mm;
a)
глатка
површина
;
б
)
калота
са
песком
φ
10 mm
Методе
управљања
граничним
слојем
око
аеропрофила
Као
што
је
речено
,
прелаз
ламинарног
режима
у
турбулентан
,
а
такође
и
одвајање
граничног
слоја
,
доводи
до
знатног
пораста
отпора
,
односно
губитка
енергије
.
Због
тога
су
се
развиле
методе
за
очување
ламинарног
граничног
слоја
на
што
је
могуће
већој
дужини
,
а
тиме
и
спречавање
одвајања
флуида
од
чврсте
површине
.
Данас
су
развијене
четири
методе
управљања
граничним
слојем
.
1.
Одсисавање
или
убрзавање
флуида
.
Одсисавањем
дела
флуида
из
зоне
граничног
слоја
са
повећаном
дебљином
(
0
p
x
∂ ∂ >
),
када
је
смањена
кинетичка
енергија
у
граничном
слоју
недовољна
за
компензацију
раста
притиска
и
губитка
услед
трења
,
удаљава
се
успорени
слој
који
своје
место
уступа
флуиду
са
већом
кинетичком
енергијом
.
На
тај
начин
спречава
се
одвајање
и
продужава
ламинаран
режим
.
Одсисавање
може
да
се
врши
и
на
више
места
,
а
ефикасност
методе
је
велика
пошто
може
да
се
добије
оптицање
тела
без
вртложног
трага
како
је
то
приказано
на
слици
15
.
Убризгавањем
флуида
(
слика
16
)
постиже
се
исти
ефекат
као
и
одсисавањем
с
обзиром
да
се
такође
уклања
успорени
слој
флуида
,
а
такође
се
саопштава
додатна
кинетичка
енергија
која
делује
на
убрзање
струјања
у
граничном
слоју
.
Оба
начина
,
због
могућности
спречавања
одвајања
у
граничном
слоју
,
користе
се
у
техници
,
а
нарочито
у
ваздухопловству
.
Ова
метода
127
теоријски
није
довољно
обрађена
,
па
због
тога
још
није
нашла
примену
коју
заслужује
.
Слика
15.
Одсисавање
успореног
флуида
Слика
16.
Убризгавање
додатног
флуида
Слика
17.
Коефицијент
отпора
за
неке
врсте
аеропрофила
2.
Обликовање
чеоног
дела
тела
ради
задржавања
ламинарног
слоја
.
Градијент
притиска
и
почетна
турбулентност
основног
протока
највише
учествују
у
развоју
турбулентног
граничног
слоја
.
Уколико
градијент
притиска
постане
раније
неповољан
,
тачка
одвајања
помера
се
уструјно
,
а
отпор
при
томе
расте
.
Задржавањем
повољног
градијента
притиска
на
што
већој
дужини
тела
условљено
је
обликом
чеоног
дела

129
Слика
18.
Карактеристике
ламинаризованог
слоја
Слика
19.
Асимптотски
гранични
слој
Табела
2.
Поређење
коефицијената
отпора
без
и
са
равномерним
одсисавањем
Re
l
10
7
10
8
10
9
c
Dlam
4,2
⋅
10
-4
1,32
⋅
10
-4
4,2
⋅
10
-5
c
Dtur
3
⋅
10
-3
2,1
⋅
10
-3
1,6
⋅
10
-2
c
Da
2,36
⋅
10
-4
2,36
⋅
10
-4
2,36
⋅
10
-4
4.
Ламинаризација
граничног
слоја
са
прекидним
одсисавањем
.
Одсисавањем
кроз
специјалне
прорезе
или
решетке
могу
се
добити
знатна
снижења
отпора
.
Спољашња
граница
слоја
тада
има
степенаст
облик
како
је
то
приказано
на
слици
20
.
Однос
између
коефицијента
130
отпора
са
и
без
одсисавања
одређује
се
експериментално
,
а
зависи
од
0
0
sr
v
U
,
где
је
0
sr
v
средња
брзина
одсисавања
,
према
дијаграму
на
слици
21
.
Слика
20.
Степенаст
гранични
слој
Слика
21.
Зависност
коефицијента
отпора
са
и
без
одсисавања
У
данашње
време
постоје
и
друге
методе
снижавања
хидродинамичког
трења
,
као
што
су
примена
пригушних
покривки
,
увођење
гаса
у
гранични
слој
,
формирање
гасних
прстенова
,
уношење
полимера
,
вибрирање
омотача
тела
и
др
.
Реално
је
очекивати
да
ће
се
у
будућности
појавити
још
ефикасније
методе
за
снижавање
хидродинамичког
трења
,
чиме
ће
хидро
механика
значајно
допринети
развоју
и
убрзању
техничког
прогреса
.

132
2
2
2
,
, , ,...
D
e
F
v f
Re Ma
ρ
φ
=
,
где
се
функција
f
2
може
заменити
коефицијентом
отпора
c
D
различитим
за
сваку
врсту
опструјавања
тела
и
зависним
од
многобројних
утицаја
.
Увођењем
динамичког
притиска
2
1
2
v
ρ
уместо
2
v
ρ
и
замена
2
са
карактеристичном
површином
тела
А
добија
се
формула
по
којој
се
најчешће
одређује
сила
отпора
2
1
.
2
D
D
F
v c A
ρ
=
Отпор
на
тело
уроњено
у
струју
флуида
састоји
се
од
отпора
трења
и
отпора
облика
.
Горња
формула
,
добијена
комбинованом
применом
димензијске
анализе
и
експериментима
измерених
података
,
важи
за
обе
врсте
,
с
тим
што
је
потребно
водити
рачуна
о
дефиницији
карактеристичне
површине
.
Отпор
трења
Својства
отпора
трења
разматраће
се
за
равну
,
танку
,
глатку
или
храпаву
плочу
у
струји
флуида
.
Овај
отпор
зависи
од
формираног
граничног
слоја
,
како
је
већ
расправљано
у
претходној
глави
.
Коефицијент
отпора
трења
дефинисан
је
за
оквашену
површину
.
A
b
= ⋅
Дужина
плоче
је
у
правцу
кретања
плоче
(
или
струје
флуида
,
ако
је
плоча
непокретна
)
а
R
е
број
,
од
кога
зависи
врста
граничног
слоја
,
је
v
ν
.
За
ламинаран
гранични
слој
са
распоредом
брзине
и
дебљином
2
3
30
3
1
4,65
2
;
;
или
;
;
2
2
x
x
v
v
y
y
y
y
U
x
Re
U
x
Re
δ
δ
δ
δ
δ
δ
=
−
=
=
−
=
коефицијент
отпора
је
1,33
D
c
Re
=
за
4
5
10
6 10
Re
<
< ⋅
(
)
.
Re
Re
=
За
вредности
R
е
<10
3
горња
формула
не
даје
задовољавајуће
вредности
.
Треба
је
заменити
са
133
1,33
4,12
D
c
Re
Re
=
+
за
R
е
<6
⋅
10
3
.
За
турбулентан
гранични
слој
преко
глатке
површине
са
распоредом
брзине
и
дебљином
1
7
1 5
0,38
;
,
v
y
U
x
Re
δ
δ
=
=
коефицијент
отпора
одређен
је
са
1 5
0,074
D
c
Re
=
за
10
5
<
R
е
<2
⋅
10
7
За
турбулентан
гранични
слој
преко
глатке
површине
,
са
логаритамском
расподелом
брзине
,
коефицијент
отпора
одређен
је
са
(
)
2,58
0, 455
log
D
c
Re
=
за
10
6
<
R
е
<10
9
.
За
турбулентан
гранични
слој
преко
храпаве
површине
коефицијент
отпора
зависи
само
од
релативне
храпавости
e
(
не
зависи
од
R
е
).
На
слици
1
приказане
су
наведене
промене
коефицијента
отпора
трења
у
зависности
од
R
е
броја
.
За
надзвучну
струју
уздуж
глатке
површине
коефицијент
отпора
се
смањује
са
повећањем
М
броја
,
како
је
показано
на
слици
.
За
ламинарну
струју
је
смањење
око
15 %
при
М
=4.
За
турбулентну
струју
(
R
е
=10
4
)
смањење
коефицијента
отпора
је
око
50 %
при
М
=4.
Слика
1.
Коефицијенти
отпора
услед
трења
равне
плоче

135
2
2
9
1
sin
,
1
4
2
o
o
p
p
v
θ
ρ
−
= −
како
је
то
приказано
на
слици
3
а
.
Сл
. 3.
Распоред
притиска
и
брзине
за
струјање
око
лопте
136
При
ламинарном
опструјавању
лопте
нормалан
вискозан
напон
и
пројекција
силе
притиска
правцем
долазеће
струје
(
x
-
осе
)
износе
0
|
3
cos ;
2
rr r=R
v
p
R
η
θ
= −
2
;
xp
o
F
Rv
πη
=
док
су
тангенцијални
напони
и
пројекцијалне
силе
на
x
-
осу
0
|
3
sin ;
2
r r= R
v
R
θ
τ
η
θ
=
4
.
x
o
F
Rv
τ
πη
=
Укупна
отпорна
сила
3
D
o
F
Dv
πη
=
,
односно
силе
вискозног
трења
троше
енергију
притиска
тим
више
што
је
гранични
слој
дебљи
,
па
се
због
тога
формира
несиметричан
распоред
притиска
дуж
осе
струјања
како
је
то
приказано
на
слици
3
б
.
Када
су
Рејнолдсови
бројеви
мали
,
нема
инерцијских
ефеката
па
се
поремећаји
услед
оптицања
тела
преносе
далеко
од
граничних
површина
тела
.
Тако
нпр
.
при
опструјавању
лопте
са
R
е
<
0,5,
тек
на
удаљењу
четири
пречника
од
лопте
(
вертикална
оса
),
брзина
је
92 %
од
брзине
неузнемирене
струје
.
Одвајање
ламинарног
граничног
слоја
одговара
слици
3
ц
којој
одговарају
и
одвајања
турбулентног
слоја
све
до
вредности
R
е
=4,3
÷
4,5
⋅
10
5
када
долази
до
спонтане
турбулизације
целокупног
граничног
слоја
.
Ова
тачка
назива
се
кризом
отпора
и
њој
одговарају
слика
3
д
и
тачка
К
на
слици
5
.
Интензивна
турбуленција
знатно
повећава
способност
задржавања
граничног
слоја
тако
да
се
одвајање
дешава
на
задњем
делу
лопте
,
а
вртложни
траг
који
је
мерило
величине
отпора
је
и
три
пута
ужи
од
оног
на
слици
3
ц
.
Криза
отпора
је
појава
карактеристична
и
за
друге
профиле
.
У
тачки
К
не
само
да
је
c
D
опало
на
отприлике
1/3
вредности
претходне
аутомоделне
области
када
је
c
D
=const. (
не
зависи
од
R
е
),
него
је
и
вредност
укупног
отпора
мања
без
обзира
што
је
отпор
сразмеран
v
2
.
Опструјавање
тела
без
одвајања
прати
гранични
слој
какав
је
показан
на
слици
4
.
На
задњем
делу
тела
нема
карактеристичног
одвајања
и
вртложног
трага
.
Формира
се
“
задебљани
”
гранични
слој
,
односно
вртложни
траг
без
изразитих
вртлога
који
захтевају
већи
утрошак
енергије
.
Вртложни
траг
типа
задебљаног
граничног
слоја
дозвољава
да
се
притисак
спољашњег
струјања
пренесе
и
у
траг
и
на
тај
начин
овакву
струјну
слику
приближава
опструјавању
идеалног
флуида
.
Ово
је

138
узгона
није
изражена
(
слика
3
б
).
Несиметричан
распоред
притиска
у
односу
на
осу
симетрије
у
правцу
струјања
обезбеђује
силу
узгона
.
Облик
тела
које
флуид
опструјава
или
које
се
креће
кроз
непокретан
или
покретан
флуид
(
све
анализе
спроводе
се
преко
релативних
брзина
)
одређује
распоред
притиска
око
њега
па
је
одомаћен
израз
отпор
облика
идентичан
са
отпором
(
несиметричног
распореда
)
притиска
.
Отпори
при
малим
и
умереним
R
е
Отпор
при
малим
R
е
.
За
R
е
<
1
опструјавање
тела
је
ламинарно
а
облик
тела
је
другоразредног
значаја
у
односу
на
димензије
тела
,
вискозност
флуида
и
брзину
струје
.
Озенова
корекција
ламинарног
опструјавања
лопте
са
уношењем
инерцијских
утицаја
даје
следећи
израз
за
коефицијент
отпора
24
3
1
16
D
c
Re
Re
=
+
1.
Re
≤
Сила
отпора
одређује
се
преко
2
1
2
D
D
F
v c A
ρ
=
2
,
4
vd
d
Re
A
π
ν
=
=
где
је
d
пречник
лопте
.
За
диск
нормалан
на
ламинарну
струју
коефицијент
отпора
је
20, 4
,
D
c
Re
=
vd
Re
ν
=
за
R
е
<
0,5,
а
за
диск
паралелан
струји
је
13,6
,
D
c
Re
=
vd
Re
ν
=
за
R
е
<
0,1,
а
сила
отпора
одређује
се
према
обе
површине
диска
,
тј
.
2
2
1
2
.
2
4
D
D
d
F
v c
π
ρ
=
Коефицијент
отпора
за
дводимензијско
ламинарно
струјање
око
кружног
цилиндра
је
(
)
8
,
2, 2
ln
D
c
Re
Re
π
=
−
vb
Re
ν
=
за
R
е
<0,1,
139
где
је
b
мања
димензија
плоче
.
За
дводимензијску
ламинарну
струју
преко
равне
плоче
која
је
паралелна
струји
,
за
врло
мале
R
е
бројеве
коефицијент
отпора
одређен
је
са
4,12
,
D
c
Re
=
v
Re
ν
=
за
R
е
<0,01,
а
сила
отпора
за
обе
површине
плоче
одређена
са
2
1
2 .
2
D
D
F
v c
b
ρ
=
Отпор
при
умереним
R
е
.
У
дијаграмима
приказаним
на
сликама
5
и
6
и
табели
1
дате
су
промене
коефицијента
отпора
са
R
е
до
R
е
≈
10
6
,
када
за
већину
тела
долази
до
одвајања
флуидне
струје
од
граничних
површина
,
те
отпорна
сила
постаје
функција
и
других
утицаја
(
нестационарности
,
М
).
Слика
5.
Коефицијент
отпора
за
лопту
и
друга
обртна
тела
У
аутомоделној
области
2
6
5 10
2
5 10 ,
Re
⋅
<
< ÷ ⋅
када
c
D
не
зависи
од
R
е
дешава
се
одвајање
граничног
слоја
за
тела
приказана
на
дијаграму
.
Тада
је
отпор
услед
вртложења
95 %
укупног
отпора
.
Већ
споменута
криза
отпора
тачку
одвајања
граничног
слоја
помера
знатно
низструјно
(
тачке
К
).

141
кризе
отпора
за
цилиндар
коначне
дужине
.
Формирање
вртлога
је
наизменично
(
слика
8
).
Види
се
да
одвајање
граничног
слоја
представља
одвајање
целокупног
протока
(
граничног
слоја
и
спољашњег
струјања
)
од
површине
тела
.
Место
упражњене
количине
флуида
попуњава
се
повратним
струјањем
и
концентрисаним
вртлозима
.
При
одвајању
граничног
слоја
параметри
струјања
се
толико
измене
да
се
слободно
може
рећи
да
гранични
слој
престаје
да
постоји
,
па
се
овакав
вртложни
траг
не
може
поредити
са
вртложним
трагом
без
одвајања
.
Вртлози
се
крећу
спорије
од
неузнемирене
струје
(
v
),
а
одвајају
се
од
тела
брзином
v
v
≈
0,2
v
.
Слика
7.
Распоред
притиска
при
опструјавању
цилиндра
Слика
8.
Наизменично
формирање
вртлога
при
опструјавању
цилиндра
Карманов
вртложни
траг
који
се
појављује
ниструјно
од
цилиндра
,
резултат
је
наизменичног
формирања
вртлога
и
њиховог
одвајања
од
цилиндра
,
због
чега
и
наизменично
дејствују
силе
на
цилиндар
.
Однос
уздужног
и
попречног
растојања
вртлога
за
кружни
цилиндар
износи
а
/b
=3,55,
а
за
равну
плочу
око
3,3.
142
Слика
9.
Карманов
низ
вртлога
Фреквенција
одвајања
вртлога
може
се
одредити
брзином
којом
се
вртлози
удаљују
од
тела
v
v
и
растојањем
вртлога
у
трагу
а
,
тј
.
v
v
f
a
=
.
Фреквенција
је
функција
кинематске
вискозности
флуида
ν
,
брзине
флуидне
струје
v
и
пречника
цилиндра
d
,
што
се
може
комбиновати
у
зависност
два
бездимензијска
параметра
:
Струхаловог
и
Рејнолдсовог
броја
.
(
)
fd
Str
Re
v
φ
=
=
Ова
зависност
приказана
је
на
дијаграму
на
слици
6
за
одређивање
коефицијента
отпора
цилиндра
.
Аналитичко
одређивање
коефицијента
вртложног
отпора
при
одвајању
граничног
слоја
још
није
добијено
сем
приказаног
случаја
опструјавања
цилиндра
.
Вртложни
низ
који
је
у
наизменичном
шаховском
распореду
,
назива
се
Карманов
двоструки
низ
вртлога
,
а
према
Карману
клоефицијент
отпора
вртложног
спрега
дат
је
са
2
1,59
0, 63
v
v
Dv
Re
Re
d
c
Re
a
Re
=
−
где
су
:
;
v
v
av
vd
Re
Re
ν
ν
=
=
d
–
пречник
цилиндра
.
Откидање
вртлога
прати
појава
импулсне
силе
према
закону
о
промени
количине
кретања
.
Уколико
се
фреквенција
одвајања
вртлога
поклопи
са
фреквенцијом
сопствених
осцилација
тела
које
флуид
опструјава
,
долази
до
резонанције
и
до
претераног
пораста
амплитуде
осцилације
тела
,
што
може
да
проузрокује
лом
тела
(
рушење
димњака
,
откидање

144
Отпор
у
компресибилној
струји
везан
је
за
велике
брзине
струјања
,
па
тиме
и
велике
Re
бројеве
,
те
се
отпор
при
великим
Re
посебно
не
анализира
.
Отпор
аеропрофила
и
сила
узгона
На
кружни
цилиндар
дужине
l
који
се
обрће
,
хоризонтална
,
једнолика
флуидна
струја
дејствује
вертикалном
силом
,
узгоном
,
која
се
одређује
са
0
L
F
lv
ρ
=
Γ
где
је
d
v l
Γ =
∫
.
Мада
се
обртни
цилиндар
данас
не
користи
за
стварање
узгонске
силе
,
његова
анализа
помаже
у
схватању
основног
принципа
појаве
узгонске
силе
код
аеропрофила
.
Аеропрофили
Према
дефиницији
аеропрофил
је
обликовано
тело
које
изазива
аеродинамичку
силу
нормално
на
његово
кретање
кроз
флуид
са
малим
отпором
.
Ова
сила
је
узгонска
F
L
и
заједно
са
отпорном
F
D
одређује
укупну
силу
F
којом
флуид
дејствује
на
аеропрофил
.
Карактеристична
дужина
аеропрофила
означена
је
са
l
,
његов
распон
са
b
,
а
носећа
карактеристична
површина
профила
је
A=bl
.
Силе
отпора
и
узгона
дефинисане
су
са
:
2
1
2
D
D
F
v c A
ρ
=
и
2
1
2
L
L
F
v c A
ρ
=
где
су
c
L
и
c
D
коефицијенти
узгона
и
отпора
и
једнаки
су
за
све
сличне
моделе
аеропрофила
при
истим
нападним
угловима
.
На
слици
11
приказане
су
карактеристичне
криве
промене
коефицијента
отпора
аеропрофила
за
различите
нападне
углове
α
.
145
Слика
10.
Силе
на
аеропрофилу
Слика
11.
Коефицијенти
отпора
аеропрофила
Постанак
циркулације
При
врло
малим
брзинама
када
аеропрофил
почиње
да
се
креће
из
стања
мировања
,
слика
опструјавања
у
близини
задње
ивице
аеропрофила
изгледаће
као
на
слици
12
а
,
са
зауставном
тачком
0
на
леђној
страни
,
на
извесној
малој
удаљености
са
задње
ивице
.
Уколико
брзина
аеропрофила
расте
,
струјнице
са
грудне
стране
нису
више
у
стању
да
обиђу
око
задње
ивице
због
великих
сила
вискозности
,
које
почињу
да
делују
услед
великог
градијента
брзине
,
те
се
струјање
одвија
од
задње
ивице
,
а
у
некој
тачки
,
између
ове
и
малопређашње
зауставне
,
ствара
се
вртлог
,
као
што
је
показано
на
слици
12
б
који
се
назива
почетни
вртлог
јер
се
формира
покретањем
аеропрофила
.
Када
овај
вртлог
порасте
до
извесне
величине
,
он
се
одваја
од
аеропрофила
и
одлази
низ
струјно
поље
.
Пошто
циркулација
око
било
какве
контуре
мора
бити
једнака
нули
,
јер
је
таква
била
и
пре
кретања
,
око
аеропрофила
мора
постојати
циркулација
Г
која
је
по
величини
једнака
,

147
опише
појавом
вертикалне
компоненте
брзине
w
која
ефективну
брзину
усмерава
наниже
,
тј
.
привидно
се
смањује
нападни
угао
на
α
i
као
што
је
приказано
на
слици
14
.
Из
ове
анализе
следи
D
Do
Di
F
F
F
=
+
где
су
:
F
Do
–
сила
отпора
на
бесконачном
профилу
услед
брзине
v
F
Di
–
индуковани
отпор
.
Такође
је
:
и
i
o
D
Do
Di
c
c
c
α α
α
′
=
+
=
+
.
Пошто
је
обично
угао
α
i
изазван
индукованом
силом
мали
,
следи
,
,
L
L
Di
L
i
F
F
v
v
F
F
α
′
′
=
=
=
;
и
на
основу
елиптичног
распореда
силе
узгона
може
се
наћи
2
и
L
L
Di
i
c
c
c
b
b
l
l
α
π
π
=
=
.
Горња
једначина
показује
да
уколико
је
распон
крила
већи
и
карактеристична
дужина
крила
мања
,
коефицијент
индукованог
отпора
се
смањује
.
Слика
14.
Индуковани
отпор
Потребне
променљиве
c
Di
,
c
D
и
c
L
у
зависности
од
нападног
угла
α
,
обично
се
дају
у
дијаграму
званом
полара
(
слика
15
).
Затупљени
облик
аеропрофила
осигурава
велики
узгон
при
малим
брзинама
и
релативно
малу
отпорну
силу
.
Међутим
,
при
великим
брзинама
отпор
знатно
расте
,
те
је
у
данашње
време
изражена
тежња
ка
смањењу
коефицијента
отпора
,
макар
то
за
собом
повлачило
и
смањење
коефицијента
узгона
.
148
Како
је
већ
било
приказано
раније
,
градијент
брзине
у
близини
чврсте
површине
много
је
мањи
у
ламинарној
него
у
турбулентној
струји
.
Ламинаран
гранични
слој
обезбеђује
мањи
отпор
трења
,
а
условљен
је
постепеном
променом
притиска
дуж
површине
.
На
слици
16
приказани
су
старији
и
новији
тип
аеропрофила
где
је
уочљиво
да
се
промена
из
ламинарног
у
турбулентан
гранични
слој
за
старији
профил
дешава
на
предњем
делу
профила
,
а
за
новији
тип
на
задњем
делу
профила
.
Отпор
доброг
аеропрофила
може
да
буде
знатно
испод
отпора
равне
плоче
чија
је
дужина
једнака
тетиви
аеропрофила
при
истом
Рејнолдсовом
броју
,
што
се
објашњава
дужим
одржањем
ламинарног
граничног
слоја
на
аеропрофилу
од
оног
на
равној
плочи
.
Слика
15.
Полара

150
На
фотографији
на
слици
18
приказана
је
појава
вртлога
при
опструјавању
тела
различитих
облика
.
Слика
18.
Опструјавање
тела
са
формирањем
вртлога
.
151
Теорија
сличности
и
димензијска
анализа
Теорија
сличности
Механика
флуида
у
другој
половини
XX
века
доживљава
изузетно
интензиван
развој
теоријских
метода
и
експерименталних
истраживања
.
Рачунарска
техника
је
знатно
проширила
области
теоријског
решавања
али
је
примена
механике
флуида
продрла
у
многе
,
за
њу
нове
,
области
,
па
број
нерешених
проблема
расте
.
Захваљујући
много
прецизнијим
мерним
инструментима
и
ефикаснијој
обради
мерних
података
,
експериментална
истраживања
су
задржала
,
ако
не
и
повећала
,
свој
значај
.
Моделска
испитивања
су
најефикаснији
метод
експерименталног
истраживања
јер
са
несразмерно
мањим
улагањима
добијају
се
параметри
прототипскбих
конструкција
.
Резултати
моделских
испитивања
могу
се
искористити
само
уколико
постоји
сличност
струјања
око
модела
и
прототипа
.
Елементе
сличности
дефинише
теорија
сличности
која
одређује
групу
услова
који
омогућавају
уопштавање
експериментално
добијених
резултата
.
Основни
став
ове
теорије
је
да
су
физички
појаве
сличне
,
ако
се
величине
једне
појаве
могу
добити
из
величина
друге
,
узетих
у
одговарајућим
просторно
-
временским
тачкама
,
једноставним
умножавањем
.
Тако
се
две
физички
сличне
појаве
разматрају
само
у
мерилу
и
интензитетима
функција
појаве
.
За
потпуну
сличност
треба
постићи
геометријску
,
кинематску
и
динамичку
сличност
струјања
флуида
.
Геометријски
слична
струјања
захтевају
модел
геометријски
сличан
прототипу
,
односно
исту
размеру
свих
одговарајућих
дужина
на
моделу
и
прототипу
,
што
обезбеђује
и
једнакост
углова
.
Геометријска
сличност
је
сличност
облика
.
Кинематска
сличност
захтева
сличност
струјања
,
тј
.
истоветну
струјну
слику
(
струјнице
,
углови
оптицања
,
вртложења
)
при
чему
се
кинематске
величине
у
истим
тачкама
модела
и
прототипа
разликују
за
константан
број
пута
.
Кинематска
сличност
подразумева
геометријску
а
означава
сличност
кретања
.

153
p
m
c
ρ
ρ
ρ
=
-
мерило
густине
.
Ова
мерила
не
зависе
од
координата
и
времена
.
Ако
се
сваки
члан
Навије
-
Стоксових
једначина
за
опструјавање
прототипа
изрази
преко
одговарајућих
мерила
и
моделских
величина
добија
се
2
2
v
v
p
v
xm
xm
xm
xm
m
m
m
xm
ym
zm
F
m
m
m
m
xm
t
m
l
m
m
m
l
m
l
c c
c c
c
c c c
v
v
v
v
p
v
v
v
c c
X
v
c
t
c
x
y
z
c
x
c
ρ
ρ
ρ ν
ρ
ρ
ρ
ρ
ρ ν
∂
∂
∂
∂
∂
+
+
+
=
−
+
∆
∂
∂
∂
∂
∂
(3)
Поређење
са
једн
. (2)
показује
да
се
сваки
члан
једн
. (3)
разликује
од
одговарајућег
члана
једн
. (2)
за
групу
бездимензијских
мерила
.
Према
другом
услову
динамичке
сличности
однос
истоимених
сила
прототипског
и
моделског
оптицања
мора
бити
исти
,
што
испуњава
услов
2
2
.
v
v
p
v
F
t
l
l
l
c c
c c
c
c c c
c c
const
c
c
c
c
ρ
ρ
ρ ν
ρ
=
=
=
=
=
(4)
где
поједини
чланови
редом
представљају
односе
:
локалних
сила
инерције
,
конвективних
сила
инерције
,
запреминских
сила
,
сила
притиска
и
вискозних
сила
.
Дељењем
релације
(4)
са
другим
чланом
добија
се
2
2
1
p
l
F
l
v t
v
v
v l
c
c
c c
c
c c
c
c c
c c
ν
ρ
= =
=
=
,
одакле
је
2
2
1;
1;
1;
1.
p
v t
v
v l
l
F
l
v
c
c c
c
c c
c
c c
c c
c
ρ
ν
=
=
=
=
=
(5)
Ако
се
у
једначину
(5)
мерила
замене
одговарајућим
величинама
,
добија
се
p p
m m
p
m
p
m
v t
v t
vt
H
H
H
l
l
l
=
=
=
=
=
-
број
хомохроности
,
2
2
2
p
m
p
m
p p
m m
v
v
v
Fr
Fr
Fr
F l
F l
Fl
=
=
=
=
=
-
Фрудов
број
,
154
2
2
2
p
m
p
m
p
p
m m
p
p
p
Eu
Eu
Eu
v
v
v
ρ
ρ
ρ
∂
∂
∂
=
=
=
=
=
-
Ојлеров
број
,
p p
m m
p
m
p
m
v l
v l
vl
Re
Re
Re
ν
ν
ν
=
=
=
=
=
-
Рејнолдсов
број
.
Одавде
следи
теорема
:
два
геометријски
слична
струјања
вискозног
нестишљивог
флуида
су
динамички
сличне
ако
су
Фрудов
,
Ојлеров
,
Рејнолдсов
и
број
хомохроности
једнаки
.
Критеријум
хомохроности
представља
однос
локалних
и
конвективних
инерцијских
сила
.
Фрудов
број
представља
однос
конвективних
инерцијских
сила
и
спољашњих
запреминских
сила
.
У
случају
када
на
струјање
утиче
само
сила
Земљине
теже
,
Фрудов
број
представља
се
са
2
Fr
v
gl
=
или
Fr
v
gl
=
.
Ојлеров
број
је
однос
конвективних
инерцијских
сила
и
сила
притиска
,
а
Рејнолдсов
број
однос
конвективних
инерцијских
сила
и
вискозних
сила
.
Наведени
критеријуми
добијени
су
поређењем
сила
које
дејствују
у
флуиду
и
конвективне
инерцијске
силе
(
у
одсуству
локалног
убрзања
конвективне
инерцијске
силе
називају
се
инерцијским
силама
).
Уколико
се
силе
које
одржавају
равнотежу
при
кретању
флуида
пореде
са
неком
другом
силом
,
добијају
се
други
критеријуми
сличности
често
сложени
од
већ
наведених
.
Табела
1.
Чешће
употребљаване
бездимензијске
карактеристике
Озн
.
Израз
Назив
Примедба
Употреба
Re
vl
ν
Рејнолдсов
број
(Reynolds)
представља
однос
инерцијских
и
вискозних
сила
критеријум
сличности
при
струјању
реалног
(
вискозног
)
флуида
Fr
2
v
gl
Фрудов
број
(Froud)
представља
однос
инерцијских
и
гравитационих
(
спољних
,
запреминских
)
сила
критеријум
сличности
при
кретању
флуида
под
дејством
гравитацијских
сила
М
v
c
Махов
број
(Mach)
представља
однос
инерцијских
и
еластичних
сила
критеријум
сличности
при
струјању
стишљивог
флуида
W
е
2
v l
ρ
γ
Веберов
број
(Weber)
представља
однос
инерцијских
и
сила
и
сила
површинског
напона
критеријум
сличности
када
на
кретање
флуида
утичу
површински
напони

156
Критеријуми
сличности
(
H
,
Fr
,
Е
u
,
R
е
)
везани
су
за
силе
које
се
јављају
у
Навије
-
Стоксовим
једначинама
па
за
проширењем
броја
критеријума
сличности
не
постоји
потреба
,
тим
пре
што
је
немогућ
e
постићи
апсолутну
сличност
два
струјања
,
односно
једнакост
свих
критеријума
сличности
.
Најчешће
се
моделска
испитивања
могу
спровести
уз
одржавање
једнакости
само
једног
критеријума
сличности
,
о
чему
сведоче
бројни
примери
,
па
је
моделска
испитивања
потребно
допунити
корекцијама
које
обухватају
утицаје
и
осталих
сила
.
Једнакост
бездимензијског
критеријума
везује
се
за
утицај
доминантних
величина
,
па
стога
важе
следећа
правила
:
o
при
кретању
течности
код
којих
су
изражене
вискозне
силе
критеријум
сличности
је
R
е
o
кретање
течности
у
пољу
Земљине
теже
треба
моделирати
једнакошћу
R
е
и
Fr
;
међутим
,
како
се
оба
критеријума
у
обично
расположивим
условима
не
могу
задовољити
,
усваја
се
Fr
критеријум
сличности
,
а
утицај
вискозности
обрачунава
посебно
o
при
истакнутим
нестационарним
струјањима
користи
се
један
од
следећих
критеријума
H
,
P
е
,
St
.
Ојлеров
број
ређе
се
употребљава
као
критеријум
сличности
код
течности
зато
што
се
при
кретању
силе
притиска
ретко
експлицитно
истичу
.
Међутим
,
при
струјању
гасова
силе
притиска
замењују
се
врло
активним
еластичним
силама
па
је
однос
инерцијских
и
еластичних
сила
основни
критеријум
сличности
у
динамици
гасова
.
Овај
однос
познат
је
као
Ма
x
ов
број
M=v/c
.
За
сличност
турбулентног
опструјавања
геометријски
сличног
модела
и
прототипа
поред
једнакости
Fr
,
Fo
и
R
е
броја
који
су
дефинисани
у
односу
на
осредњену
брзину
v
потребна
је
и
једнакост
броја
турбулентности
N
ij
дефинисаног
са
2
2
jp ip
jm im
ji
p
m
v v
v v
N
v
v
′ ′
′ ′
=
=
где
су
,
j
i
v v
′ ′
брзине
одступања
у
два
нормална
правца
.
Сврха
моделирања
је
:
o
студија
општег
спољашњег
изгледа
и
геометрије
и
релативна
пропорционалност
различитих
фигура
;
157
o
студија
струјног
т
o
ка
,
величине
и
правца
брзине
и
убрзања
,
преко
,
кроз
и
око
објекта
;
o
проучавање
распореда
притиска
и
резултујућих
сила
на
телу
или
његовим
деловима
;
o
одређивање
струјног
протока
и
калибрација
различитих
струјних
пролаза
;
o
одређивање
хидрауличког
степена
корисности
хидрауличких
машина
;
o
одређивање
енергијских
губитака
услед
трења
и
распореда
притиска
;
o
процена
квалитета
рада
прототипа
;
o
утврђивање
пута
побољшања
конструкције
или
рада
прототипа
;
o
одређивање
начина
за
смањење
цене
прототипа
.

159
Користећи
принцип
димензијске
хомогености
може
се
нпр
.
показати
да
једначина
2
1
2
p
v
gz
C
ρ
ρ
+
+
=
,
представља
могућу
везу
између
притиска
p
,
брзине
v
и
висине
изнад
неке
референтне
равни
z
за
струјање
делића
идеалног
флуида
густине
ρ
дуж
струјнице
,
и
може
да
се
одреди
димензија
константе
C
.
Први
члан
2
1
2
2
сила
MLT
p
ML T
површина
L
−
−
−
=
=
=
.
Други
члан
2
2
2
1
2
3
2
1
2
маса
дужина
M L
v
ML T
запремина
време
L T
ρ
−
−
=
=
=
.
Трећи
члан
1
2
3
2
маса
дужина
M L
gz
дужина
L
ML T
запремина време
L T
ρ
−
−
=
=
=
.
Сви
чланови
једначине
имају
исту
димензију
па
и
константа
C
мора
да
има
димензију
М
L
-1
Т
-2
.
Основу
димензијске
анализе
,
поред
принципа
димензијске
хомогености
,
чини
и
π
теорема
(Vaschy-Buckingham 1915).
Нека
се
физичка
појава
изражава
помоћу
n
величина
(
x
1
,
x
2
,...
x
n
).
Зависно
од
области
физике
где
се
та
појава
јавља
,
постојаће
p
величина
које
се
изражавају
преко
основних
јединица
(
нпр
.
у
динамици
p
=3),
док
ће
се
преосталих
m=n-p
величина
изражавати
помоћу
изведених
јединица
.
π
теорема
гласи
:
физички
закон
који
се
математички
може
представити
(
)
1
2
,
,...
0
n
f x x
x
=
(6)
може
да
се
изрази
прекоједначине
(
)
1
2
,
,...
0
m
F
π π
π
=
(7)
где
су
1
2
,
,...
m
π π
π
бездимензијски
мономи
који
су
функције
1
2
,
,...
n
x x
x
.
160
π
теорема
има
сврху
o
да
се
здруживањем
величина
са
димензијом
добију
бездимензијски
изрази
чија
бројна
вредност
не
зависи
од
изабраног
система
јединица
;
o
да
се
смањи
број
променљивих
за
број
основних
јединица
;
o
да
се
олакша
теоријско
и
експериментално
проучавање
постављених
проблема
.
Пример
.
Показати
да
губитак
притиска
по
јединици
дужине
цеви
,
при
протицању
флуида
кроз
цеви
кружног
пресека
може
да
се
изрази
једначином
2
,
p
v
vd e
l
d
d
ρ
φ
ν
∆
=
;
где
су
d
пречник
цеви
,
l
дужина
цеви
,
е
апсолутна
храпавост
цеви
,
ν
кинематска
вискозност
флуида
и
v
средња
брзина
флуида
.
Пад
притиска
по
јединици
дужине
цеви
може
да
се
изрази
једначином
a
b
c
f
g
i
p
B
d v e
l
ρ
ν
∆
=
∑
.
Димензије
појединих
величина
су
1
2
p
ML T
−
−
∆ =
,
3
ML
ρ
−
=
,
l
L
=
,
1
v
LT
−
=
,
d
L
=
,
e
L
=
,
2
1
L T
ν
−
=
.
Димензијска
једначина
гласи
1
2
3
1
2
1
(
)
(
)
(
)
a
b
c
f
g
ML T
ML
L LT
L L T
L
−
−
−
−
−
=
или
2
2
3
2
a
a
b
c
c
f
g
g
ML T
M L
L L T L L T
−
−
−
−
−
=
.
Изједначавањем
степена
уз
основе
М
, L
и
Т
добијају
се
једначине
М
: 1 =
а
,
L: -2 = -3
а
+
b
+
c
+
f
+2
g
,
Т
: -2 = -
c
-
g
.
Долази
се
до
три
једначине
са
пет
непознатих
величина
(
строго
посматрано
-
две
једначине
са
четири
непознате
).
Поред
величине
а
која
је
једнозначно
одређена
,
треба
одабрати
још
две
величине
које
ће
се
изразити
преко
преосталих
двеју
.
Приликом
овог
избора
мора
се
руководити
искуством
или
унапред
замишљеним
обликом
крајњег
израза
.
У
овом
случају
због
експлицитно
изражених
величина
ρ
,
v
и
d

162
Литература
1.
Букуров
,
Ж
.,
Механика
флуида
,
Факултет
техничких
наука
,
Нови
Сад
, 1984.
2.
Цвијановић
,
П
.,
Предавања
из
механике
флуида
са
карактеристичним
примерима
,
Стилос
,
Нови
Сад
, 1997.
163
Садржај
Увод
у
механику
флуида
..................................................................................................
Error! Bookmark not defined.
Историја
механике
флуида
..........................................................................................
Error! Bookmark not defined.
Древне
цивилизације
..............................................................................................
Error! Bookmark not defined.
Античка
времена
.....................................................................................................
Error! Bookmark not defined.
Ренесанса
.................................................................................................................
Error! Bookmark not defined.
18.
век
(
Хидраулика
)...............................................................................................
Error! Bookmark not defined.
19.
век
(
невискозно
струјање
)................................................................................
Error! Bookmark not defined.
Савремена
механика
флуида
..................................................................................
Error! Bookmark not defined.
Експериментална
механика
флуида
......................................................................
Error! Bookmark not defined.
Рачунарска
динамика
флуида
(
ЦФД
)....................................................................
Error! Bookmark not defined.
Будућност
................................................................................................................
Error! Bookmark not defined.
Примена
механике
флуида
............................................................................................
Error! Bookmark not defined.
Општи
појмови
..................................................................................................................
Error! Bookmark not defined.
Модели
флуида
.......................................................................................................
Error! Bookmark not defined.
Струјни
режими
......................................................................................................
Error! Bookmark not defined.
Врсте
струјања
........................................................................................................
Error! Bookmark not defined.
Реолошка
класификација
.......................................................................................
Error! Bookmark not defined.
Временска
зависност
..............................................................................................
Error! Bookmark not defined.
Састав
флуида
.........................................................................................................
Error! Bookmark not defined.
Физичка
својства
флуида
.................................................................................................
Error! Bookmark not defined.
Молекуларна
грађа
-
микроструктура
........................................................................
Error! Bookmark not defined.
Подела
физичких
својстава
.........................................................................................
Error! Bookmark not defined.
Притисак
..................................................................................................................
Error! Bookmark not defined.
Густина
....................................................................................................................
Error! Bookmark not defined.
Стишљивост
............................................................................................................
Error! Bookmark not defined.
Брзина
звука
............................................................................................................
Error! Bookmark not defined.
Вискозност
...............................................................................................................
Error! Bookmark not defined.
Површински
напон
,
капиларност
и
напон
паре
...................................................
Error! Bookmark not defined.
Кавитација
...............................................................................................................
Error! Bookmark not defined.
Промена
запремине
и
једначина
стања
идеалног
гаса
..............................................
Error! Bookmark not defined.
Једначина
стања
реалних
гасова
............................................................................
Error! Bookmark not defined.
Статика
флуида
.................................................................................................................
Error! Bookmark not defined.
Хидростатички
притисак
............................................................................................
Error! Bookmark not defined.
Ојлерова
једначина
за
миран
флуид
............................................................................
Error! Bookmark not defined.
Распоред
притиска
у
течностима
и
гасовима
у
пољу
Земљине
теже
....................
Error! Bookmark not defined.
Распоред
притиска
у
течностима
...........................................................................
Error! Bookmark not defined.
Промена
притиска
у
гасовима
...............................................................................
Error! Bookmark not defined.
Притисак
течности
на
равне
површине
....................................................................
Error! Bookmark not defined.
Одређивање
нападне
тачке
силе
притиска
...........................................................
Error! Bookmark not defined
Хидростатички
парадокс
........................................................................................
Error! Bookmark not defined.
Притисак
течности
на
криве
површине
.....................................................................
Error! Bookmark not defined.
Пливање
..........................................................................................................................
Error! Bookmark not defined.
Хидростатички
потисак
..........................................................................................
Error! Bookmark not defined.
Стабилност
при
пливању
.......................................................................................
Error! Bookmark not defined.

165
Потенцијално
струјање
...............................................................................................................................................
Значај
модела
потенцијалног
струјања
са
граничним
слојем
...............................................................................
Функција
потенцијала
брзине
φ
и
струјна
функција
ψ
...........................................................................................
Раванска
потенцијална
струјања
добијена
применом
функција
комплексне
променљиве
.................................................................................................................................................................
Комплексни
потенцијали
основних
струјања
.....................................................................................................
Комплексни
потенцијал
облика
степене
функције
.............................................................................................
Струјање
око
кружног
цилиндра
..........................................................................................................................
Једначина
о
промени
количине
кретања
.................................................................................................................
Општи
облик
једначине
о
промени
количине
кретања
..........................................................................................
Примери
примене
једначине
о
промени
количине
кретања
на
непокретне
елементе
.....................................................................................................................................................................
Непокретни
елементи
са
једним
улазом
и
излазом
флуида
...............................................................................
Нагло
проширење
..................................................................................................................................................
Млазник
и
дифузор
са
слободним
истицањем
....................................................................................................
Бордин
наглавак
.....................................................................................................................................................
Отпор
полутела
у
каналу
.......................................................................................................................................
Примена
једначине
о
промени
количине
кретања
на
пропелер
........................................................................
Непокретни
елементи
са
различитим
бројем
улаза
и
излаза
флуида
................................................................
Потисак
ваздушно
-
усисавајућег
погона
..............................................................................................................
Промена
момента
количине
кретања
.....................................................................................................................
Струјање
компресибилног
флуида
............................................................................................................................
Једначина
континуитета
,
брзина
слабог
еластичног
поремећаја
,
Махов
конус
. ...........................................................................................................................................................................
Брзина
слабог
еластичног
поремећаја
–
брзина
звука
........................................................................................
Махов
број
..............................................................................................................................................................
Основни
закони
компресибилног
струјања
..............................................................................................................
Одређивање
критеријума
нестишљивог
струјања
гаса
......................................................................................
Истицање
кроз
конвергентан
млазник
....................................................................................................................
Струјање
кроз
сужени
пресек
..................................................................................................................................
Струјање
кроз
конвергентно
–
дивергентан
млазник
............................................................................................
Закон
промене
брзине
зависно
од
променљиве
површине
проточног
пресека
....................................................................................................................................................................
Динамика
вискозног
флуида
......................................................................................................................................
Ламинарно
струјање
.................................................................................................................................................
Ламинарно
струјање
–
Хаген
-
Поазејево
струјање
кроз
кружну
цев
.................................................................
Турбул
e
нтно
струјање
..............................................................................................................................................
Физичке
величине
и
параметри
у
турбулентном
струјању
................................................................................
Једначине
турбулентног
струјања
течности
(
Рејнолдсове
једначине
)..............................................................
Претпоставке
о
турбулентним
напонима
..............................................................................................................
Прантлова
хипотеза
(Prandtl, 1925) ......................................................................................................................
Профил
брзине
и
пад
притиска
у
кружној
цеви
................................................................................................ 100
Гранични
слој
.............................................................................................................................................................. 106
Једначине
кретања
у
ламинарном
граничном
слоју
.............................................................................................. 111
Једначине
кретања
у
турбулентном
граничном
слоју
......................................................................................... 113
Примена
једначине
о
промени
количине
кретања
................................................................................................ 116
166
Ламинаран
гранични
слој
.................................................................................................................................... 117
Турбулентан
гранични
слој
................................................................................................................................. 119
Одвајање
у
граничном
слоју
.................................................................................................................................... 123
Методе
управљања
граничним
слојем
око
аеропрофила
..................................................................................... 126
Отпори
којима
флуид
дејствује
на
тело
................................................................................................................. 131
Отпор
трења
........................................................................................................................................................... 132
Отпор
облика
,
укупан
отпор
.................................................................................................................................. 134
Отпори
при
малим
и
умереним
R
е
..................................................................................................................... 138
Отпор
аеропрофила
и
сила
узгона
......................................................................................................................... 144
Аеропрофили
........................................................................................................................................................ 144
Теорија
сличности
и
димензијска
анализа
............................................................................................................ 151
Теорија
сличности
и
димензијска
анализа
............................................................................................................ 151
Теорија
сличности
.................................................................................................................................................... 151
Димензијска
анализа
................................................................................................................................................ 158
Литература
................................................................................................................................................................... 162
Садржај
......................................................................................................................................................................... 163

168
САДРЖАЈ
УВОД
У
МЕХАНИКУ
ФЛУИДА
............................
Error! Bookmark not defined.
Историја
механике
флуида
..................................
Error! Bookmark not defined.
Древне
цивилизације
......................................
Error! Bookmark not defined.
Античка
времена
.............................................
Error! Bookmark not defined.
Ренесанса
.........................................................
Error! Bookmark not defined.
18.
век
(
Хидраулика
).......................................
Error! Bookmark not defined.
19.
век
(
невискозно
струјање
)........................
Error! Bookmark not defined.
Савремена
механика
флуида
..........................
Error! Bookmark not defined.
Експериментална
механика
флуида
..............
Error! Bookmark not defined.
Рачунарска
динамика
флуида
(
ЦФД
)............
Error! Bookmark not defined.
Будућност
........................................................
Error! Bookmark not defined.
Примена
механике
флуида
....................................
Error! Bookmark not defined.
ОПШТИ
ПОЈМОВИ
.................................................
Error! Bookmark not defined.
Модели
флуида
...............................................
Error! Bookmark not defined.
Струјни
режими
..............................................
Error! Bookmark not defined.
Врсте
струјања
................................................
Error! Bookmark not defined.
Реолошка
класификација
...............................
Error! Bookmark not defined.
Временска
зависност
......................................
Error! Bookmark not defined.
Састав
флуида
.................................................
Error! Bookmark not defined.
ФИЗИЧКА
СВОЈСТВА
ФЛУИДА
.........................
Error! Bookmark not defined.
Молекуларна
грађа
-
микроструктура
................
Error! Bookmark not defined.
Подела
физичких
својстава
.................................
Error! Bookmark not defined.
Притисак
..........................................................
Error! Bookmark not defined.
Густина
............................................................
Error! Bookmark not defined.
Стишљивост
....................................................
Error! Bookmark not defined.
Брзина
звука
....................................................
Error! Bookmark not defined.
Вискозност
.......................................................
Error! Bookmark not defined.
Површински
напон
,
капиларност
и
напон
паре
.........
Error! Bookmark not
defined.
Кавитација
......................................................
Error! Bookmark not defined.
Основне
термодинамичек
законитости
које
се
користе
у
механици
флуида
.................................................................................
Error! Bookmark not defined.
Једначина
стања
реалних
гасова
....................
Error! Bookmark not defined.
Статика
флуида
.........................................................
Error! Bookmark not defined.
Хидростатички
притисак
....................................
Error! Bookmark not defined.
Ојлерова
једначина
за
миран
флуид
....................
Error! Bookmark not defined.
169
Распоред
притиска
у
течностима
и
гасовима
у
пољу
Земљине
теже
..
Error!
Bookmark not defined.
Распоред
притиска
у
течностима
...................
Error! Bookmark not defined.
Промена
притиска
у
гасовима
.......................
Error! Bookmark not defined.
Притисак
течности
на
равне
површине
............
Error! Bookmark not defined.
Одређивање
нападне
тачке
силе
притиска
...
Error! Bookmark not defined.
Хидростатички
парадокс
................................
Error! Bookmark not defined.
Притисак
течности
на
криве
површине
.............
Error! Bookmark not defined.
Пливање
..................................................................
Error! Bookmark not defined.
Хидростатички
потисак
..................................
Error! Bookmark not defined.
Стабилност
при
пливању
...............................
Error! Bookmark not defined.
Релативно
мировање
течности
..........................
Error! Bookmark not defined.
Релативно
мировање
при
транслацији
..........
Error! Bookmark not defined.
Релативно
мировање
при
ротацији
................
Error! Bookmark not defined.
Кинематика
флуида
и
напонско
стање
............................................................... 1
Лагранжов
метод
.............................................................................................. 2
Ојлеров
метод
................................................................................................... 3
Скаларно
поље
...................................................................................................... 4
Градијент
скаларне
функције
.......................................................................... 4
Векторско
поље
брзине
....................................................................................... 5
Струјница
.......................................................................................................... 5
Струјна
цев
....................................................................................................... 6
Дивергенција
вектора
брзине
.......................................................................... 6
Проток
............................................................................................................... 7
Једначина
континуитета
.................................................................................. 8
Ротор
брзине
..................................................................................................... 9
Вртлог
............................................................................................................. 10
Циркулација
................................................................................................... 13
Класификација
струјних
поља
...................................................................... 14
Векторско
поље
убрзања
................................................................................... 16
Тензорско
поље
................................................................................................... 16
Тензор
брзина
деформисања
делића
............................................................ 16
Тензор
напона
–
напонско
стање
.................................................................. 18
Веза
напона
са
брзинама
деформисања
....................................................... 19
Притисак
......................................................................................................... 21
Динамика
идеалног
флуида
.....................................
Error! Bookmark not defined.
O
јлерова
једначина
................................................
Error! Bookmark not defined.
Бернулијев
интеграл
Ојлерове
једначине
............
Error! Bookmark not defined.

Ovaj materijal je namenjen za učenje i pripremu, ne za predaju.
Slični dokumenti