Pag-unlad ng pag-aaral sa kinamakanang aerodiniko ng kamera ng pagsisiyasat ng motor ng eroplano batay sa malaking eddy simulation

Pag-unlad ng pag-aaral sa kinamakanang aerodiniko ng kamera ng pagsisiyasat ng motor ng eroplano batay sa malaking eddy simulation

Ang kamber ng pagsisilaw ay isa sa mga pangunahing bahagi ng isang propulsyon ng eroplano, at ang anyo ng aerodinamiko ng kamber ng pagsisilaw ay naglalaro ng mahalagang papel sa pagganap ng buong propulsyon. Upang tugunan ang mga karagdagang makatitserong pangangailangan teknikal ng propulsyon para sa kamber ng pagsisilaw, ang organisasyon ng pamamaraan ng pagsisilaw at ang mga characteristics ng pamumuhunan sa loob ng kamber ng pagsisilaw ay naging napakalikwat. Ang proseso ng pagpabagal at pagdaddaan ng presyo ng diffuser ay maaaring makaharap sa paghiwa ng pamumuhunan sa ilalim ng malakas na mapagpakiramdam na presyo gradient; ang hangin ay dumadaan sa multistage swirl device upang bumuo ng malaking-scale vortex structure, na sa isang banda ay humahaya sa atomization at paguubos ng likido fuel at bumubuo ng malakas na pulsating, unsteady mixture kasama ang fuel, at sa kabilang banda ay nagbubuo ng stationary flame sa aerodynamic recirculation zone; ang maramihang jets ng pangunahing combustion/mixing hole ay sumasangguni sa lateral flow sa loob ng flame tube upang bumuo ng counter-rotating vortex pair, na mayroong mahalagang impluwensya sa turbulent mixing. Sa pundasyon ng pamumuhunan, ang multiscale physical at chemical processes tulad ng atomization at evaporation, mixing, chemical reaction, at interaction pagitan ng turbulence at flame ay malakas na coupled, na pareho ay tumutukoy sa aerodinamiko characteristics ng kamber ng pagsisilaw. Ang high-precision modeling at pagkuha ng resulta ng mga physical at chemical processes ay laging isang maikling paksa ng pag-aaral sa loob at labas ng bansa.

Ang mga proseso ng atomization, evaporation, mixing, at combustion sa loob ng combustion chamber ay umuunlad at umaabot sa isang kapaligiran ng turbulent flow, kaya ang pamumuhunan ay ang pundasyon para sa pag-simul ng aerodynamic na pagganap ng combustion chamber. Ang pangunahing characteristics ng turbulence ay ang pagpapakita ng random na pulsation ng mga parameter ng pamumuhunan dahil sa nonlinear na proseso ng convection. Ang turbulence ay naglalaman ng maraming vortex structures. Ang mga sakop ng iba't ibang vortex sa mga length at time scales ay malaki, at habang tumataas ang Reynolds number, dumadagdag nang masusing ang mga sakop sa pagitan ng mga scale. Ayon sa proporsyon ng mga vortex structure ng turbulence na direktang inisyal na solusyon, ang simulasyon ng turbulence  ang mga paraan ay nahahati sa direktang numerikal na simulasyon (DNS), Reynolds-Averaged Navier-Stokes (RANS), malaking eddy simulasyon (LES) at mga halal na paraan ng simulasyon ng turbulensiya. Ang paraan ng RANS, na madalas gamitin sa inhinyeriya, ay naglutas ng turbolenteng mean field at gumagamit ng isang modelo upang simulan ang lahat ng impormasyong pagpapatakbo ng turbulensiya. Maliit ang dami ng pagkuha, ngunit mababa ang katumpakan. Para sa malakas na pagiging swirl at hindi pamilyar na proseso ng pamumuhunan sa kamera ng pagsisilaw, hindi makakamit ng RANS ang mga kinakailangan ng detalyadong disenyo. Tinukoy ni Pitsch na ang komplikadong pagkalkula ng LES ay nasa gitna ng RANS at DNS, at kasalukuyan ay ginagamit para sa pagkalkula ng turbolenteng pagsisilaw sa walang hanggang espasyo na may medium at mababang bilang ng Reynolds. Dahil sa maliit na skalang turbulensiya sa karatig-kuta ng kamera ng pagsisilaw at ang mataas na bilang ng Reynolds ng pamumuhunan, ang dami ng mga grid na kinakailangan para sa LES pagkalkula ng isang solong ulo ng kamera ng pagsisilaw lamang ay nasa daang-milya o bilyones. Ang ganitong mataas na paggamit ng mga yunit ng pagkalkula ay limita ang malawak na paggamit ng LES sa simulasyon ng kamera ng pagsisilaw.

