I simuleringar av flödesfält runt rotorblad används grundläggande fysikaliska storheter som vätskans densitet (ρ), flödeshastighetsvektor (υ), statiskt tryck (p), temperatur (T) och total energi per massenhet (E). SU2-programvaran är ett avancerat verktyg för multiphysiksimuleringar och optimeringar, särskilt inriktat på lösningar av partiella differentialekvationer på ostrukturerade nätverk. Den använder finita volymmetoden (FVM) med en kantbaserad datastruktur på ett dualt nät, där kontrollvolymer konstrueras med en vertexbaserad median-dual metod. I kärnan av SU2 ligger en Reynolds-averaged Navier-Stokes (RANS) lösare, här kombinerad med Spalart–Allmaras turbulensmodell. Tidsintegrationen kan ske implicit för både stationära och icke-stationära problem med en dual-tidsstegsstrategi som ger andra ordningens noggrannhet i både tid och rum. Denna metodik har nyligen visat sig särskilt effektiv för simulering av rotorflöden.

Partikelspårning är en central komponent i simuleringar av moln med superkylda vattendroppar. En Lagrangisk partikelspårningskod utvecklad vid Politecnico di Milano möjliggör modellering på enskild partikelnivå, vilket gör det möjligt att inkludera fenomen som partikelspridning vid kollisioner. Partikelrörelser beräknas utifrån ett system av ordinära differentialekvationer där aerodynamiskt drag, gravitation och partikeltröghet är de huvudsakliga krafterna. Partikelpartikelinteraktioner och deformeringar betraktas dock inte i nuläget, vilket innebär en viss förenkling. Dragkoefficienten och Reynolds tal beräknas genom en kombination av modeller för att bäst överensstämma med experimentella data, vilket ger en realistisk beskrivning av partikelns rörelse genom fluiden. Tekniker för att spåra partiklar genom rörliga nät och över icke-konforma gränssnitt gör det möjligt att simulera komplexa molnflöden runt rotorblad.

Isbildning modelleras med hjälp av programvaran PoliMIce, som implementerar avancerade fysikaliska modeller för isbildningsprocessen. Problemet betraktas som en envägs Stefan-problem som diskretiseras längs kroppens ytor. Masskonservering säkerställs vid cellgränser genom beräkning av mängden ofruset vatten som passerar in och ut ur cellen under påverkan av flöde och tröghetskrafter. En lokal exakt lösning av det ostadiga Stefan-problemet används för temperaturprofilen i islagret, särskilt under glaze-isförhållanden. Isens form beräknas genom en iterativ process där ytan uppdateras med hänsyn till icke-linjär isackretion.

Isavlossning behandlas i en tredimensionell modul inom PoliMIce, där rotorns isyta delas in i flera isblock. När de centrifugala krafterna på isblocken överstiger kombinationen av isens kohesiva och adhesiva krafter sker avlossning. Adhesionsmodeller för isens skjuvspänning vid avlossning bygger på experimentella studier där skjuvspänningen kopplas till temperatur genom en empirisk relation. Kohesionsmodellen visar att sprickbildning initieras vid en kritisk dragspänning och växer tills en kritisk förskjutning uppnås, vilket leder till brott och isavlossning. Detta sker med försumbar viskoplastisk deformation under ackretion.

SRB-II experimentet, utvecklat vid Anti-icing Material International Laboratory i Kanada, har tillhandahållit omfattande data för att studera isbildning och isavlossning under kontrollerade in-flight-icing-förhållanden för små helikoptrar. Testerna omfattade variationer i lufttemperatur, flyghastighet, vätskehalt och droppstorlek i molnet samt rotorvarvtal. Isens påverkan på rotorprestanda dokumenterades genom mätning av effektbehovet och avlossningens karaktär. Adhesionsskjuvspänningen registrerades vid olika temperaturer från rime- till glaze-isförhållanden, vilket visade stark temperaturberoende effekt på avlossningsfenomen.

