Kristallanish jarayoni laboratoriyadan ishlab chiqarish miqyosiga koʻchirilganda aralashtirish, harorat va turbulentlik endi tankning har bir nuqtasida bir xil boʻlib qolmaydi. Devor yaqinidagi sovuqroq hududlar, aralashtirgich atrofidagi yuqori turbulentlik va oqimi past turgʻun hududlar turli yadrolanish va kristall oʻsish sharoitlarini yaratishi mumkin. Natijada ayni kimyoviy modda tankning turli hududlarida turli oʻlcham va shakllarda kristallanishi mumkin.

Tadqiqotchilar bu muammoni hisoblash suyuqliklar dinamikasidan olingan maʼlumotlarni ikki oʻlchamli populyatsiya balansi modeli bilan birlashtiruvchi fizikadan ilhomlangan koʻp boʻlmali doira orqali hal qildilar. Oʻn litrli, ikki aralashtirgichli kristallizator yettita gidrodinamik hududga ajratildi. Har bir hudud uchun harorat, oqim almashinuvi, turbulent energiya dissipatsiyasi, oʻta toʻyinganlik, yadrolanish va kristallarning uzunlik-kenglik yoʻnalishlaridagi oʻsishi alohida hisoblandi.

Model standart momentlar usulining hisoblash tezligini yuqori aniqlikdagi chekli hajm usulining batafsil kristall oʻlchami taqsimoti bilan birlashtiradi. Bir boʻlmali 400 × 400 aniqlikdagi taqqoslashda hisoblash vaqti 262 soniyadan 36,2 soniyaga tushdi. Yettita boʻlmali 100 × 100 aniqlikdagi taqqoslashda toʻliq chekli hajm yechimi 1.245,6 soniya davom etgan boʻlsa, gibrid model ayni amalni 53,1 soniyada tugatdi; bu taxminan 23,5 baravar tezlashishni anglatadi.

Natijalar eritmaning oʻrtacha konsentratsiyasi tank boʻylab deyarli bir xil qolgan taqdirda ham devordan kelib chiqadigan harorat farqlari mahalliy oʻta toʻyinganlik va oʻsish tezligini oʻzgartirishi mumkinligini koʻrsatdi. Turbulent energiya dissipatsiyasi yadrolanish tenglamasiga kiritilganda aralashtirgichga yaqin hududlar ustun yadrolanish markazlariga aylandi. Bu farqlar yakuniy kristall oʻlchami taqsimoti va tomonlar nisbatini sezilarli darajada oʻzgartirdi.

Nega kristallanishni masshtablash qiyin?

Kristallanish dori vositalari, agrokimyoviy va maxsus kimyoviy moddalarni ajratish hamda tozalashda qoʻllanadigan asosiy jarayonlardan biridir. Mahsulot faqat kimyoviy tozaligi bilan baholanmaydi. Kristallarning kattaligi, oʻlcham taqsimoti, shakli, tomonlar nisbati va kristall shakli keyingi filtrlash, yuvish, quritish, tashish va formulatsiya bosqichlariga ham taʼsir qiladi.

Ignasimon yoki plastinasimon anizotrop kristallar ayniqsa muhimdir. Uzun va ingichka kristallar filtr qatlamlarida yuqori qarshilik yaratishi, oquvchanlik muammolariga olib kelishi yoki mexanik ishlov paytida sinishi mumkin. Kristall faqat bitta xarakterli diametr bilan ifodalanadigan bir oʻlchamli modellar uzunlik va kenglik turli tezliklarda oʻsadigan bunday tizimlarni yetarlicha tushuntira olmaydi.

Laboratoriya miqyosida aralashma ideal holatga yaqinlashishi mumkin. Tank kattalashgan sari esa:

• Devor va markaz oʻrtasida harorat farqlari yuzaga kelishi mumkin.
• Aralashtirgich atrofida turbulent energiya dissipatsiyasi ortib, uzoq hududlarda kamayishi mumkin.
• Aralashtirish vaqti bilan kristall oʻsish vaqti bir-biriga yaqinlashishi mumkin.
• Yadrolanish faqat oʻrtacha oʻta toʻyinganlikka emas, mahalliy toʻqnashuv va yemirilish sharoitlariga bogʻliq boʻlib qolishi mumkin.
• Tank boʻyicha oʻrtacha qiymatlar mahsulot hosil boʻladigan mahalliy sharoitlarni yashirishi mumkin.