Ang pagsisimula ng mga model at paraan ng mataas na katuturang pagkalkula batay sa mga framework ng Very Large Eddy Simulation (VLES) at Hybrid RANS-LES Method ay isang mahalagang trend sa numerikal na simulasyon. Ang VLES method na inilapat ni Han at iba pa ay naglutas ng problema ng mababang epektibidad sa pagkuha ng resulta na dulot ng pag-i-filter ng grid scale at pag-solve ng pagtugma ng turbulensiya scale sa tradisyonal na LES, at naitatag ang kumpiyansang pag-modelo sa pagitan ng mga multi-scale na characteristics ng turbulensiya, transiente na pag-unlad ng characteristics, at grid resolution. Sa pamamagitan ng VLES, adaptibong pinapabago ang proporsyon sa pagitan ng solusyon ng turbulensiya at pag-modelo batay sa talaksan ng real-time na pag-unlad ng vortex structure, siginifikanteng pinaikli ang gastos sa pagkuha ng resulta habang tinuturing ang katuturan ng pagkalkula.

Gayunpaman, kumpara sa tradisyonal na LES, ang teorya at katangian ng VLES ay hindi pa madalas na pinag-aralan at ginamit. Ang artikulong ito ay sistematikong nagpapakita ng teorya ng pag-model ng VLES at ng mga epekto ng kanyang pamamahagi sa iba't ibang pisikal na sitwasyon na may kinalaman sa combustion chambers, na sumusubok ng malawak na pamamahagi ng VLES sa larangan ng simulasyon ng combustion chamber ng eroplano.

Large Eddy Simulation Method

Inipakita sa Figure 1 ang impluwensya ng mga pamamaraan ng simulasyon ng turbulensiya sa paggamit ng mga resource ng kompyuter at sa mga modelo. Ang mga paraan ng RANS, LES, at VLES ay lahat ayumang naiuunlad sa pamamagitan ng pag-modelo ng turbulensiya. Dapat tandaan na si Pope ang unang nagbigay ng malinaw na definisyon ng VLES, na tumutukoy sa 'computational grid scale ay sobrang makapal kaya ang direkta namong inilutas na enerhiya ng turbulensiya ay mas mababa sa 80% ng kabuuang enerhiya ng turbulensiya.' Sa parehong panahon, ang ibinigay ni Pope [6] bilang kahulugan ng LES ay 'ang computational grid ay napakadikit kaya ang direkta namong inilutas na enerhiya ng turbulensiya ay mas mataas sa 80% ng kabuuang enerhiya ng turbulensiya.' Gayunpaman, dapat tandaan na ang ipinakilala sa artikulong ito na VLES ay isang bagong pamamaraan ng pagkuha ng datos na binago at inilimbag batay sa dating paraan. Bagaman ang mga pangalan ay pareho, ang bagong paraan ng VLES ay esesyal naiba mula sa VLES na tinukoy ni Pope. Maaaring makita mula sa larawan, ang tradisyonal na mga mode ng turbulensiya ay RANS, URANS, hybrid RANS/LES, LES, at DNS ayon sa presisyon ng pagsukat. Sa ilalim ng bagong framework ng modelo, ang mga mode ng turbulensiya ay hinati sa RANS, VLES, at DNS ayon sa presisyon ng pagsukat. Ito'y nangangahulugan na ang paraan ng VLES ay nag-uuna ng maramihang tradisyonal na mga mode ng turbulensiya, at ang mga modelong iba't iba ay adaptibong sumusunod at maaaring mag-convert nang maayos batay sa lokal na katangian sa tunay na pagsukat.

Pagpapalipad ng tipikal na mga proseso ng pisikal sa loob ng combustion chamber

Malaking Eddy Simulation ng Malakas na Swirling Flow

