A motor beömlőjében két kulcsfontosságú terület van, amelyek esetében fontos meghatározni a konvekciós hőátadási koefficienseket: az egyik a piccolo csőből származó forró levegő sugárainak és a belső géptető élének kölcsönhatása, a másik pedig a szabad áramlás és a géptető kiemelkedő éle közötti kapcsolat. Downs és James (1988) modellje szerint a külső áramlási mezőt egy henger elülső részére körüljáró áramlás és két párhuzamos lap mentén folyó áramlás kombinációjaként ábrázolták. A modell célja az volt, hogy az áramlások közötti eltéréseket minimalizálja, ezáltal csökkentve a numerikus eredmények és a tesztadatok közötti eltérést.

Riley (1991) később egy CFD (számítógépes áramlástani) eszközt alkalmazott a hőátadási koefficiensek javítása érdekében, így a kapott szilárd felületi hőmérséklet közelebb került a tesztadatokhoz. Henry (1989) és később Henry (1992) egy numerikus kódot dolgozott ki a jegesedésmentesítő rendszerek paramétereinek előrejelzésére. A kód a hőmérsékleti egyensúlyi eljárások módosításával próbálta pontosabban modellezni a jegesedésmentesítő rendszerek működését. Ezen kívül Al-Khalil (1991) matematikai modellt fejlesztett ki, amely figyelembe veszi a folyadékfolyás és a felületek közötti rivulátumok hatásait, melyek szintén szerepet játszanak a hőmérsékleti egyensúly fenntartásában.

Al-Khalil és kollégái (2001) továbbfejlesztették a szimulációt, és beépítették az ANTICE kódot, amely a LEWICE kód áramlásoldó és csepptrajektória-számító rutinfunkcióit használja. Az ANTICE kód célja a hőmérsékleti egyensúly szimulálása volt, figyelembe véve a légi járművek szárnyfelületén eloszló hőáramot, valamint a folyékony vízfilm lebomlását rivulátumokká. A modell alkalmazása során figyelembe vették a különböző levegőáramlatokat és a cseppimpulzusokat is. Morency et al. (1999a, b) két változatban is megvalósították a CANICE nevű numerikus kódot, amely a jegesedésmentesítő rendszerek működésének szimulációjára szolgált.

A CANICE A változata az áramlás és hőátadás kísérleti hőátadási együtthatóját alkalmazza, míg a CANICE B a hőmérsékleti gradientek hatásait is figyelembe veszi, miközben a turbulens hőátadás integrális elemzését végzi. A modell alkalmazása során a szimulációk figyelembe veszik a laminaris és turbulens áramlás közötti átmenetet, valamint a hőhatások változását a hőmérsékleti réteg mentén. A hőmérsékleti egyensúly szimulációja azokat a komplex folyamatokat modellezi, amelyek az áramlási paraméterek és a felületi hőmérsékletek kölcsönhatásain keresztül alakítják a rendszer működését.

A különböző modellek és numerikus kódok alkalmazása során figyelembe kell venni a folyékony víz filmrétegének lefolyásakor fellépő rivulátum-formálódás hatásait is. Az előzőekben ismertetett modellek több változata figyelembe veszi ezt a dinamikát, és próbálja meghatározni a folyamatot, amely során a vízfilm egyes részei eltérnek a felületektől, így rivulátumokká alakulnak. Ezen rivulátumok viselkedésének modellezése az áramlásdinamika szempontjából fontos, mivel azok hatással vannak a hőátadás és a hőmérsékleti eloszlás előrejelzésére.

A legújabb fejlesztések, mint a Gosset (2017) által alkalmazott film-átrendeződési modell, amely az OpenFoam CFD kódot használja a folyadékfilm szétesésének szimulálására, lehetőséget biztosítanak arra, hogy pontosabban előre jelezzük a vízfilm viselkedését és annak hatásait a hőátadásra. A film szerkezete, például az, hogy milyen irányban és milyen sebességgel változik, szoros kapcsolatban áll a hőmérsékleti és áramlási paraméterekkel, ami segíthet a jegesedésmentesítő rendszerek hatékonyságának pontosabb meghatározásában.

A hőátadás szimulációk során a numerikus modellek alkalmazásánál fontos, hogy a különböző áramlási régiók (lamináris, átmeneti, turbulens) hatásait helyesen modellezzük. A megfelelő kódok és algoritmusok alkalmazása lehetővé teszi a hőátadás pontos előrejelzését, ami alapvető a hővédelmi rendszerek, például a jegesedésmentesítő rendszerek fejlesztésében.

A legfontosabb, hogy az áramlás és a hőmérsékleti hatások közötti kölcsönhatásokat nem lehet egyszerűsíteni túlzottan, mivel ezek a rendszerek rendkívül komplexek, és a legapróbb változások is jelentős hatással lehetnek a végső eredményekre. A különböző modellek folyamatos fejlesztése, illetve az újabb kísérleti adatok figyelembevétele elengedhetetlen ahhoz, hogy a jegesedésmentesítő rendszerek és a hőátadási modellek minél pontosabb előrejelzéseket adjanak a valós környezetben.