Eng batafsil yechim hisoblash suyuqliklar dinamikasi va populyatsiya balansi tenglamalarini tankning har bir katagida bir vaqtning oʻzida yechishdir. Biroq bir necha soatlik kristallanish jarayonida ikki ichki oʻlchamga ega kristall taqsimotini millionlab oqim kataklari bilan birga hisoblash hatto bitta ishlash sharoiti uchun ham soatlar yoki kunlar olishi mumkin.

Tadqiqotning asosiy yondashuvi

Taklif etilgan PIMC-PBM yondashuvi toʻliq aralashgan bir hududli model bilan toʻliq CFD-PBE yechimi orasida oraliq darajani yaratadi. CFD faqat reaktorning tez gidrodinamik va qisqa muddatli issiqlik xatti-harakatini aniqlash uchun ishlatiladi. Keyin bu maʼlumot uzoq kristallanish jarayonini ancha past xarajat bilan hisoblaydigan boʻlmali modelga uzatiladi.

Doira toʻrtta asosiy komponentdan iborat:

1. Reaktorning tezlik va turbulentlik maydonini hisoblaydigan CFD modeli,
2. Reaktorni jismoniy jihatdan mazmunli hududlarga ajratadigan optimallashtirish usuli,
3. Har bir hududdagi haroratni hisoblaydigan past tartibli issiqlik surrogat modeli,
4. Kristall uzunligi va kengligini birgalikda kuzatadigan ikki oʻlchamli populyatsiya balansi modeli.

Oʻn litrli reaktor yettita hududga ajratildi

Tadqiqotda ikki parrakli tuzilishga ega, 10 litr ish hajmli aralashtiriladigan davriy kristallizator modellashtirildi. Taxminan 1,5 million katakdan iborat O-grid toʻri ishlatildi. 26-sahifadagi reaktor geometriyasi va toʻr tasviri pastki va yuqori aralashtirgich parraklari, suyuqlik-havo interfeysi hamda devor yaqinidagi zichlashtirilgan hisoblash toʻrini koʻrsatadi.

CFD natijalariga koʻra reaktor quyidagi yettita hududga ajratildi:

Boʻlma chegaralari qoʻlda yoki faqat vizual baholash bilan belgilanmagan. Tadqiqotchilar NSGA-II koʻp maqsadli optimallashtirish algoritmidan foydalanib, har bir boʻlma ichidagi tezlik taqsimoti va turbulent energiya dissipatsiyasini imkon qadar bir jinsli qilishga harakat qildilar.

Ikki maqsad funksiyasi tezlik va energiya dissipatsiyasining kvartillar oraligʻiga asoslangan edi:

\[ f_1(\mathbf{d})=\operatorname{IQR}(u) \]

\[ f_2(\mathbf{d})=\operatorname{IQR}(\varepsilon) \]

Bu yerda u tezlik kattaligini, ε turbulent kinetik energiya dissipatsiyasi tezligini, d esa boʻlma geometriyasini belgilaydigan oʻlchamlarni koʻrsatadi. 27-sahifadagi yettita Pareto grafigida har bir boʻlma uchun tezlik bir jinsliligi va energiya dissipatsiyasi bir jinsliligi oʻrtasidagi muvozanat hamda tanlangan yechim qizil yulduzlar bilan koʻrsatilgan.

Eng yuqori sezgirlik pastki aralashtirgich joylashgan toʻrtinchi boʻlmada kuzatildi. Boʻlma chegaralari aralashtirgich hosil qilgan qayta aylanish tuzilmalarini yaxshiroq qamrab oladigan tarzda oʻzgartirilganda hudud ichidagi energiya dissipatsiyasi taqsimoti sezilarli toraydi.

CFD modeli qanday qurildi?

Oqim hisoblari OpenFOAM muhitida bajarildi. Suyuqlik siqilmaydigan va Nyuton suyuqligi deb qabul qilindi. Aralashtirgich aylanishi parraklarning vaqtga bogʻliq harakatini toʻgʻridan-toʻgʻri yechish oʻrniga Koʻp tayanchli sanoq tizimi usuli bilan ifodalandi.