Ang combustion chamber ng isang eroplano engine ay madalas na gumagamit ng mga anyo ng organisasyon ng field tulad ng multi-stage swirl at malakas na swirl. Ang swirl flow ay ang pinakabasik na anyo ng pamumuhunan sa combustion chamber. Dahil ang swirl ay dominant sa parehong direksyon ng pamumuhunan at tanggente, mas mabibigat ang turbulent na pagpupulse ng swirl kaysa sa tradisyonal na pipe flow, channel flow at jet flow. Kaya't, ang numerikal na pagpapalipad ng swirl ay nagdadala ng malaking hamon sa pamamaraan ng pagpapalipad ng turbulensiya. Gumamit si Xia at iba pa ng VLES method upang magkalkula ng klasikong halimbawa ng malakas na swirling flow sa tube; ginawa ni Dellenback at iba pa [14] ang mga eksperimento sa field sa halimbawang ito at may detalyadong datos ng eksperimento. Ang Reynolds number ng pamumuhunan sa binibilang na halimbawa ay 1.0 × 105 (batay sa diyametro ng bilog na tube) at ang swirl number ay 1.23. Dalawang set ng strukturadong mga grid ang ginagamit sa pagsukat. Ang kabuuan ng mga sparse grid (M1) ay tungkol sa 900,000 at ang kabuuan ng mga encrypted grid (M2) ay tungkol sa 5.1 milyon. Ang mga resulta ng statistical moment na nakuhang mula sa pagsukat ay patuloy na pinaparehas sa mga eksperimental na resulta upang suriin ang katuturan ng pagsukat ng VLES method.

Ang paghahambing ng mga resulta ng pagkalkula ng iba't ibang paraan at ang mga eksperimental na resulta ng pagsasanay ng distribusyon ng paligid na pang-promedio at bilis ng pagpupulso sa iba't ibang posisyon ng pababa sa ilalim ng malakas na pagigilig ay ipinapakita sa Figura 4. Sa larawan, ang horizontal at vertical coordinates ay walang dimensiyong distansya at walang dimensiyong bilis, kung saan ang D1 ay ang diyametro ng bilog na pipa ng pasukan at ang Uin ay ang promedyong bilis ng pasukan. Maaaring makita mula sa larawan na ang patag ng pamumuhunan ay nagpapakita ng tipikal na Rankin-tulad na kompyund na ulo na anumang pababagong lumilipat papuntang isang solong vortex ng katiging katawan. Pag-uugnay ang mga resulta ng pagkalkula at eksperimento, maaaring matukoy na ang VLES method ay may mataas na katiyakan sa paghula ng bilis ng paligid ng malakas na pagigilig, na nasaayos ang paghula ng mga eksperimental na sukatan. Ang tradisyonal na RANS method ay may napakalaking pagkakaiba sa pagkalkula ng pagigilig na pamumuhunan at hindi maipapahayag ang espasyal na pag-unlad ng patag ng pagigilig at pagpupulso ng turbulensiya. Sa halip, ang VLES method ay may napakataas na katiyakan sa paghula ng patag ng promedyo ng bilis, patag ng pagpupulso ng bilis at espasyal na pag-unlad sa ilalim ng kompleks na malakas na pagigilig, at maaaring patuloy na magbigay ng mataas na katiyakan ng pagkalkula kahit sa isang relatibong masinsing resolusyon ng grid. Para sa paghula ng promedyo ng bilis ng paligid, ang mga resulta ng pagkalkula ng VLES method ay halos konsistente sa dalawang grupo ng masinsing at masinsing resolusyon ng grid.

Malaking Eddy Simulation ng Turbulent Combustion

Upang pag-aralan ang kaya nang mangyari ng VLES na paraan sa pagsusuri ng mga problema ng turbulent combustion [15-16], inilapat ang isang modelo ng turbulent combustion na batay sa VLES method na nakakabit sa flamelet generated manifolds (FGM). Ang pangunahing ideya ay ang ipinapalagay na mayroong lokal na isang-dimensional na laminar flame structure ang turbulent flame, at ang ibabaw ng turbulent flame ay ang ensemble average ng isang serye ng laminar flame surfaces. Kaya nito, maaaring imapekta ang mataas na dimensional na component space sa isang mababang dimensional na paternong pagsisikap na binubuo ng ilang karaniwang mga variable (mixture fraction, reaction progress variable, etc.). Sa kondisyon na kinokonsidera ang detalyadong reaksyon mechanism, tinataas ang bilang ng transport equations na kinakailangang sulisin, kaya napapababa ang komputasyonal na gastos.

