Numeriska lösningar med hjälp av Detached Eddy Simulation (DES) och Improved Delayed Detached Eddy Simulation (IDDES) har visat sig vara begränsade när det gäller att korrekt förutsäga aerodynamiska koefficienter som lyftkraft och moment vid höga anfallsvinklar, där dessa koefficienter uppvisar icke-linjära beteenden. Specifikt fångar dessa metoder inte tillfredsställande den tvärgående (spanwise) komponenten i flödet över de stallade delarna av en svept vinge, vilket leder till felaktiga förutsägelser av tryckfördelningen på vingens ovansida. Detta är särskilt tydligt vid flow separation där flödet delvis lämnar vingytan och visar komplexa, tredimensionella egenskaper.

I motsats till detta har konventionella Reynolds-Averaged Navier-Stokes (RANS) metoder, trots sin förenklade turbulensmodellering, visat sig vara överlägsna i att kvantitativt fånga denna specifika dynamik. Detta antyder att de hybridmodeller som DES och IDDES bygger på fortfarande kräver vidareutveckling för att hantera de komplexa flödesfenomen som uppstår vid stall på svepta isbelagda vingar.

En kritisk aspekt är övergången mellan RANS- och LES-regionerna i dessa hybridmodeller. I områden där flödet är delvis separerat och där den spanwise flödeskomponenten är betydande, är turbulensmixningen komplex och inte tillräckligt väl fångad av nuvarande modeller. Det behövs därför förbättrade övergångsmodeller som kan hantera denna tre-dimensionella aerodynamik mer precist.

Utmaningen för numeriska simuleringar av isbelagda svepta vingar ligger i flödets komplexitet och den stora variationen i både geometrin (isackretionernas komplexa form) och de aerodynamiska förhållandena. Experimentella data från NASA:s Glenn Research Center har visat att även vid relativt små anfallsvinklar, redan från cirka 6 grader, uppstår delvis flödesseparation längs vingens längd, vilket skapar höggradigt ostadiga och tredimensionella flödesmönster. Dessa egenskaper gör det svårt för turbulensmodeller att ge tillförlitliga förutsägelser.

Det är viktigt att förstå att den aerodynamiska påverkan av isbildning på svepta vingar är ett nationellt säkerhetsproblem som är dåligt kartlagt, delvis på grund av brist på detaljerade experimentella data med realistiska isformer på moderna vinggeometrier. Genom att använda laserskanning och digitalisering av experimentellt ackumulerade isformer har forskare kunnat skapa noggranna CAD-modeller som kan användas för vidare numeriska och experimentella studier.

Det är dessutom centralt att CFD (Computational Fluid Dynamics) inte bara är ett verktyg för att simulera flöden, utan också ett potentiellt kostnadsbesparande hjälpmedel i design- och certifieringsprocessen för flygplan. Med förbättrade numeriska metoder kan problem som uppstår i tidiga designfaser identifieras och åtgärdas, vilket minskar behovet av omfattande och kostsamma vindkanal- och flygtester.

Den nuvarande forskningen och experimentella arbetet kring DES och IDDES visar dock att dessa metoder ännu inte är tillräckligt robusta för att på ett tillförlitligt sätt simulera de mycket komplexa och ostadiga flöden som uppstår vid stall på svepta isbelagda vingar. Detta beror till stor del på svårigheterna i att korrekt hantera övergången mellan RANS- och LES-regionerna i närvaro av starka spanwise flödeskomponenter och stora flödesseparationer.

Därför är det nödvändigt att vidareutveckla övergångsmodeller och turbulensmodeller som kan integrera turbulent mixing bättre och därigenom fånga den tredimensionella karaktären av stallflödet. Detta kräver inte bara avancerade numeriska tekniker utan också ökad förståelse för flödets fysik vid komplexa aerodynamiska tillstånd.

Vikten av dessa insikter sträcker sig bortom enbart akademiska intressen. Den direkta kopplingen till flygsäkerhet och de höga kostnaderna för iscertifiering gör det nödvändigt att finna pålitliga och effektiva numeriska metoder för att kunna förutsäga och kontrollera effekterna av isbildning på kommersiella flygplansvingar. Den pågående forskningen från internationella konsortier visar vägen framåt, där experimentella data och avancerade CFD-verktyg kombineras för att skapa en mer heltäckande förståelse av dessa komplexa fenomen.