Aylanuvchi sanoq tizimidagi impuls tenglamasi tadqiqotda quyidagicha berilgan:

\[ \nabla\cdot(\rho\mathbf{u}_r\mathbf{u}) +\rho(\boldsymbol{\omega}\times\mathbf{u}) = -\nabla p+\rho\mathbf{g}+\mu\nabla^2\mathbf{u} \]

Bu yerda u mutlaq tezlikni, u_r aylanuvchi sanoq tizimiga nisbatan nisbiy tezlikni, ω burchak tezligini, p bosimni va μ dinamik qovushqoqlikni ifodalaydi.

Turbulentlik standart k-ε RANS modeli bilan hisoblandi. Devorlarda sirpanmaslik sharti va devor funksiyalari, suyuqlik-havo interfeysida esa sirpanish sharti ishlatildi. Barcha oʻzgaruvchilar uchun ikkinchi tartibli fazoviy diskretlashtirish qoʻllanib, sonli diffuziyani kamaytirish maqsad qilindi.

Oqim maydoni hisoblari 64 protsessor yadrosida bajarildi va har bir CFD sharoiti taxminan 20 daqiqa davom etdi. Haroratning vaqtga bogʻliq CFD namunalari esa taxminan 24 soat devor vaqtini talab qildi.

Uzoq muddatli issiqlik uzatish qanday tezlashtirildi?

Toʻrt soatlik yoki undan uzun kristallanish jarayonining barchasini CFD bilan yechish oʻrniga uch daqiqalik qisqa issiqlik CFD namunalari ishlatildi. Devor harorati 273 K va 373 K boʻlgan sovitish va isitish sharoitlari namuna olinib, past tartibli issiqlik modeli yaratildi.

Har bir boʻlmaning energiya balansi quyidagi tenglama bilan ifodalandi:

\[ \rho C_pV_i\frac{dT_i}{dt} = \sum_{j\neq i}F_{j\rightarrow i}C_p(T_j-T_i) + hA_i(T_w-T_i) \]

V_i boʻlma hajmini, F_j→i boʻlmalar orasidagi hajmiy oqimni, T_i boʻlma haroratini, T_w devor haroratini va h samarali issiqlik uzatish koeffitsiyentini koʻrsatadi.

Yettita gidrodinamik boʻlma issiqlik hisobini yanada kamaytirish uchun ikkita asosiy guruhga birlashtirildi:

• Devorga qoʻshni 1-3-raqamli boʻlmalar,
• Markaz va aralashtirgich taʼsiridagi 4-7-raqamli boʻlmalar.

Samarali issiqlik uzatish koeffitsiyenti haroratning birinchi tartibli funksiyasi sifatida modellashtirildi:

\[ h=aT_i+b \]

Parametrlar avval genetik algoritm yordamida global qidirildi, soʻng chiziqli boʻlmagan eng kichik kvadratlar usuli bilan takomillashtirildi. Ishonch oraliqlari Levenberg-Marquardt yondashuvi bilan hisoblandi.

Kristallar ikki oʻlchamda kuzatildi

Model kristallarni bitta ekvivalent diametr bilan emas, uzunlik l va kenglik w koordinatalari bilan ifodalaydi. Har bir boʻlmadagi son zichligi n_z(l,w,t) bilan koʻrsatiladi.

Soddalashtirilgan ikki oʻlchamli populyatsiya balansi tenglamasi quyidagi tuzilishga ega:

\[ \frac{\partial n_z}{\partial t} + \frac{\partial(n_zG_{l,z})}{\partial l} + \frac{\partial(n_zG_{w,z})}{\partial w} = \delta(l-l_n)\delta(w-w_n)B_z + \sum_{k\neq z} \frac{F_{k\rightarrow z}}{V_z} (n_k-n_z) \]

Tenglamaning birinchi hadi vaqtga bogʻliq oʻzgarishni, ikkinchi va uchinchi hadlari kristallarning uzunlik va kenglik yoʻnalishidagi oʻsishini, Dirac delta hadi yangi yadrolarning muayyan boshlangʻich oʻlchamda hosil boʻlishini, soʻnggi had esa kristallarning boʻlmalar oʻrtasida tashilishini ifodalaydi.