Det är av yttersta vikt att förstå att isbildning och avlossning är resultatet av komplexa samspel mellan aerodynamiska, termiska och mekaniska faktorer. Numeriska modeller måste därför inkludera detaljerade representationer av fluidflödet, partikelrörelser, värmeöverföring och materialegenskaper hos isen och ytan. De empiriska samband som fastställts i experiment som SRB-II är nödvändiga för att validera och kalibrera modellerna och för att kunna förutsäga isbildningens och isavlossningens dynamik under varierande flygförhållanden. En fullständig förståelse kräver också insikt i begränsningar och förenklingar i dagens modeller, exempelvis antagandet att partikelinteraktioner och viskoplastisk deformation är försumbara, vilket kan påverka noggrannheten i vissa situationer.

Hur olika turbulensmodeller påverkar beräkningar av laminar-turbulent övergång i flygplansvingar

I simuleringar av luftflöde över flygplansvingar spelar turbulensmodeller en avgörande roll för att korrekt förutsäga laminar-turbulent övergång, särskilt vid utvärdering av olika strömningsbeteenden som uppstår på ytan av vingen. Dessa modeller måste ta hänsyn till specifika faktorer som turbulensnivåer, hastighetsfördelningar och tryckgradienter som påverkar det lokala flödet. I fallet med T3AM-modellen, som presenteras i Figur 16, observerades att det fanns en liten skillnad i beräknad Cf-fördelning (skjuvspänningskoefficient) i övergångsområdet, vilket kan förklaras av begränsningar i programmet TEXSTAN. Detta kodsystem kan inte hantera turbulensnivåer lägre än 1,0 %, och därför modifierades den ursprungliga turbulensnivån från 0,87 % till 1,0 %. Detta innebär att de förändringar som observerades kan ligga inom ramen för felmarginalen.

Vid analys av andra fall som T3C2 och T3C5, där tryckgradienten var variabel, visade det sig att den accelererande och deaccelererande flödeshastigheten på den plana plattan orsakade avvikelser i resultaten från experimentella data, särskilt i de områden där flödet ändrar sin hastighet. Här var användningen av integrala analysmetoder (såsom TEXSTAN) mindre exakt jämfört med differentialkoder, som gav bättre överensstämmelse med experimentella resultat. När separationsfenomen inträffade, såsom i TEXSTAN-modellen, ledde det till felaktiga förutsägelser av övergången, eftersom modellen inte kan hantera bubbelströmningsområden.

För att bättre förutsäga dessa övergångsfenomen, utvecklades olika turbulensmodeller som till exempel de som använder modifierade Abu-Ghannam och Shaw (AGS) modeller för att förutsäga övergångens startpunkt under låg turbulens (naturlig övergång). Modifikationer i denna modell gjorde det möjligt att beräkna den kritiska turbulensnivån och den relevanta Reynoldsnummer för att definiera övergångens början mer exakt. AGS-modellen använder en funktion för att beräkna turbulensnivån baserat på en tanh-funktion, där det specifika värdet av denna modell kan påverka noggrannheten i övergången i praktiska tillämpningar.

För modeller som involverar separationsbubblor krävs ytterligare överväganden. Modellen som föreslås av Mayle för att hantera övergångar orsakade av separationsbubblor kan ge bättre resultat än de som använder standard AGS-modellen, särskilt i situationer där separation inträffar tidigare än den turbulensbaserade övergången. Enligt Mayles modell kan även ett korrekt identifierat separationsläge förbättra noggrannheten vid förutsägelser om övergångslängd, vilket gör att övergångsberäkningarna blir mer realistiska i de fall där separationsfenomen spelar en roll.

En annan viktig aspekt är Stratford’s separationskriterium, som är känd för sin precision vid identifiering av separationsbubblor. Denna metod ger mer exakta resultat än integrala metoder, särskilt i situationer där tryckgradienten är negativ, vilket är vanligt vid separation. Stratford-modellen har visat sig vara tillförlitlig för att bestämma separationspunktens läge, vilket är viktigt för att noggrant förutsäga övergångens startpunkt. Däremot kan denna modell ge ett resultat där separationen förutspås något tidigare än den verkligen inträffar, vilket kan leda till mindre precis förutsägelse av övergångens placering i vissa situationer.

