Beräkningsmetoden för att simulera och visualisera virvelzoner (VZ) i avgasflöden är komplex och omfattar flera steg. För att förstå hur dessa zoner påverkar flödet vid inloppet till en rektangulär avgaskanal, måste man först skapa en detaljerad beräkningsmatris och sedan använda den för att lösa de algebraiska ekvationerna som beskriver flödet.

Först genereras en matris GpqG_{pq}, där p=1,2,,N1p = 1, 2, …, N-1 och q=1,2,,Nq = 1, 2, …, N, som beskriver effekten av de bifogade virvlarna (fyrkantiga virvelramar och raka rektangulära virvlar med enhetsintensitet) på beräkningspunkterna. Detta görs för att ta hänsyn till hur dessa virvlar påverkar flödet genom kanalen. Därefter skapas en matris GN1qG_{N1q}, där q=1,2,,Nq = 1, 2, …, N, som beskriver effekten av enhetsintensitetsvirvlar på de beräknade punkterna längs den aktiva tvärsnittet i avgassystemet. Dessa matriser används för att beräkna flödeshastigheten vid varje punkt på tvärsnittet av kanalen och upprätta en grundläggande förståelse för hur virvlarna sprider sig.

Efter att ha skapat dessa matriser, påbörjas den iterativa processen för att bestämma den fria virvelströmsytan, det vill säga den yta som genereras från den skarpa kanten på flänsen. I denna process definieras de initiala värdena för antalet hörnpunkter i de fria virvelpolylinjerna som γnew=0\gamma_{\text{new}} = 0 och γold=0\gamma_{\text{old}} = 0, samt för cirkulationen av de fria krökta ramverken. Den iterativa metoden fortskrider genom flera nivåer av yttre och inre loopar, där man successivt beräknar cirkulationen för de fria virvlarna och uppdaterar dessa värden tills precisionen γnewγold\gamma_{\text{new}} - \gamma_{\text{old}} är tillräckligt låg, vanligtvis under ε=105\varepsilon = 10^{ -5}.

Denna process leder till lösningen av ett system av simultana linjära algebraiska ekvationer, som används för att beräkna de okända intensiteterna av virvelementen från randvillkoren för den normala hastighetskomponenten. Efter att ha löst dessa ekvationer erhålls de nya värdena för cirkulationen av de fria virvelpolylinjerna och ramarna. Genom att iterera denna process förfinas beräkningarna av virvelströmmens yta och flödesfältet.

När dessa beräkningar är slutförda kan man definiera hastighetsfältet och plotta strömlinjerna som beskriver flödet. Detta görs genom att använda en formel som tar hänsyn till både de fria virvlarna och randvillkoren för flödet vid inloppet till kanalen. På detta sätt kan man visuellt och kvantitativt uppskatta hur virvelzonerna påverkar flödet i kanalen och definiera deras karaktäristiska dimensioner.

För att ytterligare förstå hur dessa virvelzoner påverkar flödet används experimentella metoder för att visualisera gränserna för virvelzonerna vid inloppet till den kvadratiska avgaskanalen. En experimentell anordning gör det möjligt att observera och mäta dessa zoner i ett brett intervall av Reynolds-tal (från 6 · 10³ till 9 · 10⁴), vilket täcker både början på den turbulenta regeln och det utvecklade turbulenta flödet. Här används en ultraljudsfuktare för att generera ånga som visualiserar flödet vid separationen från den skarpa kanten på kanalen. Dessa experimentella observationer bekräftar teoretiska beräkningar och bidrar till en mer komplett förståelse av virvelzonernas form och påverkan på flödet.

Den datorbaserade programvara som utvecklats för att beräkna luftflödets hastighet och gränserna för strömskillnadens yta vid ingången till den rektangulära avgaskanalen ger en värdefull metod för att både visualisera och kvantifiera flödesdynamiken. Programmet gör det möjligt att sätta upp randvillkor för den normala hastighetskomponenten, skapa strömlinjer och bestämma hastigheten i flödesområdet. Det ger också möjlighet att konstruera fria strömsystem och bestämma deras karakteristiska dimensioner.

Det är viktigt att förstå att även om dessa beräkningar och experimentella observationer ger en detaljerad bild av virvelzonernas natur, är det en dynamisk process som kräver kontinuerlig förbättring och uppdatering av både matematiska modeller och experimentella metoder. Den noggranna simuleringen av flödet genom avancerade iterativa metoder gör det möjligt att optimera konstruktioner av avgaskanaler och ventilationssystem där virvelzonernas påverkan är betydande.

Hur påverkar luftflödet öppningar i ventilationssystem?