Boshlangʻich urugʻ ikki oʻlchamli Gauss taqsimoti bilan belgilandi. Kristallarning oʻrtacha uzunligi 200 µm, oʻrtacha kengligi 100 µm; ikki yoʻnalishdagi taqsimot parametrlari esa 150 µm qilib tanlandi.

Model sinish, aglomeratsiya va erishni eʼtiborsiz qoldiradi. Shu sababli natijalar oʻsish va ikkilamchi yadrolanish ustun boʻlgan kristallanish sharoitlari bilan cheklangan.

Momentlar usuli hisoblashni qanday kamaytirdi?

Har bir boʻlmada barcha uzunlik-kenglik kataklarini vaqt davomida yechish oʻrniga avval son taqsimotining aralash momentlari hisoblandi:

\[ \mu_{ij,z}(t) = \int_0^\infty\int_0^\infty w^il^j n_z(l,w,t)\,dl\,dw \]

Past tartibli momentlar kristallar soni, oʻrtacha oʻlcham va umumiy kristall hajmini oz sonli differensial tenglama bilan ifodalaydi. Shunday qilib asosiy kristallanish dinamikasini standart momentlar usuli bilan tez yechish mumkin.

Biroq faqat momentlar ishlatilganda butun taqsimot, koʻp choʻqqili tuzilma yoki dum xatti-harakati yoʻqolishi mumkin. Shu sababli tadqiqotchilar jarayon dinamikasini SMOM bilan yechib, yakuniy ikki oʻlchamli kristall oʻlchami taqsimotini yuqori aniqlikdagi chekli hajm usuli bilan qayta qurdilar.

Yuqori aniqlikdagi chekli hajm usuli

HRFVM usulida uzunlik va kenglik sohasi muntazam hisoblash toʻriga ajratildi. Kristallar oʻsgan sari son zichligi kichik oʻlchamli kataklardan kattaroq oʻlchamli kataklarga koʻchdi.

Birinchi tartibli yuqori oqim sxemasi barqarorlikni taʼminlar ekan, ikkinchi tartibli tuzatish sonli diffuziyani kamaytirdi. Keskin taqsimot frontlarida fizik boʻlmagan tebranishlarning oldini olish uchun oqim cheklagichi ishlatildi:

\[ \phi(\theta)=\frac{|\theta|+\theta}{\theta+1} \]

Silliq hududlarda usul ikkinchi tartibli aniqlikka yaqinlashadi, keskin gradientlarda esa barqarorroq birinchi tartibli xatti-harakatga qaytadi.

Oʻta toʻyinganlik, oʻsish va yadrolanish

Har bir boʻlmaning oʻta toʻyinganligi mahalliy konsentratsiya va haroratga bogʻliq eruvchanlikdan hisoblandi:

\[ \sigma_z=\frac{C_z-C_s(T_z)}{C_s(T_z)} \]

Eruvchanlik ikkinchi tartibli harorat polinomi bilan ifodalandi:

\[ C_s=a+bT_z+cT_z^2 \]

Kristallning uzunlik va kenglik yoʻnalishlaridagi oʻsish tezliklari alohida munosabatlar bilan berildi:

\[ G_{l,z}=k_{gl}\sigma_z^{g_l} \]

\[ G_{w,z}=k_{gw}\sigma_z^{g_w} \]

Bu ajratish kristallning ikki yoʻnalishda bir xil tezlikda oʻsmasligi va tomonlar nisbatining jarayon davomida oʻzgarishiga imkon beradi.

Ikkilamchi yadrolanish esa oʻta toʻyinganlik, mavjud kristall hajmi va mahalliy energiya dissipatsiyasi bilan bogʻlandi:

\[ B_z=k_s\sigma_z^sM_{t,z}^{\,j}D_{\varepsilon,z} \]

Energiya dissipatsiyasi taʼsiri yoqilganda aralashtirgich atrofidagi yuqori turbulent hududlarda kristall-kristall toʻqnashuvlari va sirt yemirilishi bilan bogʻliq yadrolanish kuchayishi mumkin.

Hisoblash samaradorligi

Birinchi taqqoslash bir boʻlmali va toʻliq aralashgan kristallanish modeli ustida bajarildi. Toʻliq HRFVM bilan gibrid SMOM-HRFVM ayni toʻr aniqliklarida taqqoslandi.