Ang tiyak na proseso ng pagsasagawa ay ito: gumawa ng talahanayan ng datos ng laminar ng FGM batay sa fraction ng pagmiksa at mga variable ng progreso ng reaksyon, konsidera ang interaksyon sa pagitan ng turbulent combustion sa pamamagitan ng pag-aasumpsiyon ng paraan ng probability density function upang mag-integrate ng talahanayan ng laminar data, at kung gayon ay makakuha ng talahanayan ng datos ng turbulensya. Sa numerikal na pagkuha, inisolve ang mga ekwasyon ng transportasyon ng fraction ng pagmiksa, mga variable ng progreso ng reaksyon, at ang katumbas na variance, at nakukuha ang impormasyon ng field ng pagkabubo sa pamamagitan ng paghahanap sa talahanayan ng datos ng turbulensya.

Ginamit ang modelo ng turbulenteng pagsisilaw na batay sa VLES at FGM upang ipagawa ang numerikong pagkuha ng datos sa methane/air turbulent jet flame (Flame D) na tinukoy ng Sandia laboratory sa Estados Unidos, at ginawa ang kuantitatibong paghahambing sa mga eksperimental na datos. Ang materyales ng combustible ng halimbawa ng Sandia Flame D (Reynolds number ay 22400) ay isang buong mikstura ng methane at hangin na may proporsyon ng 1:3, ang bilis ng pagnyusd ng fuel ay halos 49.9 m/s, at ang bilis ng pagkabulaklak ay halos 11.4 m/s. Ang duty flame ay isang mikstura ng sinunog na methane at hangin, at ang materyales ng pagkabulaklak ay pure hangin. Ginagamit sa pagkuha ng datos ang strukturadong grid, at ang kabuuang bilang ng mga grid ay halos 1.9 milyon.

Ipinapakita sa Figura 5 ang pagkalat ng pangkalahatang bahagi ng masa ng mga iba't ibang komponente sa hilagang-aksyal. Ang horizontal at vertical na koordinadong ipinapakita sa larawan ay walang-dimensyon na distansya (D2 ay ang diyametro ng inlet jet tube) at walang-dimensyon na bahagi ng masa, sa kabila ng kung ano man. Maaaring makita mula sa larawan na ang paghahambing ng mga pangunahing komponenteng kinikilala ng paraan ng VLES sa proseso ng pagsusunog ay umuugnay nang husto sa mga eksperimental na resulta. Ipinapakita sa Figura 6 ang trend ng pagkalat ng temperatura sa iba't ibang posisyon sa ilalim ng espasyong pagkakasunod-sunod ng pagmamasid. Maaaring makita mula sa larawan na ang trend ng pagkalat na inihanda ng pamamaraan ng VLES ay pangkalahatang konsistente sa mga eksperimental na resulta, at lamang ang pinakamataas na balor ng temperatura ay kaunting mas mataas kaysa sa eksperimental na balor. Ipinapakita sa Figura 7 ang pagkalat ng agampan, temperatura, at resolusyong kontrol na punsiyon na hinati ng VLES, kung saan ang linyang solid ay kinukuha bilang Zst=0.351. Maaaring makita mula sa larawan na ang rehiyon ng core jet ay nagpapakita ng malakas na turbulenteng pagpulso, at habang lumalago ang patlang ng ilalim ng ilalim, ang sukat ng struktura ng bulok ay paulit-ulitumigtaas. Gayunpaman, mula sa Figura 7 (b) at (c), sa karamihan sa mga rehiyon ng kimikal na reaksyon, ang resolusyong kontrol na punsiyon ay nasa pagitan ng 0 at 1, na nagpapahayag na ang lokal na grid na resolusyon ay makakakuha ng malaking-eskalang turbulensiya at lamang simulan ang maliit na eskala ng turbulensiya sa pamamagitan ng modelo. Sa oras na ito, ang VLES ay gumaganap bilang isang halos malaking eddy simulation mode. Sa rehiyon ng jet shear layer at sa labas na hangganan ng apoy sa ilalim ng ilalim, ang resolusyong kontrol na punsiyon ay malapit sa 1, na nagpapahayag na ang tinatapos na filter scale ng computational grid ay mas malaki kaysa sa lokal na eskala ng turbulensiya. Sa oras na ito, ang VLES ay gumaganap bilang isang unsteady Reynolds average mode. Sa karagdagan, maaaring makita na ang paraan ng VLES ay maaaring mag-realize ng pagbabago ng maraming solusyon ng turbulensiya ayon sa tunay na panahon na katangian ng pag-unlad ng struktura ng bulok, at maaaring makipag-ugnayan nang husto sa hindi-pantay na proseso ng pagsusunog sa malubhang apoy.

