MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА...

21
ГРАЂЕВИНСКИ ФАКУЛТЕТ УНИВЕРЗИТЕТА У БЕОГРАДУ КАТЕДРА ЗА ХИДРОТЕХНИКУ И ВОДНО ЕКОЛОШКО ИНЖЕЊЕРСТВО ДОМАЋИ РАД ИЗ УЖЕ НАУЧНЕ ОБЛАСТИ: МЕХАНИКА ФЛУИДА MОДЕЛИРАЊЕ ТУРБУЛЕНЦИЈЕ НА ПРИМЈЕРУ РИБЉЕ СТАЗЕ ВАРИЈАНТА 2 МЕНТОРИ: ПРОФ. ДР ДУШАН ПРОДАНОВИЋ ДОЦЕНТ ДР ДАМЈАН ИВЕТИЋ СТУДЕНТ: МИЛИЦА ТРИФКОВИЋ Фебруар, 2020.

Transcript of MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА...

Page 1: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

ГРАЂЕВИНСКИ ФАКУЛТЕТ

УНИВЕРЗИТЕТА У БЕОГРАДУКАТЕДРА ЗА ХИДРОТЕХНИКУ И

ВОДНО ЕКОЛОШКО ИНЖЕЊЕРСТВО

ДОМАЋИ РАДИЗ УЖЕ НАУЧНЕ ОБЛАСТИ: МЕХАНИКА ФЛУИДА

MОДЕЛИРАЊЕ ТУРБУЛЕНЦИЈЕ НА ПРИМЈЕРУ РИБЉЕ СТАЗЕВАРИЈАНТА 2

МЕНТОРИ:

ПРОФ. ДР ДУШАН ПРОДАНОВИЋ

ДОЦЕНТ ДР ДАМЈАН ИВЕТИЋ

СТУДЕНТ:

МИЛИЦА ТРИФКОВИЋ

Фебруар, 2020.

Page 2: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

САДРЖАЈ

1. ОПИС ПРОБЛЕМА

2. ПОСТУПАК ИЗРАДЕ МОДЕЛА (ОДАБИР ПАРАМЕТАРА)

3. РЕЗУЛТАТИ

4. ЗАКЉУЧЦИ

iRIC Nays Cube

Win Flume

Page 3: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

ОПИС ПРОБЛЕМАРазвој турбуленције у рибљој

стази датих геометријских

карактеристика

Стаза се састоји од:

• прилазне дионице

• сужења (препрека) и

• низводног базена

Page 4: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

ПОСТУПАК ИЗРАДЕ МОДЕЛА

Састоји се у неколико корака приказаних алгоритмом лијево

Page 5: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

ПОСТУПАК ИЗРАДЕ МОДЕЛА

Састоји се у неколико корака приказаних алгоритмом лијево

Page 6: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

ПОСТУПАК ИЗРАДЕ МОДЕЛА

Page 7: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

ПОСТУПАК ИЗРАДЕ МОДЕЛА

НА ТАЧНОСТ РЕЗУЛТАТА УТИЧЕ :

❖ промјена густине мреже (геометрија канала)

❖ дубина низводно (гранични услов)

❖ временски корак прорачуна

❖ положај препреке у каналу ...

Page 8: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

РЕЗУЛТАТИ

ПРОМЈЕНА НИВОА / ДУБИНЕ ВОДЕ

iRIC Nays Cube

Т = 0 s Т = 30 s Т = 60 s

Предност програма - приказ резултата у времену и простору

Page 9: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

РЕЗУЛТАТИ

ПРОМЈЕНА НИВОА / ДУБИНЕ ВОДЕ

iRIC Nays Cube

Распоред дубина у сужењу (лијево) за модификовану скалу

Page 10: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

РЕЗУЛТАТИ

ПРОМЈЕНА НИВОА / ДУБИНЕ ВОДЕ

iRIC Nays Cube

Дубина воде у каналу се устали већ након 40s

Page 11: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

РЕЗУЛТАТИ

ПРОМЈЕНА БРЗИНЕ

• ПО ШИРИНИ КАНАЛА

iRIC Nays Cube

Page 12: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

РЕЗУЛТАТИ

ПРОМЈЕНА БРЗИНЕ

• ПО ДУЖИНИ КАНАЛА

iRIC Nays Cube

Page 13: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

РЕЗУЛТАТИ

ПРОМЈЕНА БРЗИНЕ

• ПО ДУБИНИ КАНАЛА

iRIC Nays Cube

Page 14: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

РЕЗУЛТАТИ

ПРОМЈЕНА СТРУЈНЕ СЛИКЕ (СТРУЈНИЦЕ)