Hogyan zajlik a túlhűtött cseppek fagyása: a jégképződés és a fázisátalakulások részletes modellje

A túlhűtött vízcseppek fagyása során a jégképződés két fő hipotézist követhet: az egyik szerint a jégkezdemény a csepp felületén alakul ki, amely hőmérsékletileg hidegebb a belső részénél, így először ott éri el a nukleációs hőmérsékletet. Ez a folyamat egy gömbhéjszerű jégréteg kialakulásához vezet, amely idővel befelé terjed. A másik elképzelés homogén nukleációt feltételez, ahol a folyékony és szilárd fázisok keveredése egyenletes, opálos megjelenésű jég- és vízkeveréket eredményez a csepp egész területén.

Hindmarsh és munkatársai által javasolt képlet (2003) alapján a képződött jég térfogata a következő összefüggéssel becsülhető: a jég térfogata arányos a csepp térfogatával, a folyadékfázis fajhőjével és sűrűségével, valamint a fagyáspont és a nukleációs hőmérséklet közti különbséggel, miközben figyelembe veszi a jég sűrűségét és a fázisátalakulás latent hőjét. A gömbhéj-modell esetében az interfész kezdeti helyzetét a jég héj térfogatából számítják ki, míg homogén eloszlás esetén a latent hőt egy módosított értékkel, az Lm-mel helyettesítik, amely figyelembe veszi a víz-jég keverék jellemzőit.

A jégképződés harmadik szakasza a szilárdulás, amely során egy időben változó, egy-dimenziós hővezetési egyenlet írja le a folyamatot a csepp sugara mentén. A határfeltételekben szerepelnek a konvektív hőátadás, az evaporatív tömegátadás és a sugárzásos hőátadás együttes hatásai, miközben a szilárd-folyékony határfelület hőmérséklete a fagyáspontot követi. A szilárdulási front helyzete a jégképződés kezdeti modelljétől függően változik: a héjszerű jégképződésnél az interfész kezdeti helyzete az előzőleg számított R_ini, homogén eloszlásnál pedig a csepp külső sugara. Ezen szakasz leírására dimenziómentes paramétereket használnak, például a Stefan számot, amely az anyag fázisátalakulásának jellemzője.

A negyedik szakaszban a már megszilárdult jég tovább hűl, mely folyamatot hasonló, de a rögzített térfogaton értelmezett hővezetési egyenlet írja le, ismét konvektív, evaporatív és sugárzási hatásokkal. A jég hőmérséklet-eloszlása a szilárdulás befejező időpontjától indul, és a környező levegő hőmérséklete, valamint az anyag fizikai tulajdonságai (pl. fajhő, hővezetési tényező) befolyásolják.

A jégképződés és -szilárdulás modelljeinek pontos megértése elengedhetetlen a légkörfizikában, például a felhőképződés, csapadékképződés és jegesedési folyamatok tanulmányozásához. Az egyes fázisok dinamikájának részletes numerikus szimulációja segít feltárni a túlhűtött cseppek viselkedését repülőgépek szárnyán, jégvirág-képződésnél, vagy fagyos eső kialakulásakor.

Fontos szem előtt tartani, hogy a jégképződés során felszabaduló vagy elnyelődő latent hő jelentős befolyásoló tényező, amely nemcsak a hőmérsékletváltozást, hanem a jégkristályok növekedési ütemét és szerkezetét is meghatározza. A hő- és tömegátadási együttes hatásai komplex kölcsönhatásokat hoznak létre, melyeket a valós légköri körülmények és a cseppméret változékonysága tovább nehezít. A modelleknek ezért figyelembe kell venniük a sugárzásos hőveszteséget, a vízpára jelenlétét és az áramlási viszonyokat is, mivel ezek a paraméterek jelentősen módosíthatják a jégképződés dinamikáját.

Az anyagtudomány szempontjából a jég és víz közötti fázisátmenetnek mikroszkopikus jellege is kritikus, különösen a nukleáció kezdeti lépéseiben, ahol a helyi hőmérséklet és felületi energiák befolyásolják a kristálymagok kialakulását. Ezen folyamatok makroszkopikus megértése segítheti a jég elleni védekezési stratégiák fejlesztését, például a repülésbiztonság vagy az elektromos vezetékek jegesedésének megelőzése terén.

Hogyan modellezhetjük a repülés közbeni jéglerakódást a Level-Set módszer alkalmazásával?

A repülés közbeni jéglerakódás modellezése alapvető fontosságú terület, amely többféle fizikai problémát ötvöz. Jéglerakódás akkor jön létre, amikor egy repülőgép vízcseppekkel telített felhőbe kerül, és ezek a cseppek a szárnyak, szárnyfülek, illetve más felületek mentén fagyva jégformákat képeznek. Ez a jelenség jelentős hatással van az aerodinamikai jellemzőkre, mivel a jég kialakulása csökkenti a felhajtóerőt és növeli a légellenállást. A repülőgépek tervezése és biztonsági tanúsítása szigorú jégmentesítő tesztekhez kötött, amelyek célja annak biztosítása, hogy a gép minden körülmények között biztonságosan működjön.