Simulasyon ng malaking eddy para sa buong proseso ng atomization

Karamihan sa fuel na ginagamit sa combustion chamber ng isang aircraft engine ay liquid fuel. Ang liquid fuel ay pumapasok sa combustion chamber at dumarating sa mga proseso ng primary at secondary atomization. Mayroong maraming hamon sa simulasyon ng buong proseso ng atomization ng liquid fuel, kabilang ang pagkuha ng gas-liquid two-phase topological interface configuration, pagdeform at pagbreak ng liquid column, pagbreakdown evolution ng liquid bands at liquid filaments sa droplets, at interaksyon ng turbulent flow at droplets. Ang Huang Ziwei [19] ay nagdevelop ng isang model ng simulasyon para sa buong proseso ng atomization base sa VLES method na nakakonekta sa VOFDPM hybrid atomization calculation method, narealize ang full-process numerical simulation ng fuel atomization mula sa continuous liquid hanggang discrete droplets.

Ginamit ang isang bagong nilikha na modelo ng pag-simulasyon ng proseso ng atomization upang ipagawa ang mataas na katumpakan na pagsukat bilangiko ng tradisyonal na proseso ng atomization ng ulo ng likido, at ginawa ang detalyadong pagsusuri nito kasama ang mga resulta ng eksperimento sa pampublikong literatura [20] at ang mga resulta ng pagkuha ng sukat ng malalaking eddy simulation [21]. Sa halimbawa ng pagsukat, ang fase ng hangin ay hangin na may bilis na 77.89 at 110.0 m/s, samantalang ang fase ng likido ay tubig na may bilis na 8.6 m/s. Ang mga sumeserye na Weber number ay 100 at 200, na bawat isa. Upang mas maayos simulan ang proseso ng pangalawang pagputol, ginagamit ng modelo ng pagputol ang Kelvin-Helmholtz at Rayleigh-Taylor (KHRT) model.

Ang buong proseso ng atomisasyon na inihulaan ng VLES sa kondisyon ng Weber number 100 ay ipinapakita sa Figura 8. Gaya ng makikita mula sa larawan, isang mababaw na sheet ng likido column ay nabubuo sa simula, at pagkatapos ay nagbubreak ang likido column sa mga liquid bands at liquid filaments, at nagbubreak sa mga droplets dahil sa pwersa ng aerodinamiko, at ang mga droplets ay patuloy na nabubuo bilang mas maliit na droplets sa pamamagitan ng ikalawang pagbreakup. Ang stream velocity at spanwise vorticity distribution na binuhos ng VLES sa kondisyon ng Weber number 100 ay ipinapakita sa Figura 9. Gaya ng makikita mula sa larawan, mayroong tipikal na lugar ng low-speed recirculation sa leeward side ng likido column. Maaaring matukoy mula sa instantaneous vorticity distribution na ang leeward side ng likido column ay nagpapakita ng malakas na vortex structure, at ang malakas na galaw ng turbulento sa low-speed recirculation zone ay nag-uugnay sa pagbreak ng likido column sheet at sa pagsisimula ng mga droplets.

Ang ratio ng diametro ng initial jet sa pinakamaliit na dimensional na pamumuhunan ng likido jet nang ang kolmun ng likido ay umuubos sa ilalim ng magkakaibang Weber numbers ay ipinapakita sa Figure 10. Sa larawan, ang di ay ang pinakamaliit na dimensional na pamumuhunan ng likido jet nang ang kolmun ng likido ay umuubos, at ang D3 ay ang diametro ng initial liquid jet. Maaaring makita mula sa larawan na ang mga resulta ng pagkuha ng VLES ay magandang sumasang-ayon sa mga eksperimental na resulta, na mas mabuti kaysa sa mga resulta ng pagkuha ng malaking eddy simulation sa literatura [21].

Pagbagsak ng Pagsusunog sa Malaking Eddy Simulation

Upang tugunan ang mga kinakailangang mababang emisyon, karaniwang disenyo ang mga kamera ng pagbubunsong sibil na eroplano gamit ang premixed o bahagyang premixed lean combustion. Gayunpaman, ang lean premixed combustion ay may mahinang kasarian at madaling makipot sa thermokustikong nakasang-ayos na mode ng pagbubunsong osilasyon, na nagiging sanhi ng kawalan ng kasarian sa pagbubunsong. Ang kawalan ng kasarian sa pagbubunsong ay sobrang destruktibo at maaaring dumaan sa mga problema tulad ng flashback at solid deformation, na ito ang isang napakalaki na problema na kinakaharap ng disenyo ng kamerang pagbubunsong.

