Numeriska simuleringar av isbildning på flygplansvingar har visat sig kunna ge prediktioner som ligger mycket nära experimentella och litteraturdata, särskilt när det gäller isens tjocklek vid stagnationspunkten och på trycksidan av vingen. Skillnader på under 1 mm, vilket motsvarar mindre än 0,2 % av vingens chordlängd, har kunnat uppmätas i jämförelse med NASA:s och andra källors resultat. Eventuella avvikelser beror främst på skillnader i modellernas uppbyggnad och parametrar, och i vissa fall på kanterna av isbildningen på sug- eller undertryckssidan.

Ett centralt verktyg i dessa studier är användningen av PCE-metamodeller (Polynomial Chaos Expansion) som ersätter de mer beräkningstunga CFD-simuleringarna. Två olika PCE-metamodeller har utvecklats: en som predicerar maximal istjocklek och isens tvärsnittsarea, och en annan som ger lokala tjockleksvärden vid utvalda punkter längs vingen, framför allt nära framkanten där isbildningen är som mest dynamisk och variabel. Dessa modeller har genomgått rigorös verifiering med R²-värden som ligger mycket nära 1, vilket visar på en exceptionellt hög korrelation mellan CFD-resultat och PCE-prediktioner.

Studierna har också visat att även om inputparametrarna för ytråhet är jämnt fördelade, blir distributionen av isens tjocklek och area icke-uniform, vilket speglar komplexiteten i den fysiska processen. Till exempel koncentreras den maximala istjockleken kring 11 mm medan tvärsnittsarean visar större variationer beroende på de slumpmässiga ytråhetsparametrarna. Detta innebär att statistiska metoder och större provmängder är viktiga för att få en mer fullständig bild av isbildningsscenarier och dess variationer.

Valet av mätpunkter (prober) i simuleringarna har särskild betydelse, eftersom isbildningen är starkt beroende av lokala förhållanden som värmeöverföring och friktion i luftens gränsskikt nära vingens framkant. Att placera fler prober i denna region minskar osäkerheter och ökar simuleringens precision. Genom att använda data direkt från CFD-nätets noder undviks interpolationsfel som annars kan introducera ytterligare osäkerheter.

Metamodellernas validitet har också verifierats genom så kallad leave-one-out (LOO) korsvalidering, med mycket låga felmarginaler, betydligt bättre än vad som är standard i liknande CFD-tillämpningar. Detta ger förtroende för att dessa modeller kan användas för att snabbt och effektivt förutsäga isbildning under varierande förutsättningar, vilket är viktigt för utvecklingen av isavvisande teknologier och säkerhetsbedömningar inom flygindustrin.

Det är avgörande att förstå att även de mest avancerade numeriska metoderna är starkt beroende av kvaliteten på indata och modellens fysikaliska antaganden. Variationer i parametrar som ytråhet, temperaturfördelning och luftflödesförhållanden påverkar starkt utfallet, och därför måste statistisk analys och osäkerhetskvantifiering vara en integrerad del av simuleringarna. Att arbeta med ett stort antal simuleringar och metoder för att hantera variationer, såsom Latin hypercube sampling, är nödvändigt för att få en realistisk bild av isbildningens omfattning och spridning.

Dessutom måste resultaten tolkas med hänsyn till modellens begränsningar och de fysiska skillnaderna mellan simulerade och verkliga förhållanden. Exempelvis kan små skillnader i isens kantformation på sug- och trycksidan ge upphov till skillnader i aerodynamiska egenskaper som inte fullt ut fångas av simuleringarna. Att utveckla och förbättra dessa modeller kräver därför kontinuerlig validering mot experimentella data och uppdatering av fysikaliska parametrar.

Slutligen är det viktigt att inse att numerisk simulering av in-flight isbildning inte bara är en fråga om att förutsäga istjocklek, utan även om att förstå hur isens morfologi påverkar flygplans prestanda och säkerhet. Variationer i isbildningen påverkar luftflödets separation, friktion och därmed aerodynamiska krafter. Detta kräver multidisciplinär förståelse och nära samarbete mellan experiment, simulering och teoretisk modellering för att utveckla pålitliga prediktiva verktyg.

