Pichoq elementlari momentum nazariyasi - Blade element momentum theory

Ushbu maqolada bir nechta muammolar mavjud. Iltimos yordam bering uni yaxshilang yoki ushbu masalalarni muhokama qiling munozara sahifasi. (Ushbu shablon xabarlarini qanday va qachon olib tashlashni bilib oling) Bu maqola uchun qo'shimcha iqtiboslar kerak tekshirish. Iltimos yordam bering ushbu maqolani yaxshilang tomonidan ishonchli manbalarga iqtiboslarni qo'shish. Ma'lumot manbasi bo'lmagan material shubha ostiga olinishi va olib tashlanishi mumkin. Manbalarni toping: "Blade element momentum nazariyasi"  – Yangiliklar  · gazetalar  · kitoblar  · olim  · JSTOR (2013 yil may) (Ushbu shablon xabarini qanday va qachon olib tashlashni bilib oling) Bu maqola aksariyat o'quvchilar tushunishi uchun juda texnik bo'lishi mumkin. Iltimos uni yaxshilashga yordam bering ga buni mutaxassis bo'lmaganlarga tushunarli qilish, texnik ma'lumotlarni olib tashlamasdan. (2013 yil mart) (Ushbu shablon xabarini qanday va qachon olib tashlashni bilib oling) (Ushbu shablon xabarini qanday va qachon olib tashlashni bilib oling)

Pichoq elementlari momentum nazariyasi ikkalasini ham birlashtirgan nazariya pichoq elementlari nazariyasi va momentum nazariyasi. Bu pervanel yoki shamol turbinasi pichog'idagi mahalliy kuchlarni hisoblash uchun ishlatiladi. Pichoq elementlari nazariyasi momentum nazariyasi bilan birlashtirilib, rotorda induksion tezlikni hisoblashda ba'zi qiyinchiliklarni engillashtiradi.

Ushbu maqola BEMni erga asoslangan shamol turbinalariga qo'llashni ta'kidlaydi, ammo printsiplar pervanellarga ham tegishli. Streamtube maydoni pervanel tomonidan qisqartirilgan bo'lsa, u shamol turbinasi bilan kengaytiriladi. Har qanday dastur uchun juda soddalashtirilgan, ammo foydali yaqinlashish Rankine-Froude "momentum" yoki "aktuator disk" modelidir (1865,1889). Ushbu maqolada "Betz limiti" ning er usti shamol turbinasi samaradorligiga qo'llanilishi tushuntirilgan.

Rivojlanish Frudning pichoq elementi momentum nazariyasi (1878) shaklida bo'lib, keyinchalik Glauert tomonidan takomillashtirilgan (1926). Betz (1921) aktuator diskining (NACA TN 83, "Vida pervanesi nazariyasi" va NACA TM 491) oqimiga to'satdan aylanishini hisobga olish uchun "Rankine-Froude aktuator-disk" nazariyasiga taxminiy tuzatish kiritdi. , "Pervanel muammolari"). Pichoq elementlari momentum nazariyasida burchak impulsi modelga kiritilgan, ya'ni uyg'onish (rotor bilan o'zaro ta'sirdan keyingi havo) burchak impulsiga ega. Ya'ni, rotor bilan o'zaro ta'sirlashganda havo z o'qi atrofida darhol aylana boshlaydi (quyidagi diagramaga qarang). Burchak momentumini hisobga olish kerak, chunki shamoldan energiya chiqaradigan moslama bo'lgan rotor shamol bilan o'zaro ta'siri natijasida aylanadi.

Rankine-Froude modeli

"Betz limiti" hali pervanellarni ta'kidlab, aylanish oqimini hisobga olish uchun Betzning hissasidan foydalanmayapti, Rankine-Froude-ni qo'llaydi " aktuator disk "turg'un shamol turbinasining maksimal samaradorligini olish nazariyasi. Quyidagi tahlil havoning eksenel harakati bilan cheklangan:

Bizda streamtube bizda chapdan o'ngga suyuqlik oqadi va rotorni ifodalovchi aktuator disk mavjud. Rotor cheksiz darajada ingichka deb taxmin qilamiz.^[1] Yuqoridan, oqim trubasining boshida, aktuator diskida suyuqlik oqimi normal ekanligini ko'rishimiz mumkin. Suyuqlik rotor bilan o'zaro ta'sir qiladi va shu bilan energiyani suyuqlikdan rotorga uzatadi. Keyin suyuqlik oqim bo'ylab oqishni davom ettiradi. Shunday qilib biz tizimimiz / streamtube-ni ikkita qismga ajratishimiz mumkin: oldingi akuator disk va post-aktuator disk. Rotor bilan o'zaro ta'sir qilishdan oldin suyuqlikdagi umumiy energiya doimiydir. Bundan tashqari, rotor bilan ta'sir o'tkazgandan so'ng, suyuqlikdagi umumiy energiya doimiy bo'ladi.

Bernulli tenglamasi aniq energiya doimiy bo'lgan suyuqlik oqimida mavjud bo'lgan turli xil energiya shakllarini tavsiflaydi, ya'ni suyuqlik rotor kabi boshqa biron bir narsaga energiya o'tkazmasa. Energiya statik bosim, tortishish potentsiali energiyasi va kinetik energiyadan iborat. Matematik jihatdan bizda quyidagi ibora mavjud:

${ displaystyle { frac {1} {2}} rho v ^ {2} + P + rho gh = mathrm {const.}}$

qayerda ${ displaystyle rho}$ suyuqlikning zichligi, ${ displaystyle v}$ suyuqlikning oqim yo'nalishi bo'yicha tezligi, ${ displaystyle P}$ statik bosim energiyasi, ${ displaystyle g}$ tortishish kuchi tufayli tezlanish va ${ displaystyle h}$ erdan balandlik. Ushbu tahlil maqsadida, biz suyuqlik quyida chapdan o'ngga qarab tortishish jarayonida potentsial energiya o'zgarmas bo'ladi, deb hisoblaymiz, shunda biz quyidagilarga egamiz:

${ displaystyle { frac {1} {2}} rho v ^ {2} + P = mathrm {const.}}$

Shunday qilib, agar bizda oqim chizig'ida ikkita nuqta bo'lsa, 1-nuqta va 2-nuqta, va 1-nuqtada oqim chizig'i bo'ylab suyuqlikning tezligi ${ displaystyle v_ {1}}$ va 1dagi bosim ${ displaystyle P_ {1}}$, va 2 nuqtada suyuqlikning oqim chizig'i bo'ylab tezligi ${ displaystyle v_ {2}}$ va 2dagi bosim ${ displaystyle P_ {2}}$, va 1 va 2 nuqtalar orasidagi suyuqlikdan energiya olinmagan bo'lsa, unda biz quyidagi ifodaga egamiz:

${ displaystyle { frac {1} {2}} rho v_ {1} ^ {2} + P_ {1} = { frac {1} {2}} rho v_ {2} ^ {2} + P_ {2}}$

Endi dastlabki sxemamizga qaytamiz. Pre-aktuator oqimini ko'rib chiqing. Suyuqlikning tezligi yuqori oqimda ${ displaystyle v _ { infty}}$; suyuqlik keyinchalik rotorga yaqinlashganda kengayadi.^{[iqtibos kerak ]} Ommaviy konservatsiyaga muvofiq, massa oqim tezligi doimiy bo'lishi kerak. Ommaviy oqim tezligi, ${ displaystyle { dot {m}}}$, maydon yuzasi orqali ${ displaystyle A}$ quyidagi ifoda bilan berilgan:

${ displaystyle { dot {m}} = rho Av}$

qayerda ${ displaystyle rho}$ zichligi va ${ displaystyle v}$ suyuqlikning oqim yo'nalishi bo'yicha tezligi. Shunday qilib, agar massa oqim tezligi doimiy bo'lsa, maydonning ko'payishi oqim yo'nalishi bo'ylab suyuqlik tezligining pasayishiga olib kelishi kerak. Bu suyuqlikning kinetik energiyasi pasayayotganligini anglatadi. Agar oqim kengayib borayotgan bo'lsa, lekin energiya uzatmasa, u holda Bernulli amal qiladi. Shunday qilib kinetik energiyaning pasayishiga statik bosim energiyasining ortishi qarshi turadi.