Xabar qilingan moment farqlari asosan 0,514 foizdan past boʻlib qoldi. Shunday qilib tezlik yutugʻiga asosiy momentlar aniqligida katta yoʻqotishsiz erishildi.

Yetti boʻlmali hisoblashda farq yanada yaqqol boʻldi:

Tadqiqotchilar umumiy tezlik yutugʻini muammo oʻlchamiga qarab bir va ikki tartib oraligʻida taʼriflaydilar. Jadvaldagi bevosita taqqoslanadigan eng kuchli natija taxminan 23,5 baravar tezlashishdir.

Konsentratsiya teng boʻlsa ham harorat farqi oʻsishni oʻzgartirdi

Fazoviy bir jinsli boʻlmaslik modelga bosqichma-bosqich qoʻshildi. Birinchi holatda boʻlmali tuzilma ideal aralashish xatti-harakatiga yaqin edi. Ikkinchi holatda devordan kelib chiqadigan harorat farqlari, keyingi holatlarda esa energiya dissipatsiyasiga bogʻliq yadrolanish taʼsirlari qoʻshildi.

30-sahifadagi grafiklarda barcha boʻlmalarning erigan modda konsentratsiyasi egri chiziqlari deyarli ustma-ust tushadi. Buning sababi oʻrganilgan asosiy sharoitlarda boʻlmalararo aralashish kristallanish kinetikasidan tezroq boʻlishidir.

Bunga qarshi ravishda ayni grafikdagi oʻta toʻyinganlik va kenglik yoʻnalishidagi oʻsish tezligi egri chiziqlari ikki guruhga ajraladi. Devorga yaqin 1-3-raqamli boʻlmalar pastroq harorat sababli yuqoriroq oʻta toʻyinganlik va oʻsish tezligiga erishadi. Markazdagi 4-7-raqamli boʻlmalar esa issiqroq qoladi.

Bu natija muhim masshtablash mexanizmini koʻrsatadi: butun tankdan olingan oʻrtacha konsentratsiya oʻlchovi bir jinsli koʻrinsa ham, mahalliy harorat farqlari kristall oʻsishini bir jinsli boʻlmagan holga keltirishi mumkin.

Turbulentlik yadrolanish joyini belgiladi

Energiya dissipatsiyasi yadrolanishga taʼsir qilmagan dastlabki holatlarda boʻlmalarning yadrolanish egri chiziqlari bir-biriga yaqin edi. Energiya dissipatsiyasiga bogʻliqlik kuchaygani sari egri chiziqlar yaqqol ajraldi.

31-sahifadagi yadrolanish grafiklari aralashtirgichga yaqin va energiya dissipatsiyasi yuqori boʻlgan boʻlmalar ustun yadrolanish markazlariga aylanganini koʻrsatadi. Turbulentligi pastroq hududlarda esa yadrolanish katta darajada bostirilgan.

Bu mexanizm yakuniy mahsulot taqsimotiga qarama-qarshi ikki yoʻl bilan taʼsir qiladi:

• Yuqori turbulent hududlarda koʻproq yangi yadro hosil boʻlib, kichik kristallar populyatsiyasi ortadi.
• Yadrolanishi past hududlarda mavjud kristallar oʻsishda davom etib, kattaroq zarrachalarni hosil qiladi.

Kristall oʻlchami va tomonlar nisbati qanday oʻzgardi?

Nominal kinetika va 0,0111 °C/s sovitish tezligida haroratning bir jinsli boʻlmasligi oʻz-oʻzicha yakuniy tomonlar nisbatini juda oz oʻzgartirdi. Energiya dissipatsiyasiga bogʻliq yadrolanish kuchaygani sari oʻrtacha kristall oʻlchamlari sezilarli kattalashdi.

Birinchi va beshinchi holatlar taqqoslanganda hajm boʻyicha vaznlangan kenglik taxminan 24,5 foiz, uzunlik esa taxminan 25,5 foiz oshdi. Tomonlar nisbatidagi oʻzgarish cheklanganroq boʻlib, 2,23 dan 2,25 gacha koʻtarildi.

32-sahifadagi ikki oʻlchamli taqsimotlar oʻrtacha tomonlar nisbati oʻxshash koʻrinsa ham, zarracha zichligining uzunlik-kenglik tekisligidagi oʻrni va taqsimot shakli oʻzgarganini koʻrsatadi. Shu sababli faqat oʻrtacha tomonlar nisbatiga qarash mahsulot taqsimotidagi barcha farqlarni ochib bermaydi.

