MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА...
Transcript of MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА...
![Page 1: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/1.jpg)
ГРАЂЕВИНСКИ ФАКУЛТЕТ
УНИВЕРЗИТЕТА У БЕОГРАДУКАТЕДРА ЗА ХИДРОТЕХНИКУ И
ВОДНО ЕКОЛОШКО ИНЖЕЊЕРСТВО
ДОМАЋИ РАДИЗ УЖЕ НАУЧНЕ ОБЛАСТИ: МЕХАНИКА ФЛУИДА
MОДЕЛИРАЊЕ ТУРБУЛЕНЦИЈЕ НА ПРИМЈЕРУ РИБЉЕ СТАЗЕВАРИЈАНТА 2
МЕНТОРИ:
ПРОФ. ДР ДУШАН ПРОДАНОВИЋ
ДОЦЕНТ ДР ДАМЈАН ИВЕТИЋ
СТУДЕНТ:
МИЛИЦА ТРИФКОВИЋ
Фебруар, 2020.
![Page 2: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/2.jpg)
САДРЖАЈ
1. ОПИС ПРОБЛЕМА
2. ПОСТУПАК ИЗРАДЕ МОДЕЛА (ОДАБИР ПАРАМЕТАРА)
3. РЕЗУЛТАТИ
4. ЗАКЉУЧЦИ
iRIC Nays Cube
Win Flume
![Page 3: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/3.jpg)
ОПИС ПРОБЛЕМАРазвој турбуленције у рибљој
стази датих геометријских
карактеристика
Стаза се састоји од:
• прилазне дионице
• сужења (препрека) и
• низводног базена
![Page 4: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/4.jpg)
ПОСТУПАК ИЗРАДЕ МОДЕЛА
Састоји се у неколико корака приказаних алгоритмом лијево
![Page 5: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/5.jpg)
ПОСТУПАК ИЗРАДЕ МОДЕЛА
Састоји се у неколико корака приказаних алгоритмом лијево
![Page 6: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/6.jpg)
ПОСТУПАК ИЗРАДЕ МОДЕЛА
![Page 7: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/7.jpg)
ПОСТУПАК ИЗРАДЕ МОДЕЛА
НА ТАЧНОСТ РЕЗУЛТАТА УТИЧЕ :
❖ промјена густине мреже (геометрија канала)
❖ дубина низводно (гранични услов)
❖ временски корак прорачуна
❖ положај препреке у каналу ...
![Page 8: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/8.jpg)
РЕЗУЛТАТИ
ПРОМЈЕНА НИВОА / ДУБИНЕ ВОДЕ
iRIC Nays Cube
Т = 0 s Т = 30 s Т = 60 s
Предност програма - приказ резултата у времену и простору
![Page 9: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/9.jpg)
РЕЗУЛТАТИ
ПРОМЈЕНА НИВОА / ДУБИНЕ ВОДЕ
iRIC Nays Cube
Распоред дубина у сужењу (лијево) за модификовану скалу
![Page 10: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/10.jpg)
РЕЗУЛТАТИ
ПРОМЈЕНА НИВОА / ДУБИНЕ ВОДЕ
iRIC Nays Cube
Дубина воде у каналу се устали већ након 40s
![Page 11: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/11.jpg)
РЕЗУЛТАТИ
ПРОМЈЕНА БРЗИНЕ
• ПО ШИРИНИ КАНАЛА
iRIC Nays Cube
![Page 12: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/12.jpg)
РЕЗУЛТАТИ
ПРОМЈЕНА БРЗИНЕ
• ПО ДУЖИНИ КАНАЛА
iRIC Nays Cube
![Page 13: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/13.jpg)
РЕЗУЛТАТИ
ПРОМЈЕНА БРЗИНЕ
• ПО ДУБИНИ КАНАЛА
iRIC Nays Cube
![Page 14: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/14.jpg)
РЕЗУЛТАТИ
ПРОМЈЕНА СТРУЈНЕ СЛИКЕ (СТРУЈНИЦЕ)
iRIC Nays Cube
Т = 0 s Т = 30 s Т = 60 s
![Page 15: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/15.jpg)
РЕЗУЛТАТИ
ПОРЕЂЕЊЕ РЕЗУЛТАТА
Win Flume
• Формирана геометрија
КОНТРОЛНИ ПРЕСЈЕК - СУЖЕЊЕ
![Page 16: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/16.jpg)
РЕЗУЛТАТИПОРЕЂЕЊЕ РЕЗУЛТАТА 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… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/17.jpg)
РЕЗУЛТАТИПОРЕЂЕЊЕ РЕЗУЛТАТА 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… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/18.jpg)
РЕЗУЛТАТИ
ПРОЦЈЕНА КОЕФИЦИЈЕНТА ЛОКАЛНОГ
ГУБИТКА НА СУЖЕЊУ
На основу стечених знања из области механике
флуида и хидраулике отворених токова, закључује
се да се у сужењу, које је један вид прага, јавља
критична дубина.
Дакле, из услова да је Фрудов број једнак јединици
(Fr =1), у пресјеку 2-2 се добија критична дубина
hkr = 0.162 m.
![Page 19: MОДЕЛИРАЊЕТУРБУЛЕНЦИЈЕ НА ...hikom.grf.bg.ac.rs/stari-sajt/web_stranice/KatZaHidr/Predmeti/MehFl… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/19.jpg)
РЕЗУЛТАТИ
ПРОЦЈЕНА КОЕФИЦИЈЕНТА ЛОКАЛНОГ
ГУБИТКА НА СУЖЕЊУ
Е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… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/20.jpg)
ЗАКЉУЧЦИ
❖ Брзине у каналу се крећу у опсегу од 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… · ЗАКЉУЧАК - водити рачуна о мјесту](https://reader036.fdocument.pub/reader036/viewer/2022071223/60855545d027f656ca7f5eea/html5/thumbnails/21.jpg)
ГРАЂЕВИНСКИ ФАКУЛТЕТ
УНИВЕРЗИТЕТА У БЕОГРАДУКАТЕДРА ЗА ХИДРОТЕХНИКУ И
ВОДНО ЕКОЛОШКО ИНЖЕЊЕРСТВО
ХВАЛА НА ПАЖЊИMОДЕЛИРАЊЕ ТУРБУЛЕНЦИЈЕ НА ПРИМЈЕРУ РИБЉЕ СТАЗЕ
ВАРИЈАНТА 2
МЕНТОРИ:
ПРОФ. ДР ДУШАН ПРОДАНОВИЋ
ДОЦЕНТ ДР ДАМЈАН ИВЕТИЋ
СТУДЕНТ:
МИЛИЦА ТРИФКОВИЋ
Фебруар, 2020.