Hur påverkar osäkerhet i driftsförhållanden effektiviteten hos luftprofilers uppsamling av partiklar?

Denna studie analyserar effektiviteten i partikeluppsamling på olika luftprofiler under varierande driftsförhållanden med hjälp av numeriska metoder. Fokus ligger på ett tvådimensionellt NACA23012-luftprofil testat i en kryogen vindtunnel vid NASA Glenn IRT-facilitet. Beräkningarna bygger på SU2-solvern som rekonstruerar det aerodynamiska fältet runt luftprofilen genom att lösa RANS-ekvationerna (Reynolds-Averaged Navier-Stokes) med en finita-volymmetod. Discretiseringen använder en kantbaserad struktur på en dual grid med kontrollvolymer skapade via en median-dual vertex-baserad metod, vilket ger en noggrann och stabil lösning.

Konvektiva flöden approximativt löses med en begränsad andragradsmuscl-scheme och Roe-typens Approximate Riemann Solver, där Venkatakrishnan-Wang-flödesbegränsaren säkerställer numerisk stabilitet. Viskösa flöden beräknas genom ett standardiserat medelvärde av gradienter, medan källtermer rekonstrueras styckevis-konstant inom varje kontrollvolym. Gradienter erhålls via en viktad minsta kvadraters metod, vilket ger exakta lokala värden. Lösningsprocessen drivs till jämvikt med en implicit Euler-metod och konvergens påskyndas med en multigrid-algoritm, där beräkningarna stoppas när resttermen för densitet är åtta storleksordningar mindre än initialt.

Osäkerheten i analysen beaktar fem parametrar: fri ström Mach-tal (M), fri ström tryck (P), fri ström temperatur (T), anfallsvinkel (AoA) och median volymdiameter (MVD) hos dropparna. Dessa parametrar är definierade med uniform sannolikhetsfördelning inom givna osäkerhetsintervall. Osäkerhetskvantifieringen (UQ) genomförs med en Polynomial Chaos-approximator av tredje ordningen, vars koefficienter bestäms via ordinär minsta kvadratmetod baserad på 168 simuleringar genererade med Latin Hypercube Sampling. Surrogatmodellen valideras med ytterligare 50 datapunkter och visar hög noggrannhet med root mean squared error (RMSE) under 5×10⁻⁴ nära stagnationspunkten, något högre vid impingementsgränserna.

Resultaten visar att den mest betydande osäkerhetsbidragaren är dropparnas storlek (MVD), medan anfallsvinkeln (AoA) har betydelse i specifika områden på luftprofilen, särskilt mellan stagnationspunkten och impingementsgränserna. Höga ordningens (HO) effekter, som uppstår vid samtidig variation av flera parametrar, är mest framträdande nära impingementsgränserna, vilket indikerar icke-linjära samverkansfenomen.

En jämförelse mellan surrogatmodellens förutsägelser och experiment visar generellt god överensstämmelse med viss diskrepans vid maximala värden och impingementsgränser, vilket är kritiska områden för isbildning och partikeluppsamling. Liknande metoder och slutsatser gäller för en andra teststudie med justerade osäkerhetsintervall och varierande median droppstorlek, där surrogatmodellen också visar god precision.

Det är väsentligt att förstå att numeriska simuleringar av komplexa flödes- och uppsamlingsprocesser kräver robust hantering av osäkerheter för att realistiskt kunna förutsäga systembeteendet. Osäkerhet i ingångsparametrar, särskilt droppstorlek och anfallsvinkel, kan i hög grad påverka den lokala och totala effektiviteten i partikeluppsamlingen på luftprofiler. Dessutom kan samverkande effekter mellan parametrar leda till oväntade förändringar i resultaten, vilket understryker behovet av avancerade metoder som Polynomial Chaos för osäkerhetsanalys.

Vidare är det viktigt att beakta att denna typ av simuleringar ofta kräver noggranna experimentella data för validering, samt att de numeriska metoderna måste kalibreras och verifieras för att säkerställa trovärdiga resultat. Osäkerheter i själva mätdata och i modellparametrar kan därför bidra till den totala osäkerheten i förutsägelserna. Denna förståelse är avgörande för tillämpningar som flygplansdesign, där isbildning och partikeluppsamling påverkar säkerhet och prestanda.