Sekin sovitish uzunroq kristallarni hosil qildi

Sovitish tezligi 0,0111 °C/s dan 0,0011 °C/s gacha tushirilganda kristallarda oʻsish uchun koʻproq vaqt boʻldi. Uzunlik yoʻnalishidagi oʻsish kenglik yoʻnalishidagidan kuchliroq boʻlgani sababli tomonlar nisbati taxminan 2,91-2,95 oraligʻiga chiqdi.

Eng kuchli gidrodinamik bir jinsli boʻlmaslikka ega beshinchi holatda sekin sovitishda hajm boʻyicha vaznlangan kenglik 269,7 µm, uzunlik 785,1 µm va tomonlar nisbati 2,91 deb hisoblandi.

Bu topilma sovitish dasturi faqat oʻrtacha kristall oʻlchamini emas, oʻsishning ikki oʻq orasidagi taqsimotini ham oʻzgartirishi mumkinligini koʻrsatadi.

Aralashish va kristallanish vaqtlari raqobatlashganda nima boʻldi?

Tadqiqotchilar kristall oʻsish konstantalarini 90 baravar oshirib, kristallanish vaqtini sunʼiy ravishda qisqartirdilar. Aralashtirish vaqti oʻzgarmaganligi sababli ikki jarayon oʻrtasidagi raqobat bevosita oʻrganildi.

Asosiy kinetikada boʻlmalararo oqim erigan modda konsentratsiyasini tez tenglashtirardi. Tezlashtirilgan kinetikada esa kristallar oʻta toʻyinganlikni aralashish tenglashtira olishidan tezroq sarflay boshladi. Natijada bitta konsentratsiya egri chizigʻi oʻrniga yettita alohida boʻlma egri chizigʻi hosil boʻldi.

Devorga yaqin sovuqroq hududlarda oʻsish va erigan modda sarfi tezlashdi. Markaziy hududlar esa boʻlmalararo almashinuv bilan suyuldi. Shunday qilib masshtablashda hal qiluvchi omil faqat aralashtirish tezligi yoki faqat kristall oʻsish tezligi emasligi; balki bu ikki vaqtning bir-biriga nisbati ekani koʻrsatildi.

Tezlashtirilgan oʻsishda yadrolanish tezliklari boʻlmalar oʻrtasida katta farqlar koʻrsatgan boʻlsa-da, jarayon oʻsish ustunligida qoldi. Shu sababli hajm boʻyicha vaznlangan yakuniy kristall taqsimotlari bir-biriga yaqin qoldi. 35-sahifadagi taqsimotlarda tomonlar nisbati asosiy holatda 1,44, bir jinsli boʻlmagan holatlarda taxminan 1,50 deb koʻrinadi.

Bu muhandislik nuqtayi nazaridan nimani anglatadi?

Tadqiqot kristallanishni masshtablashda faqat tankning oʻrtacha harorati, aralashtirgich aylanishi yoki oʻrtacha oʻta toʻyinganlikka qarash yetarli boʻlmasligi mumkinligini koʻrsatadi. Ayni oʻrtacha konsentratsiya ortida har xil harorat, oʻsish va yadrolanish tarixlari boʻlishi mumkin.

Taklif etilgan doira ayniqsa quyidagi maqsadlarda baholanishi mumkin:

• Laboratoriya va ishlab chiqarish reaktorlari oʻrtasidagi morfologiya oʻzgarishini bashorat qilish,
• Aralashtirgich holati va tezligining kristall oʻlchamiga taʼsirini oʻrganish,
• Devorni sovitish dasturlari va issiqlik sikllarini loyihalash,
• Ignasimon yoki plastinasimon kristallarning tomonlar nisbatini nazorat qilish,
• Koʻplab ishlash sharoitlarini raqamli muhitda tezkor tekshirish,
• Hisoblash xarajati past kristallanish raqamli egizaklarini ishlab chiqish.

Shu bilan birga model tajribaviy tekshiruvsiz bevosita ishlab chiqarish boshqaruviga ulanmasligi kerak. Har bir kimyoviy birikma uchun eruvchanlik, oʻsish, yadrolanish va energiya dissipatsiyasiga bogʻliqliklar tajriba yoʻli bilan aniqlanishi lozim.