Det är också avgörande att ha en realistisk bild av modellernas begränsningar och att kontinuerligt verifiera och validera numeriska resultat mot experimentella data. Detta för att säkerställa att modellerna inte bara fungerar under idealiserade förhållanden utan även för de mest krävande och dynamiska flygförhållandena som uppstår vid isbildning och stall på svepta vingar.

Hur beräknas isbildning och rotorbladsdynamik vid isbildning på rotorfordon?

Beräkningen av isbildning på rotorfordon utgår från flera sammankopplade moduler som tar hänsyn till aerodynamik, termodynamik och isens tillväxt över tid. Termodynamiska moduler är centrala för att bestämma isens ansamling på ytan, där Messinger-modellen, ursprungligen utvecklad 1953, är en grundläggande princip. Modellen baseras på termodynamikens första lag och balanserar energin som går in och ut ur ett kontrollvolym, där värmeförlusten från ytan används för att producera den latenta värmen som krävs för isbildning. Denna energibalans kopplas till massbalansen för att beräkna isackumuleringen.

Messinger-modellen existerar i flera varianter, där bland annat den utökade versionen integrerar Stefan-ekvationen, vilken beskriver fasövergångar mellan vatten och is med hjälp av fyra kopplade ekvationer. En annan variant, shallow water icing model (SWIM), som används i FENSAP-ICE, löser ett ofullständigt system av partiella differentialekvationer genom att anta olika ytförhållanden som rimfrost, blandad is-vatten eller bara vatten för att iterativt bestämma yttemperatur och ismassa. Termodynamiska modeller arbetar ofta med begreppet frysningsfraktion, som beskriver förhållandet mellan is och inkommande massa och itereras tills lösningen är konsistent med antagna ytförhållanden.

Isens tillväxtberäkning utförs genom att multiplicera isackumulationshastigheten med exponeringstiden, varefter ytan modifieras för att inkludera den nya istjockleken. Denna process kräver ofta flera omgenerationer av nätet för att exakt modellera isens form och påverkan på luftflödet.

Rotorns aerodynamik analyseras separat från flygkroppen, eftersom rotorn rör sig mycket oregelbundet och påverkar isbildningsfenomenet på ett komplext sätt. Rotorns rörelse och isbildning är båda intrinsikalt dynamiska och gör simulering av rotorisk isbildning till en komplex multidisciplinär utmaning. Metoder för rotorsimulering har utvecklats för att beskriva flödesfältet runt rotorbladen och integrera aerodynamiska data med isbildningsmodeller.

Analytiska modeller såsom momentumteorin, ursprungligen utvecklad för marina propellrar, betraktar rotorn som en tunn skiva och ger en första approximation av rotorprestanda genom att beskriva luftflödet och tryckskillnader som genererar dragkraft. Denna teori är dock begränsad eftersom den inte tar hänsyn till energiförluster som inducerad effekt.

Blade Element Theory (BET) delar upp rotorn i flera tvådimensionella tvärsnitt, där lyft- och dragkrafter beräknas för varje delsektion. BET möjliggör detaljerad analys av bladegenskaper som bladvridning och profil, men har brister i att hantera tredimensionella effekter, ömsesidiga interaktioner mellan bladsektioner och komplexa flödesfenomen som spetsvirvlar.

Blade Element Momentum Theory (BEMT) kombinerar BET med momentumteori för att nyttja fördelarna hos båda modellerna och hantera både lokala och globala effekter i flödet genom rotorn. Denna hybridmodell delar upp rotorbladet i infinitesimala ringar där massflödet och den genererade dragkraften beräknas lokalt med hjälp av momenteteorin. BEMT är en effektiv metod för att analysera rotorns totala prestanda under hover- och axiella rörelser.

Vid beräkningarna är det viktigt att ta hänsyn till att både isbildning och rotorens rörelse är starkt dynamiska processer som kräver iterativa och ofta multidisciplinära lösningar för att uppnå realistiska och tillförlitliga simuleringar. Interaktionen mellan de termodynamiska processerna och den komplexa aerodynamiken runt roterande blad ställer höga krav på både modellens fysikaliska trovärdighet och beräkningsresurser.

Viktigt att förstå är att de modeller som används alltid bygger på antaganden och förenklingar, vilket innebär att de numeriska resultaten måste tolkas med förståelse för modellernas begränsningar. Dessutom kräver noggranna simuleringar ofta anpassningar av modellen för att hantera specifika flödes- och isbildningsförhållanden samt dynamiska effekter kopplade till rotorns rotation och de förändrade ytförhållandena när isen växer.