iRIC Nays Cube

Т = 0 s Т = 30 s Т = 60 s

Page 15: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

РЕЗУЛТАТИ

ПОРЕЂЕЊЕ РЕЗУЛТАТА

Win Flume

• Формирана геометрија

КОНТРОЛНИ ПРЕСЈЕК - СУЖЕЊЕ

Page 16: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

РЕЗУЛТАТИПОРЕЂЕЊЕ РЕЗУЛТАТА Win Flume

Након формирања геометрије

у програму Win Flume слиједи

прорачун КРИВЕ ПРОТОКА

За тачку непосредно испред

препреке у програму iRIC Nays

Cube је очитана дубина:

X Depth

30 0.233712

Уколико се мјерење дубине врши непосредно испред

препреке, неслагање резултата је неизбјежно!!!

0.15

0.17

0.19

0.21

0.23

0.25

0.07 0.08 0.09 0.1 0.11 0.12 0.13 0.14 0.15

H [

m]

Q [m3/s]

КРИВА ПРОТОКА

iRIC Nays Cube

Win Flume

Page 17: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

РЕЗУЛТАТИПОРЕЂЕЊЕ РЕЗУЛТАТА Win Flume

За тачку довољно узводно од

препреке (не мјесту гдје су

струјнице праве и паралелне) у

програму iRIC Nays Cube је

очитана дубина:

X Depth

30 0.243112

Уколико се мјерење дубине врши узводно од препреке

(не мјесту гдје су струјнице праве и паралелне) , добија

се одлично слагање резултата у наведена два програма!!!

0.15

0.17

0.19

0.21

0.23

0.25

0.07 0.08 0.09 0.1 0.11 0.12 0.13 0.14 0.15

H [

m]

Q [m3/s]

КРИВА ПРОТОКА

iRIC Nays Cube

Win Flume

ЗАКЉУЧАК - водити рачуна о мјесту мјерења дубине

Page 18: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

РЕЗУЛТАТИ

ПРОЦЈЕНА КОЕФИЦИЈЕНТА ЛОКАЛНОГ

ГУБИТКА НА СУЖЕЊУ

На основу стечених знања из области механике

флуида и хидраулике отворених токова, закључује

се да се у сужењу, које је један вид прага, јавља

критична дубина.

Дакле, из услова да је Фрудов број једнак јединици

(Fr =1), у пресјеку 2-2 се добија критична дубина

hkr = 0.162 m.

Page 19: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

РЕЗУЛТАТИ

ПРОЦЈЕНА КОЕФИЦИЈЕНТА ЛОКАЛНОГ

ГУБИТКА НА СУЖЕЊУ

Е1 = Е2 + Δ𝐸1−2

П1 +𝑣12

2𝑔= ℎ𝑘𝑟 +

𝑣22

2𝑔+ 𝜉𝑠𝑢ž

𝑣22

2𝑔

0.32 +0.22

2𝑔= 0.162 +

1.12

2𝑔+ 𝜉𝑠𝑢ž

1.12

2𝑔

𝜉𝑠𝑢ž = 1.594

НАПОМЕНА : ради се о грубој процјени

коефицијента локалног губитка на сужењу

јер параметри брзина и очитане дубине

узводно уносе дозу неизвјесности и грешке

у резултат.

Page 20: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

ЗАКЉУЧЦИ

❖ Брзине у каналу се крећу у опсегу од 0.00 до 2.95 m/s

❖ Прорачун је спроведен за висину препреке 20 cm

❖ Долази до појаве вртложног струјања непосредно узводно од препреке и у базену низводно

❖ Повољније рјешење са степеничастом препреком - лакши прелаз за рибе

❖ Постигнуто добро слагање прорачуна у програму Win Flume и iRIC Nays Cube

❖ Тачно рјешење не постоји - само његова варијанта

Page 21: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту

ГРАЂЕВИНСКИ ФАКУЛТЕТ

УНИВЕРЗИТЕТА У БЕОГРАДУКАТЕДРА ЗА ХИДРОТЕХНИКУ И

ВОДНО ЕКОЛОШКО ИНЖЕЊЕРСТВО

ХВАЛА НА ПАЖЊИMОДЕЛИРАЊЕ ТУРБУЛЕНЦИЈЕ НА ПРИМЈЕРУ РИБЉЕ СТАЗЕ

ВАРИЈАНТА 2

МЕНТОРИ:

ПРОФ. ДР ДУШАН ПРОДАНОВИЋ

ДОЦЕНТ ДР ДАМЈАН ИВЕТИЋ

СТУДЕНТ:

МИЛИЦА ТРИФКОВИЋ

Фебруар, 2020.