Tadqiqot usuli va topilmalari

Tadqiqot dizayni

Tadqiqot CFD orqali axborotlangan koʻp boʻlmali ikki oʻlchamli populyatsiya balansi modelini ishlab chiqish va sonli taqqoslashlar bilan baholashga asoslangan hisoblash kimyo muhandisligi tadqiqotidir.

Asosiy ishlash sharoitlari

Dasturiy taʼminot va hisoblash resurslari

• HRFVM kodi C tilida ishlab chiqilib, MATLABga MEX sifatida ulandi.
• SMOM tenglamalari MATLAB skriptlari bilan yechildi.
• PBM taqqoslashlari Intel Core i9-12900KF, 32 GB DDR5 xotirali ish stoli kompyuterida bitta oqimda bajarildi.
• CFD oqim hisoblari 64 CPU yadrosida taxminan 20 daqiqa davom etdi.
• Vaqtinchalik harorat CFD hisoblari taxminan 24 soat davom etdi.
• Boʻlmalarni optimallashtirish bitta CPU yadrosida taxminan 30 daqiqa davom etdi.
• Boʻlmalararo oqimni ajratib olish 8 CPU yadrosida taxminan 5 daqiqada tugallandi.

Asosiy sonli topilmalar

• Gibrid usul bir boʻlmali modellarda aniqlikka qarab taxminan 3,2-7,2 baravar tezlashishni taʼminladi.
• Yetti boʻlmali 100 × 100 toʻrda hisoblash vaqti 1.245,6 soniyadan 53,1 soniyaga tushdi.
• 200 × 200 va 300 × 300 koʻp boʻlmali toʻliq HRFVM yechimlari maqbul vaqt ichida tugallanmagan boʻlsa, gibrid usul mos ravishda 841,5 va 1.667,3 soniyada ishladi.
• Asosiy kinetik sharoitlarda boʻlmalararo konsentratsiya farqlari eʼtiborsiz darajada qoldi.
• Devor harorati farqlari konsentratsiya bir xil qolganida oʻta toʻyinganlik va oʻsish tezligini hududlar boʻyicha oʻzgartirdi.
• Energiya dissipatsiyasiga bogʻliq yadrolanish aralashtirgich atrofidagi hududlarni ustun yadrolanish markazlariga aylantirdi.
• Nominal tez sovitishda tomonlar nisbati 2,23-2,25 oraligʻida qoldi.
• Sekin sovitishda tomonlar nisbati taxminan 2,91-2,95 oraligʻiga koʻtarildi.
• Oʻsish konstantalari 90 baravar oshirilganda konsentratsiya maydoni yettita alohida hududiy egri chiziqqa ajraldi.
• Tezlashtirilgan oʻsishda hududiy yadrolanishdagi kuchli farqlarga qaramay, yakuniy hajm taqsimoti oʻsish ustunligi sababli cheklangan darajada oʻzgardi.

Tekshirish va taqqoslash chegarasi

Modelning sonli aniqligi SMOM momentlari bilan HRFVM natijalarini taqqoslash hamda turli toʻr aniqliklaridagi hisoblash vaqtlari orqali baholandi. Reaktorning CFD toʻri toʻr yaqinlashuvi tadqiqotidan keyin taxminan 1,5 million katak qilib tanlangani bildirilgan boʻlsa-da, toʻr yaqinlashuvi qiymatlarining batafsil jadvali taqdim etilmagan.

Tadqiqotda jismoniy reaktordan olingan tezlik, harorat, kristall oʻlchami taqsimoti yoki tomonlar nisbati oʻlchovlari bilan bevosita tekshiruv oʻtkazilmagan. Shu sababli natijalar model ichidagi sonli taqqoslash va mexanistik izchillik darajasida baholanishi kerak.

Muhim cheklovlar