Hur påverkar ytegenskaper och superkylda droppars deformering deras frysningsprocess?

I de senaste decennierna har mycket forskning lagts på att förstå och simulera droppars påverkan på olika ytor, särskilt i relation till deras dynamik vid impakt och frysningsprocesser. För att förstå och effektivt hantera dessa fenomen, till exempel för tillämpningar inom flygplansisoleringssystem och andra tekniska lösningar, är det viktigt att både teoretiska och experimentella metoder används. Denna förståelse är också avgörande för att utveckla metoder för att minska frysning eller accelerera smältning av droppar på specifika ytor.

För att undersöka dropparnas impakt och frysningsbeteende på olika hydrophobicitetsnivåer, har numeriska simuleringar blivit ett viktigt verktyg. En av de mest använda metoderna är Lattice Boltzmann Method (LBM), som har visat sig vara effektiv när det gäller att hantera komplexa ytor och gränser. Genom att tillämpa LBM och pseudopotentialteori har forskare kunnat minska kontaktiden mellan dropparna och ytorna med mer än 40 %. Trots detta är dropparnas impaktbeteende fortfarande inte helt förstått, delvis på grund av dess komplexa natur och beroende av faktorer som ytkontaktvinklar och temperaturer.

En särskilt intressant aspekt är påverkan av ytor med olika grader av superhydrofobicitet, det vill säga ytor som har en extremt låg vattnets kontaktvinkel. På sådana ytor tenderar dropparna att "rulla av" utan att stanna, vilket påverkar både impakt och frysning. Emellertid påverkas dessa fenomen ytterligare när ytorna har olika former av utbuktningar – t.ex. triangulära, fyrkantiga eller runda förhöjningar – vilket i sin tur förändrar hur droppen rör sig, stannar eller fryser.

Vid användning av LBM för att undersöka droppars dynamik på plana superhydrofoba ytor kan flera faktorer analyseras i detalj, däribland påverkan av omgivningstemperatur, We (Webertal) och Re (Reynoldstalet). Detta gör det möjligt att få en djupare förståelse för de mekanismer som styr frysning och dynamik under olika förhållanden.

Ett ytterligare intressant område av forskning rör frysningsprocessen för superkylda droppar. När en vätskedroppe är underkyld – det vill säga den är under fryspunkten men har inte frusit än – kan den frysa snabbt om den träffar en yta med en viss egenskap. Frysningens hastighet påverkas av ytkontaktvinkeln, vilket innebär att droppar på ytor med större kontaktvinklar tenderar att frysa senare. I experiment utförda av Huang et al. (2012) har det visat sig att större kontaktvinklar förlänger den totala frystiden för droppen, medan temperatur på ytan och volymen på droppen också har en betydande inverkan.

Ytor med förhöjningar påverkar frysningsprocessen ytterligare, då de kan skapa lokala förhållanden som accelererar frysning eller förändrar hur isen bildas på droppens yta. Det har observerats att isens profil på ytor med förhöjningar ofta får en spetsig form, vilket kan bidra till snabbare isbildning. Men för att bättre förstå dessa processer krävs detaljerade simuleringar som inte bara beaktar den geometriska effekten utan också förändringar i volym och temperatur under frysningsprocessen.

Vidare är användningen av numeriska modeller som kombinerar multiphase LBM och fasmätningstekniker en lovande metod för att noggrant förutsäga frysningsprocessen, särskilt för att ta hänsyn till den volymexpansion som sker när vätskedroppar fryser. Genom att undersöka dessa fenomen kan vi få en mer fullständig bild av hur dropparna deformerar sig och fryser på kalla ytor.

Förutom simuleringarna har det också utförts experiment för att studera droppars frysningsprocesser under olika förhållanden. Frysningen börjar ofta vid en punkt där droppens temperatur är tillräckligt låg, och utvecklingen av frysningszonen kan övervakas genom att studera fysiska gränsytor mellan vätska och is. Här kan den tidiga frysningsfasen påskyndas av lägre yttemperaturer och kan också visas genom en karaktäristisk form av isens struktur.