Ang pagsukat ng pamamaraan ng pagkakabu-bos na numerikal ay maaaring ibahagi sa dalawang kategorya: paraan ng paghiwa at direktang paraan ng pag-uugnay. Ang paraan ng paghahalintulad ng kaguluhan sa pagkakabu-bos ay hiwa ang solusyon ng hindi makapagpatuloy na pagkakabu-bos at akustiko. Kinakailangan ng hindi makapagpatuloy na pagkakabu-bos maraming halimbawa ng pagsukat na numerikal upang magtayo ng handa na paglalarawan ng flame. Kung ginagamit ang pamamaraan ng malaking eddy simulation, sobrang malaki ang kinakailangang yunit ng pagsasakauna. Ang direktang paraan ng pag-uugnay ay batay sa compressible solution method, at direkta nang nakukuha ang resulta ng kaguluhan sa pagkakabu-bos sa pamamagitan ng mataas na katumpakan ng hindi makapagpatuloy na pagsukat, o sa ibang salita, kumpletong ugnayan ang proseso ng pagsukat ng hindi makapagpatuloy na pagkakabu-bos at akustiko sa mga ibinigay na kondisyon ng trabaho sa loob ng parehong framework ng pagsukat.

Sa pag-aaral ng pagsisiwalat ng kawing na simulasyon sa pagkabuhay ng pagkakaburno, ni-Huang et al. [27] ay nai-develop ang isang modelo ng pagkalkula ng kawing na pagkakaburno batay sa paraan ng VLES na nakakabit sa makapal na paraan ng pagsusuri ng flame, at natugunan ang tunay na paghula ng di-madalas na proseso ng pagkakaburno sa ilalim ng akustikong ekspisyon. Ang halimbawa ng pagkalkula ay isang matatag na butas na may ganap na pre-mixed na flame ng etilen/ Hanging ng Cambridge University, na may katumbas na ratio ng 0.55 at Reynolds number na tungkol sa 17000. Ang paghahambing sa mga resulta ng pagkalkula ng VLES at mga eksperimental na resulta ng di-madalas na dinamika ng flame sa ilalim ng akustikong ekspisyon ay ipinapakita sa Figura 12. Makikita mula sa larawan na habang nagaganap ang proseso ng pagsisigla sa entrada, lumiliko ang flame sa loob at labas ng shear layers at umuunlad sa isang counter-rotating vortex pair. Sa prosesong ito, patuloy na umuunlad ang pag-unlad ng profile ng flame na hugis-sampaguita kasama ang pagbabago ng phase angle. Mabuti ang mga resulta ng pagkalkula ng VLES na bumubuhay muli sa mga karakteristikang obserbado sa eksperimento. Ang paghahambing sa amplitud at phase difference ng reaksyon ng rate ng pagpapalaya ng init sa ilalim ng 160 Hz akustikong ekspisyon na kinuha sa pamamagitan ng iba't ibang mga paraan ng pagkalkula at eksperimental na pagsukat ay ipinapakita sa Figura 13. Sa larawan, Q' at Q ͂ ang mga pulsatibong paglabas ng init at pangkalahatang paglabas ng init ng pagsisilbi, na ang A ay ang amplitud ng sinusoidal na akustikong pagsusuri, at ang ordinate ng Figura 13 (b) ay ang kakaibang fase sa pagitan ng transiyenteng signal ng paglalabas ng init ng pagsisilbi sa ilalim ng akustikong pagsusuri at ng signal ng pagsusuri ng bilis ng pasukan. Kabilang sa figure, ang katumpakan ng prediksyon ng paraan ng VLES ay maaaring kumumpara sa katumpakan ng malaking eddy simulation [28], at pareho ay magandang sumasang-ayon sa eksperimental na halaga. Bagaman ang paraan ng unsteady RANS ay nagpaprediksyon ng trend ng hindi linyang tugon, ang kinalkulang kuantitatibong resulta ay malayo sa eksperimental na halaga. Sa mga resulta ng kakaibang fase (Figura 13 (b)), ang trend ng kakaibang fase na inipredict ng pamamaraan ng VLES kasama ang amplitud ng pagdudyok ay pangunahing sumasang-ayon sa eksperimental na resulta, samantalang ang mga resulta ng malaking eddy simulation ay hindi maaaring mabuti ipinaprediksyon ang taas na trend.