Vid optimering av elektrotermiska iskydds system är en av de största utmaningarna att minska isbildningens mängd samtidigt som man håller sig inom givna energibudgetar. Ett specifikt exempel på en sådan optimering är minimisering av den så kallade "runback ice formation", eller den is som samlas längs bakre delen av vingen efter att ha frusit vid främre delar. Detta problem innefattar att anpassa värmeflödet på en iskydds yta (IPS) för att säkerställa en effektiv avsmältning av is utan att överskrida den tillåtna effekten, i detta fall 4815 W/m per enhetsspann.

Optimeringen kan formaliseras genom att minimera en norm (L-2 norm), som representerar det totala värdet av isackumulationen på hela ytan. För att undvika lokala minima i optimeringsproblemet är det vanligt att genomföra flera optimeringar med olika startpunkter i designvariabelrummet. I det aktuella fallet genomfördes tio optimeringar, och de fem bästa körningarna visades i figur 7. De visade på en snabb konvergens till ett förbättrat resultat, där samtliga optimeringsförsök konvergerade till liknande minimivärden inom mindre än 600 funktionsutvärderingar.

En av de mest intressanta observationerna från denna optimering är hur värmeflödet fördelas på olika värmeelement. I den optimerade konfigurationen reducerades värmeflödet vid främre kanten (värmare A), medan större mängder värme tillfördes längre ner på vingen där konvektiva förluster är lägre. Detta återspeglas i den optimerade värmeflödesfördelningen som visas i figur 8, där värmeflödet för den optimerade lösningen markant skiljer sig från den ursprungliga baslinjen.

En viktig aspekt av denna optimering är att man inte bara strävar efter att minimera isbildning, utan även att säkerställa en robust prestanda vid osäkra och varierande väderförhållanden. För att hantera dessa osäkerheter används en robust optimeringsmetod där flera objektiv beaktas. Ett av de viktigaste målen är att minska de mest allvarliga isbildningarna samtidigt som man maximerar sannolikheten för att inga isbildningar sker överhuvudtaget. Detta kräver att man beaktar variationer i väderförhållanden såsom molnighet, luftfuktighet och temperatur.

Vid robust optimering används ofta en Monte Carlo-sampling för att beakta osäkerheter i de ingående parametrarna. Denna metod innebär att man upprepade gånger simulerar systemet under olika osäkra förhållanden för att skapa en sannolikhetsfördelning för isbildningen. Genom att ta hänsyn till dessa osäkerheter kan man anpassa designen så att den fungerar effektivt under både de mest ogynnsamma och mest gynnsamma förhållandena.

En annan viktig aspekt som kan bidra till optimeringen av iskyddssystem är att tillämpa en modifierad version av den L-2 norm som används för att minimera isbildning. I den robusta optimeringsmodellen kan man anpassa målvariablerna för att fokusera på att minska de allvarligaste isbildningarna, till exempel genom att använda en 95:e percentil av isackumulationen, vilket skulle återspegla en design som klarar 95% av de minst gynnsamma väderförhållandena.

För att förbättra effekten av optimeringsprocessen kan dessutom vissa tekniker som surrogate-modeller eller svarsyteanalys användas för att minska antalet nödvändiga simuleringar. Dessa metoder möjliggör att man skapar förenklade modeller som kan användas för att approximera resultaten av mer komplexa simuleringar och på så sätt effektivisera hela processen.

För att uppnå långsiktig hållbarhet i dessa system bör inte bara prestanda för själva iskydds systemet beaktas utan även dess miljöpåverkan. Effektivitet och minskad energiförbrukning är avgörande för att minska både driftkostnader och miljöpåverkan. Detta innebär att det också är viktigt att utveckla system som kan anpassas till nya material och teknologier för att ytterligare förbättra effektiviteten.

Hur påverkar Biot- och Stefan-tal frysningsprocessen för överkylda droppar under flygning?

Konvergens till minst fem signifikanta siffror uppnås för temperaturvärden inom droppar med relativt låga Biot-tal, där Bi ¼ 0,1146, Bim ¼ 0,0047 och Bir ¼ 0,0012, vilket visar att temperaturfördelningen i droppen är relativt homogen. Detta förklarar varför antagandet om rumsligt uniform temperatur, som användes i experimenten av Hindmarsh et al. (2003), fungerar väl i dessa fall. Vidare analyserades två hypoteser om hur isen bildas och fördelas under fasövergången: antingen som en isring på droppens yta eller som en jämnt fördelad isvolym inuti droppen. Resultaten visade att frysningsstadiet varade ungefär 24,11 sekunder för den första hypotesen och 23,60 sekunder för den andra, vilket överensstämmer rimligt med experimentella observationer, där tiden uppskattades till cirka 20 sekunder.

Biot-talet, en dimensionlös parameter som beskriver förhållandet mellan termisk resistans för ledning och konvektion vid droppens yta, har en avgörande betydelse för värmeöverföringen och därmed för frysningsprocessens dynamik. För att testa metodens tillämpbarhetsgränser undersöktes olika Biot-värden (0,1, 1,0, 5,0 och 10,0) med hjälp av modeller baserade på reducerade representationer (CIEA) och fullständiga numeriska lösningar (GITT). Lägre Biot-tal resulterade i mer radialsymmetriska och homogena temperaturprofiler, vilket gynnar noggrannheten hos de förenklade modellerna. Vid högre Biot-värden blev approximationerna mer osäkra, särskilt för de enklare modellerna som inte inkluderade temperaturderivat, medan mer avancerade modeller kunde bibehålla god överensstämmelse med den fullständiga numeriska lösningen.

Utöver Biot-talet är Stefan-talet en annan väsentlig dimensionslös parameter som beskriver förhållandet mellan sensibel och latent värme i fasövergångsprocesser. Variationer i Stefan-talet påverkar i hög grad frysningsförloppets karaktär. En ökning av Stefan-talet från 0,11 till 0,20 innebär att den sensibla värmeöverföringen ökar i förhållande till den latenta värmeutvecklingen, vilket leder till större temperaturvariationer under frysningsprocessen och en kortare frysningsperiod. Detta är särskilt betydelsefullt i tillämpningar där temperaturen förändras snabbt och där latent värme måste hanteras noggrant för att förutsäga frysförloppets tidsram och struktur.

Den förenade analysen av Biot- och Stefan-tal visar hur dessa dimensionlösa parametrar påverkar både värmeöverföringen och fasövergångens dynamik i överkylda droppar. För lägre Biot-tal kan enklare modeller användas med hög noggrannhet, medan högre Biot-tal kräver mer komplexa representationer för att korrekt beskriva temperaturfältet och frysningsförloppet. Dessutom indikerar variationer i Stefan-talet vikten av att balansera sensibel och latent värme vid simulering och kontroll av frysningsprocesser, särskilt i industriella och naturvetenskapliga sammanhang där fasförändringar är kritiska.

Det är viktigt att förstå att även om modellerna kan approximera de termiska förloppen väl, så påverkas noggrannheten kraftigt av parametrarnas värden och modellens förmåga att beskriva rumsliga temperaturgradienter och fasgränser. Att införliva experimentella data för att kalibrera och validera dessa modeller är därför nödvändigt för att säkerställa tillförlitliga förutsägelser i praktiska tillämpningar. Vidare bör man beakta att komplexiteten i den fysiska processen kan kräva flerfaldiga fasövergångar eller icke-uniforma isbildningar, vilket kan utmana antaganden om uniformitet och kräva avancerade numeriska metoder.