докторска дисертација · хемијско инжењерство, Фишер–Тропшова синтеза, кинетичко моделирање, реакторско инжењерство, оптимизација процеса
Kinetic modeling and optimization of fixed-bed reactor for Fischer-Tropsch synthesis
др Бранислав С. Тодић — у време предаје рада докторанд Технолошко-металуршког факултета Универзитета у Београду (докторске студије уписане школске 2011/12, према биографији у раду); данас доцент на истом факултету (према јавној страници факултета)
Универзитет у Београду, Технолошко-металуршки факултет; 2015 (комисија за оцену и одбрану именована одлуком бр. 35/73 од 26. фебруара 2015, реферат од 22. априла 2015); упоредни српски наслов: Modelovanje hemijske kinetike i optimizacija reaktora sa pakovanim slojem za Fischer-Tropsch sintezu; рад је писан на енглеском језику; 4 предстране без броја + 219 страна, садржај и спискови на стр. 7–10, основни текст на странама 10–184, литература на стр. 185–199, спискови табела, слика и симбола на стр. 200–210, списак радова аутора од стр. 211, биографија на стр. 214–215, изјаве на стр. 216–219
Пун текст рада и реферат комисије у Националном репозиторијуму дисертација у Србији (НаРДуС). Осврт: „Стручни критички осврт на дисертацију „Kinetic modeling and optimization of fixed-bed reactor for Fischer-Tropsch synthesis““, 22. септембар 2026., 15 страна, с пописом литературе.
Дисертација узима за предмет Фишер–Тропшову синтезу у два њена најзахтевнија слоја, расподелу производа и кинетику: анализира утицај процесних услова на активност и селективност два катализатора, из механистичких претпоставки изводи низ детаљних Langmuir–Hinshelwood–Hougen–Watson модела кинетике, оцењује им параметре и међу њима бира, а изабрани модел спреже с једнодимензионим моделом цевног реактора с пакованим слојем и оптимизује улазне услове генетским алгоритмом. Сопствени експеримент није део тезе, што рад сам поштено изриче; аутор ради обраду, моделирање, симулацију и оптимизацију туђих лабораторијских података. Осврт налази стваран допринос кинетичком моделирању (концепт десорпције 1-олефина зависне од броја угљеникових атома, дискриминација модела по физичко-хемијским ограничењима, анализа носиоца промене селективности) и уредан научни апарат, а као слабости ослонац на један скуп података без провере на издвојеним подацима, екстраполацију кинетике на крупне честице уз претпоставку потпуне ефикасности катализатора, три неусаглашена скупа бројева за оптимизоване услове и низ уредничких омашки. Оцена: солидан и добро изведен развојно-примењени рад, уз препоруку за наведене измене у евентуалном монографском издању.
Основне примедбе
| # | Примедба | Шта рад каже или ради | У чему је проблем | Стране | Мерило |
|---|---|---|---|---|---|
| 1 | Оцена и избор модела на једном скупу података, без провере на издвојеним подацима | Слагање је приказано на подацима коришћеним за оцену параметара, без провере на издвојеном подскупу или другом катализатору и реактору; на то се надовезује реченица о применљивости модела на „all FTS reactor types“ (стр. 161), коју сам рад у истом ставу и у закључку своди на почетни дизајн. | Реч је о слагању унутар узорка; провера на независним подацима мерило је добре праксе кинетичког моделирања, па ограниченост скупа треба навести и тврдњу о применљивости доследно свести на почетни дизајн. | 161 | Froment, Bischoff, De Wilde 2011; EUROKIN (Berger и др.) 2001 |
| 2 | Кинетика оцењена на ситним честицама примењена на цевни реактор с честицама 1–3 mm | Кинетика оцењена на суспензионом реактору с честицама 44–90 µm примењена је на цевни реактор с честицама 1–3 mm уз занемарен интрачестични пренос, што рад сам ограничава: „this assumptions is justified only for very small catalyst particles“ (стр. 163); уз то „our model assumes 100 % catalyst effectiveness, which is not likely to be the case“ (стр. 170). | Екстраполација ван опсега важења претпоставке; требало би увести модел интрачестичног преноса и валидацију спрам мерених аксијалних профила, чиме би поглавље 5 прешло с процене реда величине на проверену предикцију. | 163, 170 | Froment, Bischoff, De Wilde 2011; Rawlings, Ekerdt 2012 |
| 3 | Оптимизовани услови наведени на три места с три неусаглашена скупа бројева | У тексту „The optimized inlet temperature was 477.3 K, feed H2/CO ratio 1.83 and inlet syngas flowrate 0.156 mol/s.“ (стр. 178), у Табели 5.4 однос 1,82 и проток 0,152 mol/s, а потпис слике 5.10 даје однос H2/CO = 2,13, што противречи изричитом закључку: „the optimal feed reactant ratio was lower compared to base case with 1.83 compared to 2.0“ (стр. 178). | Омашка при препису потписа слике нарушава поверење у пренети оптимум; три вредности треба ускладити, а потпис слике 5.10 исправити. | 178 | унутрашња сагласност рада |
| 4 | Снага статистичке дискриминације модела прецењена, без интервала поверења по параметру | Пријављени број података „was 696 (24 mass balances with 29 responses)“ (стр. 139) прецењује количину независне информације, јер су одзиви унутар истог биланса међузависни; „Weighting factors were not used in the objective function, because of lack of availability of sufficient number of replicate balances“ (стр. 137); за моделе с деветнаест, односно двадесет један параметром не наводе се интервали поверења ни матрица корелације, а бројилачки степен слободе F-теста код модела на CO-уметању наведен је као 11 уместо 21. | F-тест и t-тест оптерећени су оптимистичким бројем степени слободе; интервали поверења и корелација параметара стандардни су елемент исказивања нелинеарне оцене и треба их надоместити. | 137, 139 | Bard 1974; Seber, Wild 1989; EUROKIN (Berger и др.) 2001 |
| 5 | Пад притиска виши од очекиваног усвојен као активно ограничење оптимизације | „the pressure drop within the bed was higher than expected“ (стр. 170), што аутор приписује избору корелација и облику честице, али се иста вредност усваја као ограничење: „Maximum allowable pressure drop is set to 13 bar“ (стр. 176), док дијаграм алгоритма (Слика 5.9) наводи 15 bar; оптимум лежи на активним ограничењима (конверзија 60,0%, пад притиска 13,0 bar), а добитак у продуктивности C5+ према базном случају је око 9%. | Граница почива на вероватно прецењеном излазу модела и треба је преиспитати; резултат оптимизације претежно је одређен изабраним границама, што ваља изричито исказати уз основицу поређења. | 170, 176 | Edgar, Himmelblau, Lasdon 2001; Goldberg 1989 |
| 6 | Провера реакторског модела само поређењем реда величине | Провера се своди на поређење годишње продуктивности при неистоврсним основицама: „This is comparable in the order of magnitude with the industrial plant 144 000 tonnes/y“ (стр. 170), без поређења с мереним аксијалним профилима. | Тврдњу да модел даје резултате репрезентативне за велике индустријске реакторе треба умерити на ниво процене реда величине. | 169, 170 | Froment, Bischoff, De Wilde 2011; Rawlings, Ekerdt 2012 |
| 7 | Уредничке омашке поглавља 5: ознаке једначина, упуте и опис нумеричке методе | Ознака једначине (5.19) додељена је двема једначинама на стр. 168; потпис слике 5.2 упућује на Табелу 5.2 уместо 5.3; група C5+ представљена је недоследно (октан, C9, тежак n-парафин); метода нумеричког решавања именована је непрецизно, јер је коришћен солвер ode15s (имплицитни вишекорачни солвер за круте системе), а у истом реду помиње се и Runge–Kutta метода, уз словну омашку у називу; знак реакционог члана у билансима масе и топлоте користи недоследну конвенцију предзнака. | Питања интерне доследности која отежавају праћење носећег резултата; стварно коришћени солвер је исправан и примерен, па грешка погађа само опис. | 168 | Универзитет у Београду, Упутство, Прилог 4, 2011; Shampine, Reichelt 1997 |
| 8 | Погрешна референца уз опис гвозденог огледа и омашке у називу извора података | Референца уз опис припреме и реакторског система за гвоздени оглед упућује на рад о кобалтном катализатору; у називу установе која је извор кобалтних података провлаче се две словне грешке — „Center for Advanced Energy Research, University of Kenticky“ (стр. 65) — уз омашке у именима (Satterfield, Kellogg, Nowicki); притисак се исказује час у MPa, час у bar; наведени максимум селективности према CO2 од 1,2% незнатно подцењује вредност из Табеле 3.2 (1,25%). | Порекло гвоздених огледа тачно је наведено на другом месту, па је реч о погрешном броју референце, а не о непроверљивом методу; назив извора података треба исправити, а јединице уједначити. | 65 | BIPM, SI брошура, 8. изд., 2006; The ACS Style Guide, 3. изд., 2006 |
| 9 | Преоштре и прекомерно уопштене формулације у поглављу 3 | Поређење с Dinse и др. постављено је као супротност — „In contrast to the above mentioned explanation, our results, obtained at higher pressures“ (стр. 108) — иако се опсези притиска не преклапају; закључци поглавља уопштавају на кобалтни катализатор уопште, а сви подаци потичу с једног катализатора, једне шарже и једног реактора, што тело текста ограничава на „over the examined catalyst and conditions used“ (стр. 126); уз то стоји штампарска омашка 473 K где смисаоно треба 478 K: „three temperatures (478, 493, 503 K)“ (стр. 99). | Реч је пре о разлици или засићењу у опсегу него о правој противречности, па формулацију треба ублажити, а уопштавајуће исказе доследно ограничити на испитани катализатор и услове. | 99, 108, 126 | унутрашња сагласност рада |
| 10 | Омашке у приказу преноса масе и терминологији | „where later diffuses faster“ (стр. 56), иако водоник, мања молекула, дифундује брже од CO, што потврђује наредна реченица истог става; праг преноса масе отиснут је као „Catalyst pellets with diameter lower than 200 nm do not experience mass transfer issues“ (стр. 56), где смисаоно треба 200 µm, како рад сам даје на другом месту: „The particles larger than 0.2 mm are known to exhibit mass transfer difficulties“ (стр. 50); „density function theory“ (стр. 24, 130) уместо functional theory. | Јединична и терминолошка омашка без последица по резултате, јер реакторски модел претпоставља потпуну ефикасност катализатора; треба их исправити. | 24, 50, 56, 130 | унутрашња сагласност рада; The ACS Style Guide, 3. изд., 2006 |
| 11 | Научни апарат: дуплиране јединице, непотпуне референце, изостанак кода и интервала поверења | Најмање једанаест група јединица наводи исти извор под два (местимично три) броја, па списак од 387 прецењује стваран број извора за једанаест-дванаест; неколико јединица непотпуно је или малформисано; програмски код (генетски алгоритам, Levenberg–Marquardt, ode15s, оптимизација реактора) није приложен, дати су само дијаграми тока; сопствени рад аутора из 2015. на коме почива део поглавља 2 није унет у нумерисани списак. | Дупликате треба објединити, непотпуне јединице допунити, а уз оцењене параметре табелирати интервале поверења и матрицу корелације; изостанак кода отежава дословно понављање рачуна и треба га надоместити. | 185–199 | The ACS Style Guide, 3. изд., 2006, гл. 14; Универзитет у Београду, Упутство, Прилог 4, 2011; ISO 690:2010 |
| 12 | Језичке омашке у енглеском тексту | „pseudohomogenious“ (стр. 179), „either heterogeneous of pseudohomogeneous“ (стр. 162), „archived“ уместо achieved (стр. 51), варијанте речи readsorption, удвојени члан и реченица прекинута усред мисли у закључку; словне грешке у именима носећих референци у телу текста (у списку литературе исправно). | Лекторске природе; не мењају стручни смисао ни резултате. | 51, 162, 179 | The ACS Style Guide, 3. изд., 2006; Правопис српскога језика 2010 |
| 13 | Реферат: број поглавља, број извора и звање члана комисије | Реферат рачуна седам поглавља, укључујући литературу као засебно, док рад има шест нумерисаних поглавља и одвојен списак литературе; наводи 387 цитираних извора, док је стваран број различитих извора мањи због дупликата; звање једног члана комисије наводи другачије него насловне стране рада; тврдњу да је реч о првом случају употребе детаљне кинетике за оптимизацију оваквог реакторског система сам текст рада не може потврдити. | Формалне несагласности које се не тичу садржаја; тврдњу о првенству ваља читати као образложен утисак, а не као из рада доказану чињеницу. | реферат | аргумент из грађе самог рада |
Шта рад доноси
- Јасно постављен и проверљив истраживачки план: три циља наведена су оперативно — „Analysis of the effect of process conditions on FTS in terms of activity and product distribution.“ (стр. 15), „Development of a detailed kinetic model of FTS.“ (стр. 15) и „Application of detailed kinetics in fixed-bed reactor modeling and optimization.“ (стр. 16) — и дословно се пресликавају у поглавља 3, 4 и 5.
- Оперативна таксономија кинетичких модела која усмерава сав каснији избор: „Kinetic models of FTS can be classified into three major categories“ (стр. 35), уз критичко читање литературе о реадсорпционим моделима.
- Носећи изворни допринос кинетичког дела: концепт десорпције 1-олефина зависне од броја угљеникових атома, изведен из физичких основа, с математички показаном еквиваленцијом према реадсорпционом приступу — „readsorption model can only give the correct prediction of this behavior if“ (стр. 159) — и с дискриминацијом модела у којој су боље-слажући облици одбачени због кршења ограничења ентропије адсорпције.
- Самостална анализа вероватноће раста ланца на кобалтним подацима, која носиоцем промене селективности проглашава кинетику настанка метана, а не реадсорпцију: „which could not be caused by readsorption, as ethylene is the lowest olefin“ (стр. 116).
- Изразита доследност у навођењу сопствених ограничења: „Conducting experiments and collecting data for Fischer-Tropsch synthesis in a laboratory reactor setup was not a part of this thesis.“ (стр. 65); „Using only the available data it is not possible to fully validate the proposed“ (стр. 169); научни апарат уредан и потпун, с усклађеним позивањима и списком литературе у оба смера.
Образложење примедби
- Оцена и избор модела на једном скупу података, без провере на издвојеним подацимаДомет исказа аутор је примерено ограничио, а тражење независне валидационе базе делом је захтев за обимнијим истраживањем.
- Кинетика оцењена на ситним честицама примењена на цевни реактор с честицама 1–3 mmАутор ограничење изричито означава и упућује на наставак рада; спрезање детаљне кинетике с реакторским моделом је методолошки помак који прости модели не омогућавају.
- Оптимизовани услови наведени на три места с три неусаглашена скупа бројеваМеродаван је усаглашени резултат из текста и Табеле 5.4; сама оптимизација и физичко образложење оптимума остају ваљани.
- Снага статистичке дискриминације модела прецењена, без интервала поверења по параметруВишеодзивни третман је легитиман оквир, а F-вредности вишеструко премашују границе, па избор модела опстаје и уз конзервативнији број степени слободе; изостанак тежинских фактора поштено је образложен.
- Пад притиска виши од очекиваног усвојен као активно ограничење оптимизацијеАутор физички исправно образлаже положај оптимума и сам наводи да пад притиска није значајно утицао на продуктивност и селективност.
- Провера реакторског модела само поређењем реда величинеРад ограничење поштено квалификује и признаје да пуна валидација с расположивим подацима није могућа (стр. 169).
- Уредничке омашке поглавља 5: ознаке једначина, упуте и опис нумеричке методеСлике и табеле нумерисане су потпуно и потписане у складу с обавезујућим упутством Универзитета; како симулације дају пораст температуре, омашка у знаку је нотацијска, не рачунска.
- Погрешна референца уз опис гвозденог огледа и омашке у називу извора податакаОпис експерименталног дизајна упућен је на спољни извор, али рад унапред и поштено означава поделу посла (туђи оглед, сопствена обрада), узима у обзир деактивацију катализатора и приказује поновљивост.
- Преоштре и прекомерно уопштене формулације у поглављу 3Прекомерно уопштавање делом је самоисправљено у телу текста; сама анализа α1/αC2+ и аргумент против доминантне улоге реадсорпције остају самостални доприноси.
- Омашке у приказу преноса масе и терминологијиПреглед литературе је компетентан, уравнотежен и добро документован, а приписивања механизама и концепата вођена су непосредно из извора.
- Научни апарат: дуплиране јединице, непотпуне референце, изостанак кода и интервала поверењаУсклађеност позивања и списка литературе потпуна је у оба смера, примарни извори наведени су из оригинала, сва три интернет извора носе датум приступа, а сирови подаци, оцењени параметри и једначине дати су у тексту, па је реконструкција модела у начелу могућа.
- Језичке омашке у енглеском текстуУзете појединачно, ове омашке су испод прага који дира садржај, па се воде као посао лектуре, не као мана рада.
- Реферат: број поглавља, број извора и звање члана комисијеРеферат у највећем делу верно приказује рад: обим (219 страна, 12 табела, 66 слика) слаже се с текстом, изричито наводи да извођење огледа није било део тезе и тачно описује оцену параметара; оцену о научном доприносу поткрепљују објављени радови аутора.
Осврт не доводи у питање оригиналност рада и не тврди повреду академске честитости. Реч је о стручној критици тачности, метода и закључивања, с дословним наводима и странама дисертације. Као мерило су коришћени само литература, стандарди и прописи који су постојали пре предаје рукописа комисији (фебруар 2015).