Gipertovushli apparatlar atrofidagi oqimni hisoblaydigan sonli usullar kuchli zarba to‘lqinlarini ushlashi, zichlik va bosimni musbat saqlashi, barqaror ishlashi va fizik saqlanish qonunlarini imkon qadar aniq ifodalashi kerak. Biroq keng qo‘llanadigan Rusanov, HLL, HLLC va Roe kabi Godunov turidagi yechuvchilar barqaror Euler oqimlarida oqim chiziqlari bo‘ylab o‘zgarmas qolishi kerak bo‘lgan umumiy entalpiyani sonli jihatdan saqlamasligi mumkin. Zarba orqali o‘tgan oqimda hosil bo‘ladigan sun’iy entalpi o‘zgarishi keyinchalik oqim maydonining boshqa qismlariga ko‘chib, energiya va zichlik taqsimotini buzishi mumkin.
Tadqiqotda Gallice tomonidan ishlab chiqilgan Simple Riemann yechuvchilarining energiya tarkibiy qismi o‘zgartirilib, MGallice va MGallice-2D nomli entalpiyani saqlovchi bir va ko‘p o‘lchovli sxemalar yaratildi. O‘zgartirish energiya oqimidagi sonli diffuziya massa oqimidagi diffuziyaning umumiy entalpi bilan ko‘paytmasiga teng bo‘lishini ta’minlaydigan Häenel shartiga asoslanadi.
Mach 2 tezlikdagi uchburchak bo‘rtmali kanal sinovida to‘rtta sxema ham taxminan 4.000 iteratsiyada barqaror yechimga erishdi. Klassik Gallice yechuvchilari turg‘un zarbalardan o‘tayotganda kichik, biroq aniq umumiy entalpi sakrashlarini hosil qildi. O‘zgartirilgan MGallice yechuvchilari esa kirish umumiy entalpiyasi bo‘lgan 6,3 qiymatini butun hisoblash sohasida saqladi.
Mach 20 tezlikda yarim silindr atrofidagi murakkabroq sinovda bir o‘lchovli yechuvchilar karbunkul deb ataladigan zarba beqarorligiga uchradi. Ko‘p o‘lchovli yechuvchilar bu beqarorlikni oldini oldi; biroq klassik Gallice-2D zarba ortidagi entalpi yo‘qotilishi sababli tormozlanish nuqtasidagi zichlik taqsimotini noto‘g‘ri hisoblab chiqdi. MGallice-2D umumiy entalpiyani mashina aniqligigacha o‘zgarmas saqlab, kutilgan zichlik maydonini qayta tikladi.
Tadqiqot barqaror gipertovushli Euler hisoblarida umumiy entalpiyani saqlash faqat nazariy xususiyat emasligini; kuchli zarbalar ortidagi energiya va zichlik maydonini bevosita o‘zgartirishi mumkinligini ko‘rsatadi. Shunga qaramay, taklif qilingan sxemalarning entropiya barqarorligi matematik jihatdan isbotlanmagan va usul hali yopishqoq Navier–Stokes masalalarida sinovdan o‘tkazilmagan.
Gipertovushli oqimlarni sonli yechish nega qiyin?
Gipertovushli tezlikda harakatlanayotgan apparat oldida kuchli va apparatdan ajralgan yoy zarbasi hosil bo‘ladi. Zarba bo‘ylab bosim, zichlik va harorat juda qisqa masofada keskin o‘zgaradi. Apparat ortida kengayish to‘lqinlari paydo bo‘ladi, sirtga yaqin hududlarda esa tezlik va harorat gradientlari yuzaga keladi.
Haqiqiy gipertovushli oqimlar yopishqoqlik, issiqlik o‘tkazuvchanligi, kimyoviy reaksiyalar, termokimyoviy nomuvozanat ta’sirlari, turbulentlik, sirt qizishi va eroziya kabi ko‘plab jarayonlarni o‘z ichiga olishi mumkin. Tadqiqot bu butun tizimni emas, siqiluvchan Navier–Stokes tenglamalarining zarba hosil bo‘lishiga javobgar bo‘lgan yopishqoqliksiz qismini, ya’ni Euler tenglamalarini ko‘rib chiqadi.
Sonli sxema kuchli zarbani barqaror ushlagan taqdirda ham fizik jihatdan saqlanishi kerak bo‘lgan kattalikda sun’iy o‘zgarish hosil qilishi mumkin. Tadqiqotning asosiy e’tibori barqaror Euler oqimlarida umumiy entalpiyaning sonli tarzda saqlanmasligi muammosiga qaratilgan.
Umumiy entalpi nima?
Ikki o‘lchovli Euler tenglamalarida saqlanuvchi o‘zgaruvchilar zichlik, impuls va umumiy energiyadir:
\[ \mathbf{U}= \begin{pmatrix} \rho \\ \rho\mathbf{u} \\ \rho e \end{pmatrix} \]
Bu yerda ρ zichlikni, u tezlik vektorini va e birlik massaga to‘g‘ri keladigan umumiy energiyani ifodalaydi. Umumiy energiya ichki energiya va kinetik energiya yig‘indisidir:
\[ e=\varepsilon+\frac{1}{2}\mathbf{u}\cdot\mathbf{u} \]
Umumiy entalpi esa:
\[ h=e+\frac{p}{\rho} \]
munosabati bilan aniqlanadi. p bosimdir. Tadqiqotda ideal gaz tenglamasi ishlatilgan:
\[ p=\rho\varepsilon(\gamma-1) \]
γ politropik yoki solishtirma issiqliklar nisbatidir. Sonli sinovlarda γ = 7/5 deb olingan.
Barqaror Euler oqimida entalpi nega o‘zgarmas qoladi?
Ikki o‘lchovli Euler tenglamalari massa, impuls va umumiy energiyaning saqlanishini quyidagi umumiy ko‘rinishda ifodalaydi:
\[ \frac{\partial\mathbf{U}}{\partial t} + \nabla\cdot\mathbf{F}(\mathbf{U}) =0 \]
Barqaror holatda vaqtga bog‘liq had yo‘qoladi. Massa va energiya tenglamalari birga ko‘rib chiqilganda:
\[ \mathbf{u}\cdot\nabla h=0 \]
natijasi olinadi. Bu ifoda yopishqoqlik va issiqlik o‘tkazuvchanligi bo‘lmagan barqaror Euler oqimida umumiy entalpi oqim chizig‘i bo‘ylab o‘zgarmasligini ko‘rsatadi.
Agar oqim chiziqlari kirish chegarasiga bog‘langan bo‘lsa, umumiy entalpi butun hisoblash maydonida kirish qiymati h∞ bilan bir xil qolishi kerak. Zarba bosim, harorat va zichlikni o‘zgartiradi; biroq ideal barqaror Euler yechimida umumiy entalpi zarba bo‘ylab saqlanadi.
Sonli sxemalar bu xususiyatni nega yo‘qotishi mumkin?
Chekli hajm usullarida hisoblash sohasi hujayralarga bo‘linadi. Har bir hujayradagi massa, impuls va energiya o‘zgarishi hujayra yuzalaridan o‘tuvchi sonli oqimlar bilan belgilanadi. Yuzaning ikki tomonidagi holatlar orasida Riemann masalasi tuziladi va uning taqribiy yechimidan sonli oqim olinadi.
Rusanov, HLL, HLLC va Roe kabi klassik Godunov turidagi sxemalar saqlanuvchi bo‘lsa ham, massa va energiya oqimlariga ayni fizik bog‘lanishni ta’minlaydigan tarzda sonli diffuziya qo‘shmaydi. Shu sabab turg‘un zarba bo‘ylab energiya oqimidagi sun’iy diffuziya massa oqimidagi sun’iy diffuziyaning umumiy entalpi bilan ko‘paytmasiga teng bo‘lmasligi mumkin.
Natijada zarbaning bir tomonida hisoblangan umumiy entalpi boshqa tomonidagidan farq qiladi. Bu xato keyinchalik oqim bilan tashiladi va ayniqsa zarba ortidagi past tezlikli hududlarda energiya taqsimotiga ta’sir qilishi mumkin.
Entalpi xatosining muhandislikdagi ahamiyati
Tadqiqotchi ayrim Navier–Stokes hisoblarida chegara qatlamini aniqlovchi mezonlar umumiy entalpiyaga asoslanishi mumkinligini ta’kidlaydi. Zarba sababli hosil bo‘lgan sun’iy entalpi izi devor bo‘ylab oqim quyi qismiga ko‘chirilganda, fizik chegara qatlami bilan sonli energiya xatosini ajratish qiyinlashishi mumkin.
Gipertovushli apparatlarda tormozlanish nuqtasi atrofidagi zichlik, bosim, harorat va issiqlik oqimi juda muhim. Tadqiqotning Mach 20 sinovi umumiy entalpi xatosi faqat yordamchi o‘zgaruvchini emas, zichlik maydonining shaklini ham o‘zgartirishi mumkinligini ko‘rsatadi.
Simple Riemann yechuvchisi nima?
Simple Riemann yechuvchilari boshlang‘ich chap va o‘ng holatlar orasida o‘zgarmas oraliq holatlarni belgilaydi. Ushbu oraliq holatlar turli tezliklarda harakatlanadigan uzilishlar bilan bir-biridan ajraladi:
\[ W(\xi)= \begin{cases} U_l, & \xi\leq\lambda_1 \\ U_k, & \lambda_{k-1}<\xi\leq\lambda_k \\ U_r, & \xi>\lambda_m \end{cases} \]
Bu yerda ξ = x/t o‘ziga o‘xshash o‘zgaruvchi, Ul va Ur chap va o‘ng holatlar, λk esa to‘lqin tezliklaridir.
Gallice yondashuvida yechuvchi avval Lagrange koordinatalarida quriladi. Keyin Lagrange–Euler o‘zgartirishi orqali sobit fazoviy hujayralarda ishlaydigan Euler ko‘rinishiga o‘tkaziladi. Bu usul oraliq holatlarning musbat zichlik va ichki energiya shartlarini tahlil qilishni osonlashtiradi.
Bir o‘lchovli Gallice sxemasi
Bir o‘lchovli deb ataladigan sxemada hujayra yuzasidagi sonli oqim faqat yuzaning ikki tomonidagi hujayralardan hisoblanadi. Ikki o‘lchovli to‘r ishlatilsa ham, har bir yuzadagi Riemann masalasi yuzaga normal yo‘nalishda bir o‘lchovli ko‘riladi.
Oraliq holatning normal tezligi:
\[ u_f= \frac{ \lambda_lu_{n,l} + \lambda_ru_{n,r} - \Delta p }{ \lambda_l+\lambda_r } \]
deb aniqlanadi. λl va λr Lagrange to‘lqin tezliklarini, un yuzaga normal tezlik tarkibiy qismini va Δp ikki hujayra orasidagi bosim farqini ko‘rsatadi.
Maqolada Euler ko‘rinishidagi bir o‘lchovli Gallice sxemasi maxsus to‘lqin tezligi baholaridan foydalanuvchi HLLC yechuvchisiga katta darajada ekvivalent ekani ta’kidlanadi.
Ko‘p o‘lchovli Gallice-2D sxemasi
Ko‘p o‘lchovli sxemada oqim faqat yuza qo‘shnilaridan emas, ayni tugunni ulashadigan barcha hujayralardan ta’sirlanadi. 4-sahifadagi 1-rasm klassik yuza oqimlari bilan tugun atrofida hisoblangan ko‘p o‘lchovli oqimlarni taqqoslaydi.
Oraliq holatning normal tezligi tugunga tegishli umumiy tezlik vektorining yuza normaliga proyeksiyasi sifatida aniqlanadi:
\[ u_n^*=\mathbf{u}_p\cdot\mathbf{n} \]
Tugun tezligi up tugun atrofidagi barcha yuzalarning saqlanish shartini qanoatlantiradigan chiziqli tizimdan olinadi. Ushbu ko‘p o‘lchovli shablon sonli diffuziyani oqim va to‘r geometriyasiga bog‘liq tarzda yo‘naltira oladi.
Tadqiqotga ko‘ra, ko‘p o‘lchovli tuzilma kuchli yoy zarbalarida kuzatiladigan karbunkul beqarorligiga bir o‘lchovli yechuvchilarga qaraganda chidamliroq.
Musbatlik va barqarorlik shartlari
Gaz dinamikasi yechimida zichlik va ichki energiya manfiy bo‘lmasligi kerak. Tadqiqotda bir o‘lchovli yechuvchi uchun Lagrange to‘lqin tezliklari ma’lum quyi chegaralardan katta bo‘lsa, oraliq holatlar musbat qolishi ko‘rsatilgan.
Taklif qilingan umumiy shart:
\[ \lambda_l> \max \left( \rho_lc_l,\, \sqrt{\rho_l|\Delta p|},\, -\rho_l\Delta u_n \right) \]
\[ \lambda_r> \max \left( \rho_rc_r,\, \sqrt{\rho_r|\Delta p|},\, -\rho_r\Delta u_n \right) \]
ko‘rinishidadir. c tovush tezligini bildiradi. Hisoblarda har bir yuzada chap va o‘ng shartlardan olingan eng katta qiyalik ishlatilgan.
Vaqt qadami uchun Courant–Friedrichs–Lewy soni 0,5 deb tanlangan. Barqaror masalalarda mahalliy vaqt qadami qo‘llangan.
Häenel sharti
Umumiy entalpiyani saqlovchi sxemalarning asosiy sharti energiya oqimidagi sonli diffuziya massa oqimidagi diffuziyaga mutanosib bo‘lishidir:
\[ D_{\rho e}=hD_\rho \]
Bu yerda Dρe energiya oqimidagi, Dρ massa oqimidagi sonli diffuziyani bildiradi. Barqaror oqimda h = h∞ bo‘lganda energiya va massa oqimi orasidagi fizik bog‘lanish saqlanadi.
Klassik Gallice yechuvchisining oraliq energiya holatlari bu tenglikni qanoatlantirmaydi. Oraliq holatlar orqali hisoblangan energiya farqi zichlik farqi bilan umumiy entalpi ko‘paytmasiga teng kelmaydi.
MGallice o‘zgartirishi
Tadqiqotchi Gallice yechuvchisining massa va impuls tuzilishini saqlagan holda energiyaga oid oraliq holatlarni umumiy entalpi orqali qayta ta’riflaydi. O‘zgartirilgan Lagrange holat vektori:
\[ \widetilde{\mathbf{V}}_s= \begin{pmatrix} v_s \\ \mathbf{u}_s \\ h_s \end{pmatrix} \]
ko‘rinishida quriladi. v solishtirma hajmni, u tezlikni va h umumiy entalpiyani bildiradi.
Euler tizimidagi o‘zgartirilgan oraliq holatning energiya tarkibi:
\[ \widetilde{\mathbf{U}}_s^* = \rho_s^* \begin{pmatrix} 1 \\ \mathbf{u}_s+(u_n^*-u_{n,s})\mathbf{n} \\ h_s \end{pmatrix} \]
shaklida hosil qilinadi. Shu tariqa har bir to‘lqin bo‘ylab energiya tarkibidagi o‘zgarish zichlik o‘zgarishi bilan ayni umumiy entalpi orqali bog‘lanadi.
Barqaror yechimda chap va o‘ng umumiy entalpiyalar h∞ qiymatiga teng bo‘lsa:
\[ \rho_l^*e_l^*-\rho_lh_l = (\rho_l^*-\rho_l)h_\infty \]
va o‘xshash tengliklar boshqa oraliq to‘lqinlar uchun ham bajariladi. Bu tuzilma Häenel shartini bir va ko‘p o‘lchovli o‘zgartirilgan Gallice yechuvchilari uchun qanoatlantiradi.
O‘zgartirishning muhim narxi
Energiya oqimining Lagrange ko‘rinishidagi mos tarkibi o‘zgartirishdan keyin diffuziyasiz markaziy oqimga aylanadi:
\[ \overline{pu} = \frac{1}{2} \left( p_lu_{n,l} + p_ru_{n,r} \right) \]
Massa va impuls tenglamalarida yetarli sonli diffuziya mavjudligi sinovlarda to‘g‘ri yechimni ushlash uchun yetarli bo‘lgan. Biroq muallif o‘zgartirilgan MGallice yechuvchilarining entropiyani oshiruvchi yoki entropiya barqaror ekanini kafolatlaydigan shart bera olmaydi.
Shu sabab usul umumiy entalpiyani saqlasa-da, klassik sxemaning isbotlangan entropiya xususiyatining barchasini avtomatik ravishda meros qilib olmaydi. Bu tadqiqotning eng muhim nazariy cheklovlaridan biridir.
Ikkinchi tartibli MUSCL kengaytmasi
Birinchi tartibli chekli hajm sxemalari kuchli zarbalarda barqaror bo‘lsa-da, silliq hududlarda yuqori sonli diffuziya hosil qiladi. Tadqiqotda ikkinchi tartibli yechim uchun MUSCL turidagi qayta qurish ishlatilgan.
- Zichlik, tezlik va bosim gradientlari tugunlarda vaznlangan chiziqli eng kichik kvadratlar usuli bilan hisoblangan.
- Hujayra markazi gradientlari tugun gradientlari normlarining teskarisiga asoslangan vaznlar bilan birlashtirilgan.
- MGallice sxemalarida umumiy entalpi ham qayta qurilgan.
- Gradientlar R3 Nishikawa cheklovchisi bilan cheklangan.
- Cheklovchida p = 3 va hujayra o‘lchamiga bog‘liq εc ishlatilgan.
Taqqoslangan to‘rtta sxema
| Sxema | Oqim shabloni | Umumiy entalpi o‘zgartirishi |
|---|---|---|
| Gallice | Yuza qo‘shnichiligiga asoslangan bir o‘lchovli | Yo‘q |
| Gallice-2D | Tugun qo‘shnichiligiga asoslangan ko‘p o‘lchovli | Yo‘q |
| MGallice | Yuza qo‘shnichiligiga asoslangan bir o‘lchovli | Bor |
| MGallice-2D | Tugun qo‘shnichiligiga asoslangan ko‘p o‘lchovli | Bor |
Birinchi sinov: Mach 2 tezlikdagi uchburchak bo‘rtmali kanal
Birinchi masala qarama-qarshi devorlarida kichik uchburchak bo‘rtmalar bo‘lgan ikki o‘lchovli kanaldagi tovushdan tez oqimdir. Soha 2 metr uzunlikda va 1 metr balandlikda. Ikki bo‘rtma kanalning pastki va yuqori devorlariga joylashtirilgan.
| Erkin oqim Mach soni | 2 |
|---|---|
| Hisoblash sohasi | [0, 2] m × [−0,5, 0,5] m |
| To‘r | 200 × 100 to‘rtburchak hujayra |
| Boshlang‘ich va kirish holati | ρ = 1, ux = M∞√γ, uy = 0, p = 1 |
| Solishtirma issiqliklar nisbati | γ = 7/5 |
| Kirish umumiy entalpiyasi | h∞ = 6,3 |
| Tartib | Ikkinchi tartibli MUSCL |
| Yaqinlashish mezoni | Zichlik qoldig‘ining 10−8 darajasiga tushishi |
Oqim bo‘rtmalarga urilganda ularning oldida ikkita turg‘un zarba hosil bo‘ladi. Ularning ortida rivojlangan kengayish to‘lqinlari bir-biri bilan o‘zaro ta’sirlashadi va yuqori-pastki devorlardan qaytib, olmos shaklidagi to‘lqin tuzilmasini hosil qiladi.
9-sahifadagi 3-rasm barcha sxemalar taxminan 4.000 iteratsiyadan so‘ng barqaror yechimga yaqinlashganini ko‘rsatadi. Qoldiq egri chiziqlari bir-biriga yaqin bo‘lgani sababli entalpi o‘zgartirishi yaqinlashish qobiliyatini sezilarli yomonlashtirmagan.
Zichlik maydonlari bir-biriga yaqin edi
10-sahifadagi 4-rasm to‘rtta yechuvchining 0,63 bilan 1,77 oralig‘idagi zichlik izochiziqlarini ko‘rsatadi. Barcha yechuvchilar zarbalar, kengayish to‘lqinlari va olmos shaklidagi o‘zaro ta’sir hududini o‘xshash tarzda qayta hosil qilgan.
Bir o‘lchovli sxemalar biroz kamroq diffuziv ekani aytiladi. Biroq faqat zichlik maydoniga qaralganda klassik va o‘zgartirilgan yechuvchilar orasidagi asosiy farq yaqqol ko‘rinmaydi.
Entalpi maydonlari sezilarli farq qildi
10-sahifadagi 5-rasm umumiy entalpi maydonini 6,2 bilan 6,44 oralig‘ida ko‘rsatadi. Klassik Gallice sxemasida turg‘un zarbalar ustida kichik entalpi sakrashlari yuzaga keladi va bu xatolar devor bo‘ylab oqim quyi qismiga tashiladi.
Gallice-2D sxemasidagi entalpi buzilishi yanada aniqroq bo‘ladi. Aksincha, MGallice va MGallice-2D maydonlari bir xil ko‘rinadi; kirish umumiy entalpiyasi bo‘lgan 6,3 qiymati butun maydonda saqlangani bildiriladi.
Ushbu sinov zichlik maydonlari o‘xshash ko‘rinsa ham, yechuvchilarning energiyani saqlash xususiyatlari sezilarli darajada farq qilishi mumkinligini ko‘rsatadi.
Ikkinchi sinov: Mach 20 tezlikdagi yarim silindr
Ikkinchi masala radiusi 1 metr bo‘lgan yarim silindr atrofidagi barqaror gipertovushli oqimdir. Ushbu masala ham kuchli ajralgan yoy zarbasi, ham karbunkul beqarorligi sabab qiyin hisoblanadi.
| Erkin oqim Mach soni | 20 |
|---|---|
| Silindr radiusi | 1 m |
| Hisoblash sohasi | Ellips shaklidagi soha |
| To‘r | 32 × 128 to‘rtburchak hujayra |
| Boshlang‘ich va kirish holati | ρ = 1, ux = M∞√γ, uy = 0, p = 1 |
| Solishtirma issiqliklar nisbati | γ = 7/5 |
| Kirish umumiy entalpiyasi | h∞ = 283,5 |
| Dastlabki taqqoslash | Birinchi tartibli sxemalar |
Iteratsiya sonidagi hisobot nomuvofiqligi
Matnda birinchi tartibli yechimlar 105 iteratsiya davomida ishlatilgani va shu nuqtada qoldiq uch-to‘rt tartibga kamaygani aytiladi. Biroq 10-sahifadagi 6-rasmning gorizontal o‘qi taxminan 11.000 iteratsiyada tugaydi.
Shu sabab ko‘rsatilgan 105 qiymati bilan grafik orasida preprint ichidagi hisobot nomuvofiqligi mavjud. Grafik ko‘p o‘lchovli sxemalar ko‘rsatilgan oraliqda pastroq va muntazamroq kamayuvchi qoldiqlar hosil qilganini ko‘rsatadi.
Karbunkul beqarorligi
Karbunkul kuchli va to‘r chiziqlari bilan mos tushgan zarbalarning ayrim zarba tutuvchi sxemalarda fizik bo‘lmagan tarzda buzilishidir. Yarim silindr oldidagi yoy zarbasining tormozlanish chizig‘i bu beqarorlik yuzaga keladigan klassik holatlardan biridir.
Tadqiqotda Gallice va MGallice nomli bir o‘lchovli sxemalar kuchli yoy zarbasini to‘g‘ri yecha olmagani aytiladi. Entalpiyani saqlash o‘zgartirishi o‘z-o‘zidan karbunkulni yo‘qotmagan.
Gallice-2D va MGallice-2D sxemalari esa ko‘p o‘lchovli oqim tuzilmalari sabab karbunkulsiz yoy zarbalarini hosil qilgan. Bu natija zarba barqarorligi bilan umumiy entalpiyani saqlash ikki alohida sonli xususiyat ekanini ko‘rsatadi.
Entalpiyani saqlash zichlik maydonini tuzatdi
10-sahifadagi 7-rasm to‘rtta sxemaning 1 bilan 6,3 oralig‘idagi zichlik izochiziqlarini ko‘rsatadi. Klassik Gallice-2D sxemasida eng yuqori zichlik kutilgan tormozlanish nuqtasida hosil bo‘lmaydi. Bu yoy zarbasi ortidagi energiya taqsimoti fizik bo‘lmagan tarzda o‘zgarganini ko‘rsatadi.
MGallice-2D ayni ko‘p o‘lchovli oqim shablonini ishlatgan bo‘lsa ham, entalpiyani saqlovchi energiya o‘zgartirishi orqali zichlik maksimumini tormozlanish nuqtasi atrofiga qaytaradi.
10-sahifadagi 8-rasmda klassik yechuvchilar zarba ustida umumiy entalpi pasayishini hosil qilgani ko‘rinadi. Ushbu sun’iy kamayish zarba ortidagi tovushdan past hududda energiya taqsimotini o‘zgartiradi va natijada zichlik maydonini buzadi.
MGallice va MGallice-2D yechimlarida umumiy entalpi butun maydonda mashina aniqligigacha o‘zgarmas qolgani bildirilgan.
Ikkinchi tartibli yechimning ta’siri
Yarim silindr masalasi Gallice-2D va MGallice-2D sxemalarining ikkinchi tartibli MUSCL kengaytmalari bilan qayta hisoblangan. 11-sahifadagi 9-rasm zichlik va umumiy entalpi maydonlarini yonma-yon ko‘rsatadi.
Ikkinchi tartibli qayta qurish klassik Gallice-2D sxemasidagi umumiy energiya diffuziyasini kamaytirgan. Shu sabab Gallice-2D va MGallice-2D zichlik maydonlari birinchi tartibli natijalarga qaraganda bir-biriga yaqinlashgan.
Shunga qaramay, klassik Gallice-2D umumiy entalpiyani to‘liq saqlamaslikda davom etgan, MGallice-2D esa maydonni yana kirish qiymatida o‘zgarmas tutgan. Yuqori tartib xato miqdorini kamaytirgan, biroq saqlanish xususiyatining o‘rnini bosmagan.
Tadqiqotning asosiy muhandislik xulosasi
Taklif qilingan o‘zgartirish faqat grafiklarda silliqroq umumiy entalpi maydonini hosil qilmagan. Mach 20 sinovida zarba ortidagi energiya taqsimotini tuzatib, tormozlanish hududining zichlik maydonini ham o‘zgartirgan.
Bu natija uch alohida sonli talabni ko‘rsatadi:
- Sxema kuchli zarbalarni barqaror tutishi kerak.
- Ko‘p o‘lchovli zarba beqarorliklariga yetarlicha chidamli bo‘lishi kerak.
- Barqaror yechimning umumiy entalpi kabi fizik invariantlarini diskret darajada saqlashi kerak.
MGallice-2D o‘rganilgan Mach 20 sinovida ko‘p o‘lchovli zarba barqarorligi va umumiy entalpiyani saqlashni birga ta’minlagan yagona tuzilma bo‘lgan.
Turkiya uchun ehtimoliy qo‘llanish sohalari
Usul yuqori tezlikli havo apparatlari va atmosferaga kirish masalalari uchun ishlab chiqiladigan mahalliy yoki akademik CFD dasturlarida sonli oqim tanlovini baholashga ko‘chirilishi mumkin. Ayniqsa quyidagi tadqiqot yo‘nalishlariga tegishlidir:
- Gipertovushli va tovushdan tez apparatlar aerodinamikasi,
- Raketa va yuqori tezlikli parvoz apparatlarining tashqi oqim tahlillari,
- Yoy zarbasi va tormozlanish nuqtasi hisoblari,
- Tuzilmagan to‘rlardagi chekli hajm yechuvchilari,
- Zarba beqarorligiga chidamli ko‘p o‘lchovli Riemann yechuvchilari,
- Navier–Stokes yechimlarida chegara qatlami va devor issiqlik oqimi hisoblari.
Biroq tadqiqot natijalarini bevosita apparat sirt harorati, issiqlik qalqoni qalinligi yoki aerodinamik kuchlar sifatida ishlatib bo‘lmaydi. Bu kattaliklar uchun yopishqoqlik, issiqlik o‘tkazuvchanligi, haqiqiy gaz ta’sirlari, kimyoviy reaksiyalar, turbulentlik va qattiq material xatti-harakati alohida modellashtirilishi kerak.
Tadqiqot ko‘rsatgan natijalar
- Gallice turidagi Simple Riemann yechuvchilari energiya oraliq holatlarini o‘zgartirish orqali umumiy entalpiyani saqlovchi qilinishi mumkin.
- O‘zgartirish bir va ko‘p o‘lchovli yechuvchilarga qo‘llanishi mumkin.
- Mach 2 uchburchak bo‘rtmali kanal sinovida o‘zgartirilgan sxemalar umumiy entalpiyani butun maydonda saqlagan.
- Mach 20 yarim silindr sinovida ko‘p o‘lchovli sxemalar karbunkul beqarorligini oldini olgan.
- MGallice-2D umumiy entalpi yo‘qotilishini yo‘qotib, tormozlanish hududidagi zichlik taqsimotini tuzatgan.
- Ikkinchi tartibli qayta qurish klassik sxemaning entalpi xatosini kamaytirgan, biroq butunlay yo‘qotmagan.
Tadqiqot ko‘rsatmagan natijalar
- Taklif qilingan yechuvchilar barcha gipertovushli geometriyalarda eng yaxshi usul ekani ko‘rsatilmagan.
- Entropiya barqarorligi matematik jihatdan isbotlanmagan.
- Yopishqoq Navier–Stokes tenglamalarida validatsiya o‘tkazilmagan.
- Devor issiqlik oqimi yoki chegara qatlami aniqligi bevosita o‘lchanmagan.
- Haqiqiy gaz, kimyoviy reaksiya yoki termokimyoviy nomuvozanat ta’sirlari ko‘rib chiqilmagan.
- Uch o‘lchovli geometriya yoki uch o‘lchovli Riemann yechuvchisi sinovdan o‘tkazilmagan.
- Eksperimental shamol tunneli ma’lumotlari bilan taqqoslash qilinmagan.
- Hisoblash vaqti va protsessor xarajati bo‘yicha batafsil unumdorlik taqqoslanishi berilmagan.
Tadqiqot usuli va topilmalari
Tadqiqot dizayni
Tadqiqot barqaror siqiluvchan Euler tenglamalari uchun ishlab chiqilgan bir va ko‘p o‘lchovli Godunov turidagi Simple Riemann yechuvchilarini matematik o‘zgartirish va ularni ikki sonli sinov masalasida taqqoslashga asoslangan hisoblash muhandisligi tadqiqotidir.
Matematik usul
| Asosiy tenglamalar | Ikki o‘lchovli siqiluvchan Euler tenglamalari |
|---|---|
| Gaz modeli | Ideal gaz, γ = 7/5 |
| Fazoviy diskretlashtirish | Hujayra markazli chekli hajm usuli |
| Oqim yondashuvi | Godunov turidagi Simple Riemann yechuvchilari |
| Koordinata yondashuvi | Lagrange tizimida yechuvchi qurish va Lagrange–Euler o‘zgartirishi |
| Bir o‘lchovli yechuvchi | Yuzadagi ikki hujayraga asoslangan Gallice/HLLC o‘xshash sxema |
| Ko‘p o‘lchovli yechuvchi | Tugun atrofidagi hujayralardan olingan umumiy tugun tezligiga asoslangan Gallice-2D |
| Entalpi o‘zgartirishi | Oraliq energiya holatlarini umumiy entalpi orqali qayta aniqlash |
| Vaqt qadami | Mahalliy vaqt qadami, CFL = 0,5 |
| Ikkinchi tartibli kengaytma | MUSCL qayta qurish va R3 Nishikawa cheklovchisi |
Sinov natijalarining taqqoslanishi
| Sinov | Klassik Gallice | Klassik Gallice-2D | MGallice | MGallice-2D |
|---|---|---|---|---|
| Mach 2 kanal yaqinlashuvi | Barqaror | Barqaror | Barqaror | Barqaror |
| Mach 2 zichlik maydoni | O‘xshash umumiy yechim | O‘xshash umumiy yechim | O‘xshash umumiy yechim | O‘xshash umumiy yechim |
| Mach 2 umumiy entalpi | Zarbalarda kichik sakrash | Yanada aniq buzilish | Maydon bo‘ylab saqlanadi | Maydon bo‘ylab saqlanadi |
| Mach 20 karbunkul xatti-harakati | Karbunkul hosil bo‘ladi | Karbunkul hosil bo‘lmaydi | Karbunkul hosil bo‘ladi | Karbunkul hosil bo‘lmaydi |
| Mach 20 umumiy entalpi | Zarbada kamayadi | Zarbada kamayadi | Mashina aniqligida saqlanadi | Mashina aniqligida saqlanadi |
| Mach 20 zichlik maydoni | Zarba beqarorligi sabab buzilgan | Tormozlanish nuqtasida noto‘g‘ri taqsimot | Zarba beqarorligi sabab buzilgan | Kutilgan tormozlanish nuqtasi maksimumi tiklanadi |
Sonli qiymatlar
- Mach 2 kanal to‘ri: 200 × 100 to‘rtburchak hujayra.
- Mach 2 kirish umumiy entalpiyasi: 6,3.
- Mach 2 yechimlarida zichlik qoldig‘i taxminan 4.000 iteratsiyada 10−8 darajasiga tushdi.
- Mach 2 zichlik tasvir oralig‘i: 0,63-1,77.
- Mach 2 entalpi tasvir oralig‘i: 6,2-6,44.
- Mach 20 yarim silindr to‘ri: 32 × 128 to‘rtburchak hujayra.
- Mach 20 kirish umumiy entalpiyasi: 283,5.
- Mach 20 zichlik tasvir oralig‘i: 1-6,3.
- Mach 20 entalpi tasvir oralig‘i: 240-330.
- Mach 20 matnida qoldiq uch-to‘rt tartib kamaygani bildirilgan.
Rasmlarning muhandislik talqini
- 4-sahifa, 1-rasm: Bir o‘lchovli yuza oqimlari bilan tugunga asoslangan ko‘p o‘lchovli oqim shablonini taqqoslaydi.
- 6-sahifa, 2-rasm: Lagrange tizimidagi to‘rt holatli Simple Riemann yechuvchisini va to‘lqin tezliklarini ko‘rsatadi.
- 9-sahifa, 3-rasm: Mach 2 kanal masalasida to‘rtta sxemaning qoldiq egri chiziqlari taxminan 4.000 iteratsiyada yaqinlashganini ko‘rsatadi.
- 10-sahifa, 4-rasm: Mach 2 masalasidagi zichlik to‘lqin tuzilmalari to‘rtta sxemada umuman o‘xshash ekanini ko‘rsatadi.
- 10-sahifa, 5-rasm: Entalpi xatosi klassik sxemalarda zarbalardan oqim quyi qismiga tashilishini, o‘zgartirilgan sxemalarda esa yo‘qolishini ko‘rsatadi.
- 10-sahifa, 6-rasm: Mach 20 sinovining qoldiq tarixini va ko‘p o‘lchovli sxemalarning muntazamroq yaqinlashish tendensiyasini ko‘rsatadi.
- 10-sahifa, 7-rasm: Yoy zarbasi va zichlik maydonini; bir o‘lchovli yechuvchilardagi karbunkulni hamda Gallice-2D bilan MGallice-2D orasidagi tormozlanish nuqtasi farqini ko‘rsatadi.
- 10-sahifa, 8-rasm: Klassik sxemalardagi zarba sababli entalpi kamayishini va o‘zgartirilgan sxemalarning bir tekis entalpi maydonini ko‘rsatadi.
- 11-sahifa, 9-rasm: Ikkinchi tartibli qayta qurish klassik sxemadagi xatoni kamaytirishini, MGallice-2D esa umumiy entalpiyani to‘liq saqlashda davom etishini ko‘rsatadi.
Kuchli tomonlari
- Taklif qilingan o‘zgartirish bevosita fizik barqaror oqim invariantiga asoslangan.
- Bir va ko‘p o‘lchovli yechuvchilar ayni matematik doirada ko‘rib chiqilgan.
- Sxema faqat analitik tarzda berilmagan, kuchli zarbali ikki masalada sinovdan o‘tkazilgan.
- Klassik va o‘zgartirilgan yechuvchilarning zichlik va entalpi maydonlari alohida taqqoslangan.
- Mach 20 sinovi zarba barqarorligi bilan entalpiyani saqlash alohida talablar ekanini ko‘rsatgan.
- Birinchi va ikkinchi tartibli yechimlar taqqoslangan.
- MGallice-2D ning zichlik maydonidagi yaxshilanishi saqlanish xususiyati oqim yechimini o‘zgartirishi mumkinligini ko‘rsatgan.
Kritik cheklovlar
- Tadqiqot peer review’dan o‘tmagan preprintdir.
- Faqat yopishqoqliksiz Euler tenglamalari yechilgan.
- Yopishqoqlik, issiqlik o‘tkazuvchanligi va bevosita devor issiqlik oqimi yo‘q.
- Haqiqiy gaz va yuqori harorat kimyosi modellashtirilmagan.
- Entropiya barqarorligi uchun matematik kafolat berilmagan.
- Ko‘p o‘lchovli yechuvchida musbatlikni kafolatlaydigan aniq qiyalik sharti chiqarilmagan; bir o‘lchovli yechuvchidagi qiyaliklar ishlatilgan.
- Faqat ikkita sonli sinov masalasi berilgan.
- Uch o‘lchovli yechuvchi yoki uch o‘lchovli apparat geometriyasi o‘rganilmagan.
- Fizik tajriba yoki boshqa mustaqil sonli referens bilan miqdoriy xato tahlili qilinmagan.
- Mach 2 masalasida taqqoslash asosan maydon tasvirlariga tayangan.
- Mach 20 masalasida matndagi 105 iteratsiya ifodasi bilan qoldiq grafigining taxminan 11.000 iteratsiyali o‘qi orasida nomuvofiqlik bor.
- Protsessor vaqti, xotira ishlatilishi va klassik yechuvchilarga nisbatan qo‘shimcha hisoblash xarajati bildirilmagan.
- MGallice-2D ning devor issiqlik oqimi yoki chegara qatlami aniqlash unumdorligi hali sinovdan o‘tkazilmagan.
Manba va usul eslatmasi
Asl tadqiqot: “Enthalpy preserving Simple Riemann solvers for steady hypersonic flows.”
Muallif: Lucas Tallois.
Mas’ul muallif: Lucas Tallois. Tadqiqot bitta muallifga tegishli va teng hissa bayonoti yo‘q.
Muassasa: CEA-CESTA, 15 avenue des Sablières, CS 60001, 33116 Le Barp Cedex, Fransiya.
Manba turi: Original hisoblash suyuqliklar dinamikasi tadqiqoti va preprint.
Preprint sanasi: 12 Haziran 2026.
Peer review holati: Tadqiqot peer review’dan o‘tmagan. Natijalar yakuniy va peer-reviewed jurnal maqolasi topilmalari sifatida baholanmasligi kerak.
Jurnal holati: Matnda “Preprint submitted to Elsevier” iborasi bor. Muayyan jurnal nomi, qabul sanasi yoki nashr etilgan yakuniy versiya ma’lumoti yo‘q.
Platforma: SSRN.
Rasmiy havola:SSRN tadqiqot qaydi
Moliyalashtirish: Tadqiqot matnida moliyalashtirish bayonoti yo‘q.
Manfaatlar to‘qnashuvi: Tadqiqot matnida manfaatlar to‘qnashuvi bayonoti yo‘q.
Ma’lumot va kodga kirish: Ochiq ma’lumot ombori, sonli to‘r fayllari yoki manba kodi havolasi berilmagan.
Ushbu kontentdagi ilmiy usul, tenglamalar, sonli sxemalar, sinov sharoitlari, grafik talqinlari, topilmalar va cheklovlar ko‘rib chiqilgan tadqiqotga asoslanadi. Tashqi manbalar faqat sarlavha, muallif, muassasa, platforma va DOI ma’lumotlarini bibliografik tekshirish uchun ishlatilgan; tadqiqotda bo‘lmagan ilmiy natija qo‘shilmagan.
Izoh qoldiring
E-pochta manzilingiz chop etilmaydi. Majburiy maydonlar * bilan belgilangan