• Tadqiqot taqrizdan oʻtmagan preprintdir.
• Yakuniy kristall oʻlchami va tomonlar nisbati natijalari jismoniy tajriba bilan tasdiqlanmagan.
• CFD tezlik va turbulentlik maydonlarining tajribaviy oʻlchovlar bilan taqqoslovi taqdim etilmagan.
• Sinish, aglomeratsiya va erish mexanizmlari modelga kiritilmagan.
• Har bir boʻlma oʻz ichida harorat, konsentratsiya va energiya dissipatsiyasi jihatidan bir jinsli deb qabul qilingan.
• Oqim maydoni kristallanish davomida kvazibarqaror deb qabul qilingan.
• Harorat oʻzgarishlari zichlik farqi va tabiiy konveksiya orqali oqim maydonini oʻzgartirmaydi deb faraz qilingan.
• Issiqlik modeli yettita gidrodinamik boʻlmani ikkita harorat hududiga qisqartiradi.
• Yadrolanishning energiya dissipatsiyasiga bogʻliqligi empirik had bilan ifodalangan.
• Kinetik parametrlarni boshqa moddalar yoki reaktorlarga bevosita koʻchirib boʻlmaydi.
• Tadqiqotda ochiq maʼlumotlar toʻplami yoki manba kod repozitori koʻrsatilmagan.
• 6-jadval va 7-jadvalda xato ustuni sarlavhalari ayni moment belgisi bilan takrorlangan; qaysi turli momentlar taqqoslangani jadval sarlavhalaridan aniq anglashilmaydi.
• 33-sahifadagi 11-shakl izohi tasvirlangan konsentratsiya, oʻta toʻyinganlik va oʻsish grafiklariga mos kelmaydi; bu preprintda shakl izohi yoki raqamlash xatosi borligini koʻrsatadi.

Manba va usul qaydi

Asl tadqiqot: “A Physics-Inspired Multi-compartmental (PIMC) Population Balance Framework for Size and Aspect Ratio Control in Crystallization Scale-Up.”

Mualliflar: Yung-Shun Kang, Antonello Raponi va Zoltan K. Nagy.

Masʼul muallif: Zoltan K. Nagy. Teng hissa bayonoti mavjud emas.

Muassasalar: Davidson School of Chemical Engineering, Purdue University, West Lafayette, Indiana, AQSH; Mines Saint-Étienne, Univ Lyon, CNRS, UMR 5307 LGF, Centre SPIN, Saint-Étienne, Fransiya.

Muallif hissalari: Yung-Shun Kang konsepsiyalash, dasturiy taʼminot, tadqiqot, usul va dastlabki qoralamani; Antonello Raponi dasturiy taʼminot, usul va matnni baholashni; Zoltan K. Nagy konsepsiyalash, maslahat, loyiha boshqaruvi va moliyalashtirishni olishni oʻz zimmasiga olgan.

Manba turi: Asl hisoblash kimyo muhandisligi tadqiqoti va preprint.

Yil: 2026.

Taqriz holati: Tadqiqot taqrizdan oʻtmagan. Topilmalar taqrizli va yakuniy jurnal natijalari sifatida baholanmasligi kerak.

Platforma: SSRN.

Taqrizli jurnal: Hech qanday taqrizli jurnalda qabul qilingani yoki nashr qilingani haqida maʼlumot yoʻq.

DOI:10.2139/ssrn.6948349

Tadqiqotning rasmiy havolasi:SSRN tadqiqot qaydi

Moliyalashtirish: Tadqiqot AQSH Milliy ilmiy jamgʻarmasining 2132142-raqamli granti bilan qisman qoʻllab-quvvatlangan. Corteva Agrisciencening qisman moliyaviy koʻmagi va tadqiqot hamkorligi ham koʻrsatilgan.

Manfaatlar toʻqnashuvi: Matnda manfaatlar toʻqnashuvi bayonoti mavjud emas.

Maʼlumotlar va kodga kirish: Ochiq maʼlumotlar toʻplami, model fayllari yoki manba kod repozitori havolasi berilmagan.

Ushbu mazmundagi usul, tenglamalar, CFD sharoitlari, hisoblash vaqtlari, kristall oʻlchami natijalari, shakl va jadval sharhlari, topilmalar va cheklovlar koʻrib chiqilgan tadqiqotga asoslanadi. Tashqi manbalardan faqat nom, muallif shaxsi, muassasa, DOI, platforma va tadqiqotning rasmiy havolasini bibliografik tekshirish maqsadida foydalanilgan; tadqiqotda boʻlmagan ilmiy natija qoʻshilmagan.