Shunday qilib, bizda rotorgacha bo'lgan holat mavjud: yuqori oqimda, suyuqlik bosimi atmosfera bilan bir xil, ${ displaystyle P _ { infty}}$; rotor bilan o'zaro ta'sir qilishdan oldin suyuqlik bosimi oshdi va shuning uchun kinetik energiya kamaydi. Buni Bernulli tenglamasidan foydalangan holda matematik tarzda tasvirlash mumkin:

${ displaystyle { frac {1} {2}} rho v _ { infty} ^ {2} + P _ { infty} = { frac {1} {2}} rho left (v _ { infty) } (1-a) o'ng) ^ {2} + P_ {D +}}$

bu erda biz rotorda suyuqlik tezligini quyidagicha yozdik ${ displaystyle v _ { infty} (1-a)}$, qayerda ${ displaystyle a}$ eksenel induksiya faktoridir. Suyuqlikning aktuator diskining yuqori tomonidagi bosimi ${ displaystyle P_ {D +}}$. Biz rotorni cheksiz ingichka bo'lgan aktuator disk sifatida ko'rib chiqmoqdamiz. Shunday qilib, biz aktuator diskida suyuqlik tezligida hech qanday o'zgarish bo'lmaydi. Suyuqlikdan energiya olinganligi sababli, bosim kamaygan bo'lishi kerak.

Endi post-rotorni ko'rib chiqing: rotor bilan o'zaro aloqada bo'lganidan so'ng, suyuqlik tezligi hali ham saqlanib qoladi ${ displaystyle v _ { infty} (1-a)}$, ammo bosim qiymatga tushdi ${ displaystyle P_ {D-}}$; quyi oqimda suyuqlikning bosimi atmosfera bilan muvozanatga erishdi; Bu dinamik muvozanatni saqlash uchun oqim trubkasidagi oqim tezligini kamaytirishning tabiiy va dinamik ravishda sekin jarayonida amalga oshirildi (ya'ni. ${ displaystyle P rightarrow P _ { infty}}$ ancha quyi oqimda. Boshqa energiya uzatishni nazarda tutmasak, Bernulliga quyi oqim uchun murojaat qilishimiz mumkin:

${ displaystyle { frac {1} {2}} rho chap (v _ { infty} (1-a) o'ng) ^ {2} + P_ {D -} = { frac {1} {2 }} rho v_ {w} ^ {2} + P _ { infty}}$

qayerda

${ displaystyle v_ {w} =}$ Uyg'otishning quyi oqimidagi tezlik

Shunday qilib, biz rotorning old va orqa tomonlari orasidagi bosim farqining ifodasini olishimiz mumkin:

${ displaystyle P_ {D +} - P_ {D -} = { frac {1} {2}} rho (v _ { infty} ^ {2} -v_ {w} ^ {2})}$

Agar biz aktuator disk maydoni bo'yicha bosim farqiga ega bo'lsak, aktuator diskida ta'sir qiluvchi kuch bor, uni aniqlash mumkin ${ displaystyle F = Delta PA}$:

${ displaystyle { frac {1} {2}} rho (v _ { infty} ^ {2} -v_ {w} ^ {2}) A_ {D}}$

qayerda ${ displaystyle A_ {D}}$ - bu aktuator diskining maydoni. Agar rotor suyuqlikdan energiyani yutadigan yagona narsa bo'lsa, suyuqlikning eksenel momentumining o'zgarish tezligi rotorga ta'sir qiladigan kuchdir. Eksenel momentumning o'zgarishi tezligi suyuqlikning boshlang'ich va oxirgi eksenel tezliklari orasidagi farq sifatida massa oqim tezligiga ko'paytirilishi mumkin:

${ displaystyle F = { frac { mathrm {d} p} { mathrm {d} t}} = { nuqta {m}} (v _ { infty} -v_ {w}) = rho A_ { D} v_ {D} (v _ { infty} -v_ {w}) = rho A_ {D} (1-a) v _ { infty} (v _ { infty} -v_ {w})}$

Shunday qilib biz quyi oqim bo'ylab suyuqlik tezligining ifodasiga erishishimiz mumkin:

${ displaystyle v_ {w} = (1-2a) v _ { infty}}$

Ushbu kuch rotorda ishlaydi. Suyuqlikdan olinadigan quvvat - bu suyuqlikka ta'sir qiladigan kuch, quvvatni olish nuqtasida suyuqlikning tezligiga ko'paytiriladi:

${ displaystyle mathrm {Power} _ {ext} = Fv_ {D} = 2a (1-a) ^ {2} v _ { infty} ^ {3} rho A_ {D}}$

Maksimal quvvat

Aytaylik, biz suyuqlikdan olinadigan maksimal quvvatni topishga qiziqmoqdamiz. Suyuqlikdagi quvvat quyidagi ifoda bilan berilgan:

${ displaystyle mathrm {Power} = { frac {1} {2}} rho Av ^ {3}}$

qayerda ${ displaystyle rho}$ oldingi kabi suyuqlik zichligi, ${ displaystyle v}$ suyuqlik tezligi va ${ displaystyle A}$ suyuqlik oqayotgan xayoliy sirt maydoni. Yuqorida tavsiflangan stsenariyda rotor yordamida suyuqlikdan chiqarilgan quvvat bu quvvat ifodasining ba'zi bir qismidir. Biz fraktsiyani quvvatni samarali deb ataymiz, ${ displaystyle C_ {p}}$. Shunday qilib, chiqarilgan quvvat, ${ displaystyle mathrm {Power} _ {ext}}$ quyidagi ifoda bilan berilgan:

${ displaystyle mathrm {Power} _ {ext} = mathrm {Power} times C_ {p}}$

Bizning savolimiz shu: maksimal qiymati nimaga teng ${ displaystyle C_ {p}}$ Betz modelidan foydalanayapsizmi?

Suyuqlikdan rotorga o'tkaziladigan quvvat uchun olingan ifodamizga qaytaylik (${ displaystyle mathrm {Power} _ {ext}}$). Chiqarilgan quvvat eksenel induksiya faktoriga bog'liqligini ko'rishimiz mumkin. Agar biz farq qilsak ${ displaystyle mathrm {Power} _ {ext}}$ munosabat bilan ${ displaystyle a}$, biz quyidagi natijani olamiz:

${ displaystyle { frac { mathrm {d} mathrm {Power} _ {ext}} { mathrm {d} a}} = 2v _ { infty} ^ {3} rho A_ {D} times chap ((1-a) ^ {2} -2a (1-a) o'ng)}$

Agar biz quvvatni chiqarishni maksimal darajaga ko'targan bo'lsak, yuqoridagi narsani nolga o'rnatishimiz mumkin. Bu bizga qiymatini aniqlashga imkon beradi ${ displaystyle a}$ bu maksimal quvvatni chiqarishni ta'minlaydi. Ushbu qiymat a ${ displaystyle 1/3}$. Shunday qilib, biz buni topa olamiz ${ displaystyle C_ {P ~ max} = 16/27}$. Boshqacha qilib aytganda, rotor suyuqlikdagi quvvatning 59 foizidan ko'pini ololmaydi.

Pichoq elementlari momentum nazariyasi

Rankine-Froude modeli bilan taqqoslaganda, Blade element momentum nazariyasi rotorning burchak momentumini hisobga oladi. Quyidagi rasmning chap tomonini ko'rib chiqing. Bizda streamtube bor, unda suyuqlik va rotor bor. Streamtube tarkibidagi narsalar va uning tashqarisidagi narsalar o'rtasida o'zaro ta'sir yo'q deb o'ylaymiz. Ya'ni, biz izolyatsiya qilingan tizim bilan shug'ullanmoqdamiz. Fizikada ajratilgan tizimlar saqlash qonunlariga bo'ysunishi kerak. Bunga burchak momentumining saqlanishi misol bo'la oladi. Shunday qilib, oqim trubkasidagi burchak momentumini saqlash kerak. Binobarin, agar rotor suyuqlik bilan o'zaro ta'siri orqali burchak impulsiga ega bo'lsa, yana bir narsa teng va qarama-qarshi burchak momentumiga ega bo'lishi kerak. Yuqorida aytib o'tilganidek, tizim faqat suyuqlik va rotordan iborat, suyuqlik uyg'otish vaqtida burchak impulsiga ega bo'lishi kerak. Eksenel impulsning o'zgarishini ba'zi bir indüksiyon faktori bilan bog'laganimiz kabi ${ displaystyle a}$, suyuqlikning burchak momentumining o'zgarishini tangensial induksiya faktori bilan bog'laymiz, ${ displaystyle a '}$.

Quyidagi sozlamalarni ko'rib chiqing.^[1]

Biz rotor maydonini cheksiz kichik qalinlikdagi halqali halqalarga ajratamiz. Biz buni eksenel induksiya omillari va tangensial indüksiyon omillari halqa halqasi davomida doimiy deb taxmin qilishimiz uchun qilyapmiz. Ushbu yondashuvning taxmin qilishicha halqali halqalar bir-biridan mustaqil, ya'ni qo'shni halqa halqalari suyuqliklari o'rtasida o'zaro ta'sir bo'lmaydi.

Bernulli aylanadigan uyg'onish uchun

Endi Bernulliga qaytaylik:

${ displaystyle { frac {1} {2}} rho v_ {1} ^ {2} + P_ {1} = { frac {1} {2}} rho v_ {2} ^ {2} + P_ {2}}$

Tezlik - bu suyuqlikning oqim yo'nalishi bo'yicha tezligi. Oqim chizig'i ma'lum bir koordinata o'qiga, masalan, z o'qiga parallel ravishda o'tishi shart emas. Shunday qilib tezlik koordinata tizimini tashkil etuvchi o'qlar tarkibidagi qismlardan iborat bo'lishi mumkin. Ushbu tahlil uchun biz silindrsimon qutb koordinatalarini ishlatamiz ${ displaystyle (r, ~ theta, ~ z)}$. Shunday qilib ${ displaystyle v ^ {2} = v_ {r} ^ {2} + v _ { theta} ^ {2} + v_ {z} ^ {2}}$.

Izoh: Biz, aslida, barcha jihatlar bo'yicha silindrsimon koordinatalarda ishlaymiz. ${ displaystyle mathbf {F} = F_ {r} { hat { mathbf {r}}} + F _ { theta} { hat { theta}} + F_ {z} { hat { mathbf { z}}}}$

Endi yuqorida ko'rsatilgan sozlamani ko'rib chiqing. Oldingi kabi biz o'rnatishni ikkita komponentga ajratishimiz mumkin: yuqori va quyi oqimlar.

Old rotor

${ displaystyle P _ { infty} + { frac {1} {2}} rho v_ {u} ^ {2} = P_ {D +} + { frac {1} {2}} rho v_ {D } ^ {2}}$

qayerda ${ displaystyle v_ {u}}$ bu oqimning ancha yuqorisidagi oqim chizig'i bo'ylab suyuqlikning tezligi va ${ displaystyle v_ {D}}$ suyuqlikning rotor oldidan tezligi. Silindrsimon qutb koordinatalarida yozilgan bizda quyidagi ibora mavjud:

${ displaystyle P _ { infty} + { frac {1} {2}} rho v _ { infty} ^ {2} = P_ {D +} + { frac {1} {2}} rho (v_) { infty} (1-a)) ^ {2}}$

qayerda ${ displaystyle v _ { infty}}$ va ${ displaystyle v _ { infty} (1-a)}$ tezlikning z-komponentlari yuqoriga qarab va rotordan biroz oldinroq. Bu Betz modelidan oldingi oqim tenglamasi bilan bir xil.

Yuqoridagi rasmdan ko'rinib turibdiki, oqim rotorga yaqinlashganda kengayadi, bu statik bosimning oshishi va massaning saqlanishi natijasidir. Bu shuni anglatadiki ${ displaystyle v_ {r} neq 0}$ yuqori oqim. Biroq, ushbu tahlil qilish uchun ushbu ta'sir beparvo qilinadi.

Post-rotor

${ displaystyle P_ {D -} + { frac {1} {2}} rho v_ {D} ^ {2} = P _ { infty} + { frac {1} {2}} rho v_ { w} ^ {2}}$

qayerda ${ displaystyle v_ {D}}$ rotor bilan o'zaro aloqada bo'lganidan keyin suyuqlikning tezligi. Buni shunday yozish mumkin ${ displaystyle v_ {D} ^ {2} = v_ {D, ~ r} ^ {2} + v_ {D, ~ theta} ^ {2} + v_ {D, ~ z} ^ {2}}$. Tezlikning radial komponenti nolga teng bo'ladi; halqali halqali yondashuvni qo'llasak, bu to'g'ri bo'lishi kerak; Aks holda taxmin qilish, quyi oqimning bir qismida halqali halqalar orasidagi aralashuvni taklif qiladi. Disk bo'ylab eksenel tezlikda o'zgarish yo'q deb o'ylaganimiz sababli, ${ displaystyle v_ {D, ~ z} = (1-a) v _ { infty}}$. Burchak impulsini izolyatsiya qilingan tizimda saqlash kerak. Shunday qilib, uyg'otishning aylanishi to'xtab qolmasligi kerak. Shunday qilib ${ displaystyle v _ { theta}}$ quyi oqim qismida doimiy bo'ladi. Shunday qilib, Bernulli quyi oqim qismida soddalashtiradi:

${ displaystyle P_ {D -} + { frac {1} {2}} rho v_ {D, ~ z} ^ {2} = P _ { infty} + { frac {1} {2}} rho v_ {w, ~ z} ^ {2} = P_ {D -} + { frac {1} {2}} rho (v _ { infty} (1-a)) ^ {2}}$

Boshqacha qilib aytganda, rotorning yuqorisida va pastida Bernulli tenglamalari Betz modelidagi Bernulli ifodalari bilan bir xil. Shuning uchun biz Betz modelida olingan quvvatni chiqarish va uyg'onish tezligi kabi natijalardan foydalanishimiz mumkin, ya'ni.

${ displaystyle v_ {w, z} = (1-2a) v _ { infty}}$

${ displaystyle mathrm {Power} = 2a (1-a) ^ {2} v _ { infty} ^ {3} rho A_ {D}}$

Bu bizga aylanadigan uyg'otishni o'z ichiga olgan tizim uchun maksimal quvvat sarfini hisoblash imkonini beradi. Buni Betz modeli bilan bir xil qiymatni berish uchun ko'rsatish mumkin, ya'ni 0,59. Ushbu usul rotorda hosil bo'lgan moment quyidagi ifoda bilan berilganligini tan olishni o'z ichiga oladi:

${ displaystyle delta mathbf {Q} = 2 pi r delta r times rho U _ { infty} (1-a) times 2a'r omega}$

darhol quyida belgilangan zarur shartlar bilan.

Pichoq kuchlari

Havo plyonkasi atrofida suyuqlik oqimini ko'rib chiqing. Havo plyonkasi atrofidagi suyuqlik oqimi kuchlarni ko'tarish va tortish kuchini keltirib chiqaradi. Ta'rifga ko'ra, ko'tarish - bu havo plyonkasida ko'ringan suyuqlik oqimi tezligiga normal ravishda ta'sir qiladigan kuch. Drag - bu havo plyonkasida ko'rinadigan suyuqlik oqimi tezligiga ta'sir qiluvchi ta'sir qiluvchi kuchlar. Ko'rinib turgan tezlik deganda nimani tushunamiz? Quyidagi diagrammani ko'rib chiqing:

Rotor pichog'ining tezligi uch narsaga bog'liq: suyuqlikning eksenel tezligi, ${ displaystyle v _ { infty} (1-a)}$; plyonkaning tezlashishi tufayli suyuqlikning teginsel tezligi, ${ displaystyle a ' omega r}$; va rotor harakatining o'zi, ${ displaystyle omega r}$. Ya'ni, ko'rinadigan suyuqlik tezligi quyida keltirilgan:

${ displaystyle mathbf {v} = omega r (1 + a ') { hat { mathbf { theta}}} + v _ { infty} (1-a) { hat { mathbf {z} }}}$

Shunday qilib, aniq shamol tezligi bu vektorning kattaligi, ya'ni:

${ displaystyle | mathbf {v} | ^ {2} = ( omega r (1 + a ')) ^ {2} + (v _ { infty} (1-a)) ^ {2} = W ^ {2}}$

Shuningdek, biz burchakni ishlab chiqishimiz mumkin ${ displaystyle phi}$ yuqoridagi rasmdan:

${ displaystyle sin phi = { frac {v _ { infty} (1-a)} {W}}}$

Biz burchakni bilamiz deylik ${ displaystyle beta}$, keyin ishlashimiz mumkin ${ displaystyle alpha}$ shunchaki munosabat yordamida ${ displaystyle alpha = phi - beta}$; keyin biz liftni koeffitsient bilan ishlab chiqishimiz mumkin, ${ displaystyle c_ {L}}$va tortishish samarali ${ displaystyle c_ {D}}$, undan pichoqqa ta'sir qiluvchi ko'tarish va tortish kuchlarini ishlab chiqishimiz mumkin.

Qisman pichoq elementlari egallagan halqali halqani ko'rib chiqing. Halqa halqasini egallagan har bir pichoq qismining uzunligi ${ displaystyle delta r}$ (quyidagi rasmga qarang).

Pichoqlar / havo plyonkalarining har bir qismiga ta'sir qiluvchi ko'taruvchi akkord ${ displaystyle c}$ quyidagi ifoda bilan berilgan:

${ displaystyle delta L = { frac {1} {2}} rho NW ^ {2} c times c_ {L} ( alpha) delta r}$

qayerda ${ displaystyle c_ {L}}$ hujum koeffitsienti, bu hujum burchagi funktsiyasi va ${ displaystyle N}$ pichoqlar soni. Bundan tashqari, pichoqlar / havo plyonkalarining bu qismiga akkord bilan ta'sir qiluvchi tortishish ${ displaystyle c}$ quyidagi ifoda bilan berilgan:

${ displaystyle delta D = { frac {1} {2}} rho NW ^ {2} c times c_ {D} ( alpha) delta r}$

Esda tutingki, hisoblangan ushbu kuchlar odatdagi tezlik uchun tangensial va normaldir. Biz kuchlar bilan qiziqamiz ${ displaystyle { hat { mathbf {z}}}}$ va ${ displaystyle { hat { theta}}}$ o'qlar. Shunday qilib biz quyidagi diagrammani ko'rib chiqishimiz kerak:

Shunday qilib biz quyidagilarni ko'rishimiz mumkin:

${ displaystyle delta F _ { theta} = delta L sin phi - delta D cos phi}$

${ displaystyle delta F_ {z} = delta L cos phi + delta D sin phi}$

${ displaystyle F _ { theta}}$ rotor pichoqlarining aylanishi uchun javobgar bo'lgan kuch; ${ displaystyle F_ {z}}$ pichoqlarning egilishi uchun javobgar bo'lgan kuchdir.

Eslatib o'tamiz, ajratilgan tizim uchun tizimning aniq burchak momentumi saqlanib qoladi. Agar rotor burchak momentumini qo'lga kiritgan bo'lsa, undagi suyuqlik ham shunday bo'lishi kerak. Keling, uyg'ongan suyuqlik tangensial tezlikni oladi deb taxmin qilaylik ${ displaystyle v _ { theta} = 2a ' omega r}$. Shunday qilib, havodagi moment bilan beriladi

${ displaystyle | mathbf { delta {Q}} | = rho (2 pi r delta r) U _ { infty} (1-a) times (2 Omega a'r ^ {2}) }$

Burchak momentumini saqlab qolish orqali bu rotorning pichoqlaridagi momentni muvozanatlashtiradi; shunday qilib,

${ displaystyle { frac {1} {2}} rho W ^ {2} Nc (c_ {l} sin phi -c_ {d} cos phi) r delta r = rho (2 ) pi r delta r) U _ { infty} (1-a) marta (2 Omega a'r ^ {2})}$

${ displaystyle { frac {1} {2}} W ^ {2} Nc (c_ {l} sin phi -c_ {d} cos phi) = 4 pi U _ { infty} (1- a) times Omega a'r ^ {2}}$

Bundan tashqari, havodagi chiziqli impulsning o'zgarishi tezligi pichoqlarga ta'sir qiladigan tekislikdan tashqari bukilish kuchi bilan muvozanatlashadi, ${ displaystyle delta F_ {z}}$. Impuls nazariyasidan havodagi chiziqli impulsning o'zgarish tezligi quyidagicha:

${ displaystyle delta F_ {z} = rho (2 pi r delta r) U _ { infty} (1-a) times (v _ { infty} -v_ {w})}$

sifatida ifodalanishi mumkin

${ displaystyle delta F_ {z} = rho (4 pi r delta r) U _ { infty} ^ {2} a (1-a)}$

Buni tekislikdan tashqari bukish kuchi bilan muvozanatlash beradi

${ displaystyle { frac {1} {2}} W ^ {2} Nc (c_ {l} cos phi + c_ {d} sin phi) = rho (4 pi r delta r) U _ { infty} ^ {2} a (1-a)}$

Keling, quyidagi ta'riflarni beraylik:

${ displaystyle C_ {y} = c_ {l} sin phi -c_ {d} cos phi}$

${ displaystyle C_ {x} = c_ {l} cos phi + c_ {d} sin phi}$

Shunday qilib bizda quyidagi tenglamalar mavjud:

${ displaystyle { frac {1} {2}} W ^ {2} NcC_ {y} = 4 pi U _ { infty} (1-a) times Omega a'r ^ {2}}$

(1)

${ displaystyle { frac {1} {2}} rho W ^ {2} NcC_ {x} = 4 pi rho left [(a ' Omega r) ^ {2} + U _ { infty} ^ {2} a (1-a) right] r}$

(2)

Yuqoridagi rasmni tahlilidan ko'rish mumkin bo'lgan quyidagi tenglamaga murojaat qilaylik:

${ displaystyle sin phi = { frac {U _ { infty}} {W}} (1-a) rightarrow sin ^ {2} phi = left ({ frac {U _ { infty} } {W}} (1-a) o'ng) ^ {2}}$

(3)

Shunday qilib, ushbu uchta tenglama yordamida ba'zi bir algebraik manipulyatsiya orqali quyidagi natijani olish mumkin:^[1]

${ displaystyle { frac {a} {1-a}} = { frac {C_ {x} sigma _ {r}} {4 sin ^ {2} phi}}}$

Biz uchun ifoda olishimiz mumkin ${ displaystyle a '}$ shunga o'xshash tarzda. Bu bizga rotor va suyuqlik bilan nima sodir bo'lishini tushunishga imkon beradi. Ushbu turdagi tenglamalar keyinchalik takroriy usullar bilan echiladi.

BEM modellarining taxminlari va mumkin bo'lgan kamchiliklari

• Har bir halqali uzuk har qanday boshqa halqadan mustaqil deb taxmin qiladi.^[2]
• Uyg'onishning kengayishini hisobga olmaydi.
• Hisobga olinmaydi yo'qotishlar, ammo tuzatish omillarini kiritish mumkin.^{[iqtibos kerak ]}
• Hisobga olinmaydi yaw, buni amalga oshirish mumkin bo'lsa ham.
• Barqaror oqimga asoslangan (turbulent bo'lmagan).

Adabiyotlar