Slutligen är det viktigt att förstå att dessa modeller inte bara påverkar själva förutsägelsen av övergången utan också den övergripande noggrannheten för hela flödesanalysen. Integrala metoder, även om de är användbara för att beskriva allmänna trender i flödet, kan vara mindre exakta i områden med starka negativa tryckgradienter, vilket gör att användning av mer detaljerade modeller som Stratford eller Mayles modeller kan vara nödvändiga för att säkerställa ett korrekt resultat. Det är också viktigt att uppmärksamma att dessa modeller har sina begränsningar i olika flödesregimer, som exempelvis när det gäller påverkan av isbildning och dess inverkan på övergången. Här visar de flesta modeller på en otillräcklig hantering av isbildningsfenomen i samband med turbulens och separation, vilket gör att ytterligare forskning på detta område fortfarande behövs.

Hur modelleras issbildning på kablar vid frysnederbörd och vilka utmaningar finns i numeriska simuleringar?

När regndroppar faller vertikalt med terminalhastighet och träffar en cylinder, såsom en kabel, sker delvis frysning medan en del vatten fortfarande rinner längs kabelns yta under påverkan av gravitationen. Denna process skapar så kallad runback-is, där rinnande vatten fryser under kabeln och bildar små strömmar av is som syns tydligt på digitala skanningar och fotografier. Eftersom inte alla droppar fryser omedelbart på kabelns yta, når en del av det ofrusna vattnet kabelns lägsta punkter, där istappar formas och hänger nedåt. Trots svårigheter för digitala skanningstekniker att fånga detaljrikedomen i istapparna, framträder dessa tydligt på fotografier. Modellen lyckas väl med att förutsäga både isens utbredning på kabelns översida och runback-isens formation på undersidan, liksom placeringen och längden på istapparna. Även om den simulerade istjockleken tenderar att överskattas något, är överensstämmelsen mellan experimentella observationer och numeriska förutsägelser anmärkningsvärt god.

Under mer än tre decennier har morfogenetiska modeller för issbildning visat sig vara både effektiva och tillförlitliga för att simulera komplexa isbildningsmönster. Dessa modeller förenklar det omfattande multiphase-fenomenet genom att betrakta isbildningen som en successiv uppbyggnad av partiklar, där fysikaliska, termodynamiska och hydrodynamiska processer översätts till stokastiska regler och procedurer. Tack vare denna förenkling är det möjligt att genomföra omfattande beräkningar av isbildning på vanliga bärbara datorer, vilket är en av modellens stora fördelar. Dock begränsas den beräkningsmässiga effektiviteten av den relativa ineffektiviteten i beräkningar av luftflöde och värmeöverföring, vilket förhindrar snabb omberäkning när isbildningen fortskrider och förändrar sin form.

Framtidens beräkningsmöjligheter, såsom kvantdatorer, kan komma att ersätta morfogenetiska modeller med mer sofistikerade partikelsimuleringar som Smooth Particle Hydrodynamics. Detta skulle möjliggöra en mer detaljerad och direkt lösning av det fullständiga multiphase-problemet kring issbildning. Men detta är en utveckling som ligger flera decennier framåt, varför morfogenetiska modeller förväntas förbli viktiga verktyg inom isbildningsforskningen under överskådlig tid.

En central aspekt för fortsatt utveckling är vikten av att noggrant testa och validera dessa modeller mot såväl vindtunnel- som fältmätningar i verkliga flygförhållanden. Detta kräver även utveckling av nya metoder för att jämföra komplexa och dynamiskt föränderliga isformationer under varierande atmosfäriska förhållanden. Det är viktigt att förstå att isbildning inte är en statisk process, utan snarare en dynamisk och mångfacetterad företeelse där samverkan mellan nedbörd, temperatur, luftflöde och materialets egenskaper skapar varierande och komplexa strukturer.

För att fullt ut förstå och förutse issbildning är det också nödvändigt att ta hänsyn till mikroskaliga processer såsom partikelinteraktioner och övergångar mellan fasförändringar. Den termiska balansen mellan värmeavgivning vid frysning och värmeavledning genom konvektion och ledning styr mycket av isens tillväxtmönster. Därför måste både makroskopiska och mikroskopiska processer modelleras noggrant för att simuleringarna ska ge realistiska resultat.

Slutligen är det viktigt att förstå att dessa modeller inte bara är akademiska övningar utan har praktiska tillämpningar inom flygindustrin, energisektorn och infrastrukturområdet, där isbildning kan innebära säkerhetsrisker och ekonomiska kostnader. En djupare insikt i de fysikaliska och numeriska processerna bakom issbildning bidrar till att utveckla bättre strategier för förebyggande och bekämpning av isrelaterade problem.