A numerikus jégmodellek fejlesztése napjainkban különösen fontos, mivel a nagy teljesítményű számítástechnikai technológiák óriási lehetőségeket kínálnak ezen a területen. Az egyszerű konfigurációk mellett széles körű paraméterek összehasonlítására is lehetőség nyílik, ami az experimentális kutatás szűk határait jelentősen kiterjeszti. Ezen túlmenően, a numerikus modellek lehetőséget adnak arra, hogy azokat az erősen bonyolult jégformákat is előre jelezzük, amelyek gyakran előfordulhatnak a valódi repülés során, miközben költséghatékonyabbak, mint a széltunnel tesztelés. Azonban, még mindig nagy kihívást jelent a háromdimenziós jégmodellek alkalmazása a számítási hatékonyság és robusztusság terén, mivel a kézi újrahálósítási folyamatok nehézkessége miatt nem mindig sikerül helyes jégformát előállítani.

A hagyományos repülés közbeni jéglerakódási modellek a jégképződés folyamatát egy szekvenciális, közel állandó állapotú lépésekre bontják, és az egyes lépéseket külön-külön számítják. A modellek az alábbi szakaszokat követik:

  1. Légáramlás: Az első lépésben a levegő áramlásának modellezése történik. Ehhez különböző típusú áramlásmodelleket használnak, például a viszkózus áramlások esetén a Reynolds-átlagolt Navier-Stokes (RANS) egyenleteket, míg a nem-viszkózus áramlásokat Euler vagy potenciál áramlásmodellekkel számolják.

  2. Vízcseppek becsapódása: A következő lépésben a vízcseppek egy felületre csapódnak, amihez Lagrange- vagy Euler-típusú cseppmodell alkalmazható. Az ilyen típusú modellek feltételezik, hogy a vízcsepp azonnal eltűnik a felületen, és az azt követő jégformálás a folyamatot függetlenül kezeli más részekről.

  3. Fagyás és folyadékfilm képződés: Az egyes cseppek a felületen vízréteget képeznek, amely fokozatosan megfagy. A hődinamikai modellek azt számítják ki, hogy mennyi víz fagy meg, és mennyi folyik vissza a felületen.

  4. Jégformák kialakulása: Az utolsó lépésben a jégformák előállítása következik, amely az előző lépésekhez igazodó hálózatfrissítést igényel. Ez a jégképződés folyamatát különböző remeshing (újrahálósítás) eszközökkel modellezik.

Az ilyen típusú modellek fő feltételezése, hogy a légáramlás nem zavaródik meg az összes lépés alatt, és a jég kialakulása követi az eredetileg számított értékeket. Azonban ez a feltételezés nem teljesen helytálló, mivel a repülés közbeni jéglerakódás dinamikus jelenség, és a jég folyamatosan változik az idő előrehaladtával. A modern jégmodellek próbálnak ezen a problémán javítani azáltal, hogy a jeget rétegről rétegre modellezik, ahol minden egyes réteg kialakulása egy-egy időbeli lépésként értelmezhető.

A Level-Set módszer az egyik legújabb technika, amely lehetővé teszi, hogy hatékonyan modellezzük a jéglerakódás dinamikáját. A Level-Set módszer előnye, hogy csökkenti a hagyományos modellek egyik legnagyobb problémáját, a hálózatok újrahálósításának bonyolultságát. A Level-Set alapú szimulációk segítségével a jégformák rendkívül komplex geometriák, például összefonódó vagy erősen elválasztott felületek kezelésére is képesek.

A Level-Set módszer egyszerűsíti a jégképződés több lépcsőből álló szimulációját. A szilárd határfelületet különböző módszerekkel, implicit vagy explicit módon kezelhetjük. A szilárd részeket, például a szárnyat és a jeget, a Level-Set függvény segítségével definiáljuk, és az interfészt Lagrange-pontok listájaként ábrázoljuk. Az így modellezett légáramlási határfeltételek a fal mentén alkalmazott turbulenciával kerülnek figyelembevételre, míg a csepptranszport Euler-féle megközelítéssel, diszkrét formulák alkalmazásával történik.

A jéglerakódás modellezése tehát egy rendkívül bonyolult és több lépcsős folyamat, amely különböző szimulációs eszközöket igényel. A Level-Set módszer jelentős előnye, hogy képes hatékonyan kezelni a bonyolult jégformákat és dinamikusan változó geometriákat anélkül, hogy szükség lenne a hagyományos újrahálósításra. Az ilyen típusú modellek alkalmazása lehetővé teszi a repülés közbeni jégképződés pontosabb előrejelzését, és hozzájárulhat a repülőgépek biztonságának növeléséhez. Azonban fontos figyelembe venni, hogy a jéglerakódás rendkívül komplex jelenség, és a modellek mindig is csak egy megközelítést adhatnak, amelynek pontosítása és finomítása folyamatos kutatást igényel.