Det är också viktigt att förstå att när droppar fryser på en yta förändras inte bara droppens form utan också den omgivande temperaturdynamiken och fysiska egenskaper på ytan. Frysningshastigheten varierar beroende på flera faktorer, inklusive droppens storlek, den grad av underkylning, samt ytan där frysningsprocessen sker. En noggrant utformad yta kan därför bidra till att minska frysningshastigheten eller tvärtom accelerera processen, vilket är viktigt att beakta vid design och konstruktion av komponenter som ska exponeras för droppar i kalla miljöer.

Hur påverkar syntetiska jetaktuatorer isbildning och isuppbyggnad på ytor i kalla miljöer?

I studier av syntetiska jetaktuatorer (SJA) och deras effekt på isbildning visar numeriska simuleringar hur dessa aktuatortekniker påverkar tryckfält, hastighetsfält och temperaturdistributioner runt en kileformad yta. Genom att använda en SJA som aktivt manipulerar luftflödet nära ytan kan man skapa luftfickor intill kilens vägg, vilka i sin tur tvingar nedfallande droppar att flöda om och bakåt längs ytan. Denna omdirigering av dropparna minskar isackumuleringen genom att hindra direkt nedslag och ansamling på kilens främre delar.

Tryckfördelningen uppvisar konsekvent ett högtrycksområde vid kilens spets, medan zoner med lägre tryck bildas nedströms från varje SJA, vilket påverkar strömningsmönstret och bidrar till förändrad droppbana. Effekten på isbildning varierar beroende på droppstorleksfördelningen; monodispersa droppar, som är mer homogena och ofta mindre, skjuts undan effektivare än de större dropparna i en Langmuir-D-fördelning. Därigenom reduceras isbildningen mer markant vid monodispers distribution, medan större droppar i Langmuir-D fördelningen har svårare att undvikas av SJA-effekten.

Vid användning av uppvärmda SJAs, där aktuatortemperaturen kan nå upp till 100 °C medan omgivande luft är långt under fryspunkten, uppstår ett dynamiskt samspel mellan värmetillförseln och det omkringliggande kalla flödet. Detta medför att temperaturen vid kilens yta hålls långt över fryspunkten trots kylig fri strömning. SJA-cyklernas olika faser visar hur värmeflödet varierar i tid och rum, där residualvarm luft från sugfasen och intensiv jetstråle under ejektionsfasen tillsammans bidrar till att förhindra isbildning. Resultatet blir att isen nästan helt uteblir, med endast minimala isavlagringar i vissa områden.

Temperaturfältets ojämlika fördelning och de vortikulära strukturer som uppstår vid interaktionerna mellan SJA-flödet och den fria strömmen är centrala för att förstå hur effektivt isbildningen kan motverkas. Den periodiska tillförseln av värme och turbulens upprätthåller en värmebarriär som hindrar att superkylda droppar fryser fast vid ytan, utan att nödvändigtvis få dem att avdunsta.

Tryck- och hastighetskonturerna förändras också över SJA-cykeln, där sugrörelsen i sugfasen återcirkulerar varm luft in i systemet, vilket förlänger värmeeffekten. Detta skapar en termisk buffert och en tryckmiljö som bidrar till att bibehålla ytan isfri. Numeriska experiment av detta slag kan ligga till grund för optimering av designen av uppvärmda SJA-arrayer, där antalet aktuatorkällor och deras värmeeffekt kan anpassas för att uppnå önskad iskontroll.

Det är viktigt att förstå att effekten av SJA på isbildning inte enbart handlar om mekanisk bortstyrning av droppar utan också om komplex termofluiddynamik där värmeöverföring och flödesdynamik samverkar. Den tidsberoende variationen i SJA-cykeln innebär att både tryck-, hastighets- och temperaturfält förändras, vilket kräver en noggrann analys för att förutsäga och optimera iskontrollen i olika miljöförhållanden. Effekten av uppvärmda jetstrålar kan särskilt inte reduceras till en enkel termisk barriär; istället är den nära kopplad till hur värmen transporteras och blandas i det omgivande flödet.

Förutom det numeriska arbetet är det väsentligt att beakta materialets termiska egenskaper och ytkonfiguration, eftersom dessa faktorer kan påverka hur snabbt värme leds bort och därmed hur effektivt isuppbyggnad förhindras. Samspel mellan den dynamiska flödesmiljön och materialets svar är avgörande för praktiska tillämpningar. I verkliga system kan dessutom yttre faktorer såsom vindförhållanden, luftfuktighet och temperaturvariationer påverka SJA:s prestanda och bör därför integreras i en helhetsbedömning.