Методологија интегралне анализе и оптимизације аеродинамичких површина ваздухопловних конструкција
др Јелена М. Сворцан — данас редовна професорка на Машинском факултету Универзитета у Београду, Катедра за ваздухопловство (према страници факултета)
Универзитет у Београду, Машински факултет; 2014 (Комисија образована 8. маја 2014, извештај Комисије 19. маја 2014, усвојен на Наставно-научном већу 22. маја 2014, сагласност Већа научних области техничких наука 7. јула 2014, одбрана 16. јула 2014); 135 нумерисаних страна с прилозима
Пун текст рада у Националном репозиторијуму дисертација у Србији (НаРДуС); извештај Комисије у материјалима седнице Већа научних области техничких наука Универзитета у Београду. Осврт: „Стручни критички осврт на дисертацију „Методологија интегралне анализе и оптимизације аеродинамичких површина ваздухопловних конструкција““, 23. септембар 2026., 27 страна, с пописом литературе.
Дисертација гради затворен и аутоматизован поступак од параметарске геометрије преко мреже и прорачуна струјања до стохастичке оптимизације и испробава га на три примера: малој ветротурбини с вертикалном осом, окозвучној каскади профила с активним млазевима и неуронским мрежама и крилу ОНЕРА М6 које маше, пропиње се, савија и увија. Осврт признаје изведен поступак и четири праве методолошке новине, међу њима спрегу нестационарног прорачуна с једначином обртања ротора, али налази да је поставка ветротурбине подешавана према истим мерењима према којима се проверава и затим скалирана поправним фактором, што је калибрација, а не валидација; да су студије мреже и временског корака описане само речима; да су све три оптимизације изведене по једанпут без наведених параметара; да проширења нису проверена мерењима, а тачност неуронских мрежа мерена је на скупу за тренирање; и да закључак изоставља ограде које главе рада саме изричу.
Основне примедбе
| # | Примедба | Шта рад каже или ради | У чему је проблем | Стране | Мерило |
|---|---|---|---|---|---|
| 1 | Калибрација названа валидацијом | Рад каже: „Све одлуке о прорачунској мрежи и нумеричком моделу донете су након опсежног испитивања, а на основу доброг или лошег поклапања са доступним експерименталним подацима.“ (стр. 44), и то понавља на стр. 53, а затим израчунату криву скалира поправним фактором изведеним из односа максимума. | Поставка није независна од мерења с којима се пореди, а фактор по конструкцији изједначава висину максимума рачуна и мерења, па поређење проверава само облик криве и положај максимума; поређење на сликама 5.7 и 5.15 има карактер калибрације и провере сагласности, не независне валидације. | 44, 45, 53 | Oberkampf и Trucano 2002 (Sandia, SAND2002-0529); уређивачка политика часописа Journal of Fluids Engineering о нумеричкој тачности |
| 2 | Верификација прорачуна описана само речима | Студије мреже, броја итерација и временске поделе описане су без бројева — „Опет, изведена је одговарајућа студија и закључено да се може усвојити нижа вредност броја итерација за довољно фину временску поделу.“ (стр. 53) — а најављено поређење структуриране и неструктуриране мреже у наставку не постоји (стр. 52). | Дискретизациона грешка није процењена ни на један од поступака објављених пре рока, па се не може рећи да ли су мале разлике међу решењима у оптимизацијама (реда 1 %) веће од нумеричке грешке. | 52, 53 | уређивачка политика Journal of Fluids Engineering (критеријуми 4–6); Celik и др. 2008; водич NASA/NPARC |
| 3 | Оптимизације изведене по једанпут, без наведених параметара | Рој честица примењен је с 15 честица у 12 корака, а вредности коефицијената нису наведене; рад то образлаже: „У том смислу, коефицијенти који се појављују можда нису најбоље подешени“ (стр. 57) и „Основно ограничење у овој студији је дужина укупног времена прорачуна па број изведених итерација није велик.“ (стр. 60). | Без параметара, поновљених покретања из различитих семена и бројчаног исхода прве оптимизације, ток на слици 5.21 не може се разликовати од раног застоја роја, а оптимуми рада јесу оптимуми усвојеног модела у једном покретању. | 57, 58, 59, 60 | Barr и др. 1995 (смернице за извештавање о хеуристичким експериментима); Eiben и Jelasity 2002 |
| 4 | Симулација теренског рада без момента генератора | Једначина обртања ротора садржи само аеродинамички момент, а генератор је сведен на горњу границу брзине: „Угаона брзина може да расте до неке, унапред задате номиналне вредности (која је у реалности одређена генератором и/или редуктором)“ (стр. 70); закључак ипак каже: „Развијена је нова методологија анализе ветротурбине којом се врши симулација њеног реалног рада на терену при изразито променљивим брзинама ветра.“ (стр. 124). | Ротор зато убрзава ка слободном ходу, док мерена турбина у Упсали, оптерећена генератором, ради у другом опсегу; поређење с мерењем је само визуелно („Може се приметити да приказани графици сада више подсећају на резултате сакупљене током експеримента“, стр. 71), уз тренутне нумеричке вредности наспрам десетоминутних просека. | 68, 70, 71, 73, 124 | аргумент из грађе самог рада; Scheurich и Brown 2013; Danao 2012 |
| 5 | Тачност неуронских мрежа мерена на скупу за тренирање | Рад отворено каже: „Поређење различитих НМ вршено је над скупом за проверу, једнаком скупу за тренирање“ (стр. 90). | Вероватноће у табели 6.2 зато мере прилагођење мреже подацима, а не способност предвиђања; правила за мерење перформанси мрежа траже одвојен скуп за проверу, а за мале скупове постојала је унакрсна провера која не би захтевала нове прорачуне. | 90, 91 | Prechelt 1994 (PROBEN1, § 2.3); Hastie, Tibshirani и Friedman 2009 |
| 6 | Увијање крила с поларним моментом уместо торзионе константе | Увојна крутост у изразу (7.8) узета је као производ модула клизања и поларног момента инерције, а референтна оса изабрана је по минимуму поларног момента: „Ова локација одабрана је на основу минималне вредности поларног момента инерције аеропрофила“ (стр. 117); модул еластичности задат је као „нешто нижа него за алуминијум“ (стр. 115) уз вредност око четрнаест пута мању. | Производ модула клизања и поларног момента тачан је само за кружни пресек; за пуни профил торзиона константа је око 25 пута мања, па су бројчани резултати увијања ефективне, а не физичке крутости, а оса кроз тежиште није еластична оса. | 115, 117 | теорија Сен-Венанове торзије (Lagace 2001) |
| 7 | Закључак изоставља ограде из глава рада | Закључак говори о „поправног фактора којим је могуће кориговати резултате раванске анализе да би се добиле реалне, тродимензионалне, употребљиве вредности“ (стр. 124), док глава 5 о истој зависности каже: „прилично је тешко извести тачне, дефинитивне закључке“ (стр. 56). | Површ поправног фактора провучена је кроз седам тачака разнородних мерења и делимично екстраполисана, а функција циља оптимизације од ње зависи; закључак и доприноси ту ограду, као и ограде са стр. 106 и 117, не понављају. | 56, 68, 124 | аргумент из грађе самог рада |
Шта рад доноси
- Стварно изведен, затворен поступак: „Основна предност оваквог приступа је што је поступак моделирања или прорачуна могуће поновити произвољан број пута што представља потребан услов за спровођење оптимизације одабраних параметара.“ (стр. 32), с мрежама исте топологије за све конфигурације („Тако је постигнута условна једнакост мрежа“, стр. 52).
- Спрега нестационарног прорачуна струјања с једначином обртања ротора, у којој се угаона брзина рачуна из аеродинамичког момента: „Нова вредност угаоне брзине ротора рачуна се након сваког временског корака.“ (стр. 70), с приложеним кодом и тачно омеђеном новином према литератури (стр. 68).
- Сопствена реализација вишекритеријумске оптимизације ројем честица над 320 прорачуна струјања, с другим критеријумом који има инжењерски смисао (стр. 62) и исправно схваћеним Парето приступом: „Према томе, коначно одабрани дизајн мале ветротурбине са вертикалном осом обртања представља избор аутора.“ (стр. 67).
- Мерљив критеријум одвајања и активације млазева у каскади: „Млазеви се активирају само уколико долази до коначног одвајања струјања и то почевши од млаза који се налази низводно од тачке одвајања.“ (стр. 82).
- Изричито исказане границе сопствених резултата на много места, од занемарених утицаја пре поређења с мерењима (стр. 44) до неповољног исхода спрегнуте анализе, у којој је уз повећање узгона од 20 % отпор порастао за 44 % (стр. 119).
Образложење примедби
- Калибрација названа валидацијомРад поређење с мерењима јесте урадио на свим основним поставкама и занемарене утицаје унапред набројао (стр. 44); замерка је у имену поступка, не у самом поређењу.
- Верификација прорачуна описана само речимаРад компромис сам изриче: „Треба напоменути да је тренутни циљ био спровести процес оптимизације и да је то морало бити реализовано на мали уштрб тачности и у складу са доступним рачунарским ресурсима.“ (стр. 52).
- Оптимизације изведене по једанпут, без наведених параметараСкупа функција циља оправдава мали рој и мали број корака, што исте смернице допуштају; не оправдава изостанак параметара.
- Симулација теренског рада без момента генератораУгаона брзина ротора израчуната из једначине кретања јесте стварна разлика у односу на тада објављене нумеричке студије, а сам рад у глави 5 говори скромније: „нешто другачија методологија нумеричке процене перформанси ветротурбине која подражава њен реални рад на терену“ (стр. 68).
- Тачност неуронских мрежа мерена на скупу за тренирањеСпособност предвиђања проверена је на две геометрије ван скупа за тренирање, што рад наводи, а при углу изван наученог опсега поклапање је слабије, што рад такође каже.
- Увијање крила с поларним моментом уместо торзионе константеСмер увијања и квалитативни налази нису погођени, а рад упрошћења изричито набраја (стр. 121).
- Закључак изоставља ограде из глава радаОно што треба изменити у највећем је делу већ у раду тачно речено, у главама, а не у закључку.
Осврт не доводи у питање оригиналност рада и не тврди повреду академске честитости. Реч је о стручној критици тачности, метода и закључивања, с дословним наводима и странама дисертације. Као мерило су коришћени само литература, стандарди и прописи који су постојали пре предаје рукописа комисији (почетком маја 2014).