I en studie av Logachev et al. (2010c) användes CMM (Conformal Mapping Method) för att bestämma konturen för virvelzoner (VZ) vid ingången av en ventilationshuv i förhållande till avståndet från utsläppsöppningen till väggen. Avståndet varierade mellan S/b = 0 och S/b = ∞. Resultaten visade hur strömlinjer och utrustningslinjer utvecklas beroende på dessa avstånd (se figur 1.16). Forskarna utförde också en jämförelse mellan de erhållna hastighetsfälten i den initiala delen av utsläppsröret och experimentella data. För att beräkna LDC (Luftflödets dragkraftskompression) föreslogs Borda-formeln som en funktion av flödets kompression inne i röret, δ: ζ = (1/δ – 1)². För det fria inloppsröret (S/b = ∞), där δ är 0,5, var LDC ζ = 1. Vid en slitsöppning i en oändlig vägg (S/b = 0), där δ = π/(π + 2), var LDC ζ = 0,405285.

En rad studier har visat hur DVM (Digital Volume Method) och konforma avbildningar har tillämpats framgångsrikt på problem som rör luftflöden mot utsläppsöppningar med mer komplexa geometriska former, som till exempel ett utsläppshuvud bestående av en kvadratisk anslutning fritt placerad i rymden (Averkova et al., 2013c), eller utsläppshuvuden försedda med ogenomträngliga skärmar (Logachev et al., 2010a, 2010b; Averkova et al., 2013a, 2013b), eller utsläppshuvuden installerade inuti ett kåpa (Averkova et al., 2007, 2010).

Vid undersökning av sidouppgångar i ventilationssystem liknar luftflödet det som ses vid konfluens-tees, där huvudkanalens flöde blandas med flödet från den sidogående öppningen, vilket orsakar chockförluster. En viktig skillnad är att flödet vid en tee kommer från en lateral kanal, medan flödet vid en sidouppgång kommer från ett obegränsat utrymme och gör en 180° sväng. I Ryssland används ofta Designer's Handbook (Barkalov et al., 1992) för aerodynamiska beräkningar av ventilationssystem, som huvudsakligen består av anpassade data från Idel’chik (1992). För specifika fall som sidouppgångar är det dock inte helt klart varifrån dragvärdena hämtas, och dessa uppgifter återfinns inte i Idel’chik (1992).

För en sidouppgång i ventilationssystem kan dragförluster (LDC) beräknas för en rad parametrar, inklusive förhållandet mellan öppningens area och kanalens tvärsnittsarea, där värden mellan 0,1 och 0,6 är vanliga. Förhållandet mellan flödeshastigheter i öppningen och huvudkanalen varierar också, där flödeshastigheterna i öppningen kan vara mellan 0,1 och 0,5 (se figur 1.17). Forskningsresultat från Saito och Ikohagi (1994) visade att LDC för flöden genom sådana öppningar kan variera från 0 vid noll flödeshastighet till omkring 1,05 för ett öppningsflöde i intervallet 0,349 till 0,785, vilket är i linje med dragvärdena för tees enligt Maeda et al. (1960).

I studier av Maeda et al. (1960) visades också en god överensstämmelse mellan LDC-värden för flöden genom tees och för sidouppgångar. Det är värt att notera att för LDC i tees, enligt Idel’chik (1992), är detta värde oberoende av dimensionerna på sidogrenen, vilket strider mot Maeda et al. (1960).

Vidare pekar en numerisk studie av Maklakov et al. (2020) på att även om flödet genom en sidouppgång påminner om flödet genom en frontal öppning, är processen mer komplex. Hanzhonkov och Davydenko (1959) publicerade experimentella resultat på drag i runda rör, där LDC-relationen för slutliga sidouppgångar bestämdes som en funktion av öppningens area (F) i förhållande till kanalens area (F), samt avståndet från öppningen till dödänden (s). Vid små öppningsdimensioner (F/F ≈ 0,9) förändras flödets natur drastiskt, och flödet delar sig för att passera runt flödet som går in i kanalen, vilket medför en kvalitative förändring i flödets karaktär.

För att förstå dynamiken i dessa flöden är det viktigt att förstå hur storleken på öppningen påverkar strömningens beteende. Små öppningar nära dödänden kan orsaka stora förluster i flödet, medan större öppningar kan ge upphov till olika typer av flödesmönster som påverkar ventilationens effektivitet. För att optimera designen av ventilationssystem är det avgörande att noggrant välja rätt öppningsstorlek och placering i relation till kanalens utformning för att minimera dessa förluster och maximera effektiviteten i systemet.