2. Hálógenerálás

A numerikus háló elemei

2.1.ábra A numerikus cella elemei: csomópontok (kék), hálóvonalak (fekete), az egyik határfelület (sraffozott narancs), cellaközéppont (narancs).

A 2.1.ábrán a háló néhány alapösszetevője látható egyetlen elem esetében. A csomópontokat X,Y,Z koordinátáikkal, az elemek oldalfelületeit pedig felületen található csomópontok sorszámával adhatók meg. Ezeket az információkat a hálógeneráló program az outputként létrehozott hálófájlban (*.msh) tárolja. A cellákat és cellacentrumokat a megoldó a hálófájl beolvasását követően építi fel a csomópontok és felületek alapján. FLUENT megoldó esetében a hálóvonalak a csomópontokat összekötő egyenes szakaszok, az elemek oldalai a csomópont által meghatározott határfelületek.

A FLUENT rendszer minden mezőváltozónak pontosan egy diszkrét értékét tárolja a cellák középpontjának megfelelő pozíciókhoz; ezek a diszkretizációval nyert egyenletrendszer ismeretlenjei. A mezőváltozók felületi és csomóponti értékeit, amennyiben szükségesek, a cellaközépponti értékekből interpolálja a megoldó.

Lehetséges geometriai modellek

·         2D tengelyszimmetrikus

·         2D tengelyszimmetrikus, perdületes

·         2D síkáramlás

·         3D áramlás

A tengelyszimmetrikus modellek esetében mindig az x tengelyt tekintjük forgástengelynek. A cellákhoz rendelt térfogat a cellafelület x tengely körüli forgatásával előálló toroid térfogata. 2D tengelyszimmetrikus, perdületes modellek esetében a FLUENT három sebességösszetevőt határoz meg, melyek két geometriai koordinátától függhetnek.

Kétdimenziós síkáramlás esetében is értelmezett a cellák térfogata. Ilyen esetben a tartomány „z” irányú mérete mindig 1 m.

Elemtípusok

2.2.ábra Egyszerű elemtípusok

Alaptípusoknak tekinthetők kétdimenzióban a tisztán háromszögelemekből, vagy tisztán négyszögelemekből álló hálók, háromdimenzióban pedig a tetraéder, hexaéder, illetve prizma (háromszög alapú hasáb) hálótípusok. Egyes hálózási módszerek zónánként, eltérő alaptípusokból un. hibrid hálót hoznak létre. 3D modell esetében a piramis elemeket csak négyszög hálóval felbontott felület és térfogati tetra elemek közötti átmeneti rétegben alkalmaz az ANSYS hálógenerálója.

A FLUENT rendszer a 2.2.ábrán látható egyszerű elemeken kívül, 6-nál nagyobb oldalszámú (poliéder) elemeket is megenged. Poliéder elemek többféle módon létrejöhetnek:

·         Adaptálással (2. fejezet későbbi pontján részletezve);

·         Duális háló készítésével  (2. fejezet későbbi pontján részletezve);

·         Nem-konform interfészek összekapcsolásával (3. fejezetben kifejtve).

A numerikus háló mérete

Általában nem az egyenközű numerikus háló szolgáltatja a leghatékonyabb numerikus megoldást. A numerikus hibák egyenletes eloszlására kell törekedni, ezért ott kell finomabb felbontás, ahol várhatóan gyorsan változik a megoldás. A 2.3.ábrán az áramlási tér olyan pontjai láthatók, ahol nagy gradiensek jelenlétére számítani kell.

2.3.ábra Hálósűrítést igénylő pontok

·         Ha az áramlás követi a felület görbületeit, akkor görbületi középpont felé haladva a nyomás és sebesség gyorsan változik.

·         Hézagokban legalább 4 cellasor alkalmazása szükséges a sebességprofil megfelelő közelítéséhez.

·         A fali határrétegben a sebesség rohamosan változik. A fal melletti cellasorban általában falüggvényt (ld. 4. fejezet) használ használunk. E cellák  mérete erősen befolyásolja a számítás eredményét.

·         Nyírórétegekben a sebesség és a turbulens jellemzők a rétegre merőlegesen rohamosan változnak, továbbá az elején a réteg irányában haladva is.

Új alkalmazási területeken, amíg elegendő modellezési tapasztalat nem áll rendelkezésre, ajánlatos próbahálózást végezni, majd az első futási eredmények tükrében továbbfejleszteni a hálót.

A sűrítést igénylő részek lehetnek pontszerűek, vonalszerűek, vagy rétegszerűek. Az elemtípust és a sűrítés módszerét ettől függően kell megválasztani, figyelembe véve a sűrítés anizotrópiáját.

A negatív nyomáscsúcsok jelenléte miatt sűrítésre lehet szükség a szilárd felületek erősen görbült részeinél főként, ha az áramlás követi a szilárd test vonalát (azaz nem válik le a határréteg). Ezt a szempontot az ANSYS hálógeneráló szoftvere automatikusan is képes figyelembe venni (Advanced Sizing Function).  

Az átlagos hálómérethez képes szűk hézagok esetében általában azt szeretnénk elérni, hogy a modell pontosan szolgáltassa a hézagon átáramló térfogatáram értékét, mely a globális áramképet jelentősen befolyásolhatja. Ilyen réseket találunk például áramlástechnikai gépekben a lapátok és a ház között: a lapátrésben átáramló folyadék leválasztja a határréteget a lapát szívott oldalán ezzel jelentősen módosítva a gép hatásfokát és jelleggörbéjét. Amennyiben a résben turbulens áramlás is kialakulhat, célszerű legalább 4 réteg cellát alkalmazni, lamináris áramlás esetében akár 3 réteg is elegendő lehet a résáram megfelelő pontosságú számításához. 

Turbulens áramlásban a fal melletti cellákban a leíró egyenletek közelítése helyett félempirikus kiegészítő modellek alkalmazása szokásos, – ld. falfüggvények a „Turbulencia” c. fejezetben – ezért például 4 cellaréteg esetében a tényleges numerikus módszer csak 2 rétegben működik.

Az áramlási sebesség a falon, a fal sebességével tekinthető azonosnak (tapadás törvénye). A zavartalan áramlás és a fal közötti határréteg vastagságát az áramlásra jellemző Reynold-szám határozza meg: Re növekedtével a sebességváltozás egyre rohamosabb a falnál. Turbulens falfüggvények alkalmazása esetében a falfüggvény fajtája határozza meg a fa melletti cella falra merőleges maximális (továbbá esetleg minimális) méretét.

A falra merőleges irányú sűrítést a fenti példában prizmatikus elemekkel értük el. E hálótípus előnyeit: az anizotróp (csak falra merőleges irányú) sűrítés és a háló áramvonalasításának lehetőségeit a továbbiakban részletezzük.

Nyírórétegnek tekintjük azokat a rétegeket az áramlási térben, ahol a sebesség réteggel párhuzamos komponense jelentősen eltér a réteg két oldalán. Nyíróréteg leggyakrabban a határréteg leválásával jöhet létre, pl. tompa testek mögött. A fali hálóhoz hasonlóan a nyírórétegek is anizotróp sűrítést igényelnek, továbbá áramlás irányú sűrítés is szükséges a nyíróréteg elején, a turbulens jellemzők gyors változása miatt. Meg kell jegyeznünk, hogy a nyírórétegek egyedi felbontására bonyolult geometriájú gyakorlati esetekben (például egy városrész épületei körüli külső áramlás vagy egy gépjármű motorház belső áramlásának modellezésekor) nincs lehetőség, ilyenkor egyenletes sűrű hálót szokás létrehozni a szilárd testek tágabb környezetében.

A háló torzultsága

A fluxusok numerikus közelítésének pontosságát (rendjét) az határozza meg, hogy a megoldás diszkrét közelítésének hibája a cellaközéppontok közötti távolság hányadik hatványával arányos. Magasabb rendű sémák alkalmazása esetén a háló sűrítésével a hibák gyorsabban csökkennek feltéve, hogy a háló kellő mértékben „ortogonális”, továbbá a hálóméret kevéssé változik a szomszédos cellák között, tehát a háló eléggé „sima”.

Az alábbi jellemzők alapján szoktuk a háló minőségét értékelni:

2.4.ábra Torzultság (balra) és hirtelen méretugrás (jobbra)

1.       Ortogonalitás (Ortogonality):

A szomszédos cellák középpontját összekötő szakasznak a felületi normálissal bezárt szöge.

2.       Torzultság (Skewness):

Az elem belső szögeitől függő mennyiség: szabályos sokszög elemek esetében 0, ha valamelyik belső szög sokkal kisebb a többinél, akkor a skewness értéke 1-hez tart.

3.       Simaság (hirtelen méretváltozás):

A szomszédos cellák méretének aránya.

Akár egy háromszög háló is lehet a fenti értelemben ortogonális. Ilyenkor a felülettel határos cellákban tárolt mezőváltozók különbségét a középpontok távolságával osztva pontosan számítható a konduktív fluxus felületre merőleges komponense, mely éppen a cellafelületi integrálok számításához szükséges komponens. 

A háló ortogonalitásának fenti mértéke nem cellát, hanem cellapárt vagy cellafelületet jellemez, ezért nehezen kiértékelhető és megjeleníthető. Emiatt a gyakorlatban a háló torzulását legtöbbször a skewness paraméterrel értékeljük. A megengedhető legmagasabb értékek hexa elemek esetén 0.85, tetra elemek esetén 0.95.

2.5.ábra A hálóminőség ellenőrzése ANSYS Mesher rendszerben

A háló minőségére vonatkozó elemzéseket az ANSYS Meshing hálógeneráló programban a háló paraméterlapján a „Statistics” paraméterek között a „Mesh Metric” kiválasztásával végezhetjük. Az egyes minőségi osztályokba tartozó elemeket a diagram oszlopaira kattintva jeleníthetjük meg. A kritikus cellák alakjából és helyéből általában sejthető a probléma oka is.

A háló torzulásának okai lehetnek a helytelen hálóstruktúra (hexa elemek alkalmazása esetén), a hálóméret gyors változása, rövid élek vagy hegyes felületek jelenléte a geometriában, továbbá a szűk rések vagy erősen tagolt geometriai modellek viszonylag durva hálóval történő felbontása. Egyes esetekben a cellák torzulása csak egy-két pontban jelentkezik. Ilyenkor célszerű megpróbálkozni a hálóméretek kismértékű megváltoztatásával. Szükség lehet a geometriai modell egyszerűsítése is (CAD Clean-up) melynek lehetséges módjait a későbbiekben részletezzük. 

A hirtelen méretváltozás problémája falak és nyírórétegek hálózásakor szinte minden esetben fellép. Ha a cellák magassága hirtelen változik, akkor a numerikus sémák elmélet pontossága nem valósul meg (pl. Second Order Upwinding módszer esetén a diszkretizációs hiba nem csökken az elemméret négyzetével arányosan), ezért törekedni kell a cellánként 30%-nál kisebb méretnövekedésre. Tolerálható legnagyobb méretnövekedés: 100%. Kompromisszumot kell keresnünk a háló simasága és a háló elemszámának csökkentése között.

A háló áramvonalasítása

A diszkretizációs hiba egy fontos formája a fals diffúzió:

A 2.6.ábrán látható példában egy passzív megmaradó mennyiség (pl. egy kis koncentrációjú szennyező) transzportját vizsgáljuk a hálóvonalakhoz képest ferde irányban, a transzportegyenlet megoldásával, állandó áramlási sebesség és 0 értékű vezetési tényező feltételezésével. A pontos megoldás egy éles kontúrokkal rendelkező piros csík lenne kék mezőben.

2.6.ábra A fals diffúzió megjelenése különböző súrúségű numerikus hálókon. A 2D modelltartomány mérete balról jobbra haladva (mindkét sorban): 10x10, 20x20, 40x40, 80x80.

A kezdeti háló felbontását egymást követően háromszor megduplázzuk. A kék, zöld és piros színek 0%, 50%, és 100% értéket jelentenek. Látható, hogy a First Order Upwinding módszer alkalmazása esetén a csóva elfogadhatatlan mértékben elmosódik. Ez a numerikus diffúzió okozza általában a térbeli diszkretizálásból eredő legjelentősebb hibát. A numerikus diffúzió mértéke jelentősen meghaladhatja a leíró egyenletekben foglalt fizikai diffúzió mértékét. Hasonló probléma merül fel minden transzportegyenletben, ezért a numerikus modellek túlbecsülik a viszkozitás és a hővezetés intenzitását is. Second Order Upwinding séma alkalmazásával a hiba igen nagymértékben csökkenthető.

Ha a fenti példában a csóva vízszintes irányú lett volna, akkor a numerikus diffúzió egyáltalán nem lépett volna fel még First Order Upwinding sáma alkalmazása esetén sem, tehát áramvonalasított hálóval is hatékonyan küzdhetünk a fals diffúzió ellen. Ehhez prizmatikus, vagy hexaéder elemekből álló hálót kell létrehozni, melynek cellahatárfelületei áramlás irányában rendezett módon, csak áramlással párhuzamosan és áramlásra merőleges irányban helyezkedjenek el. Fontos lehet a háló áramvonalasítása lekerekített testek esetében, ahol a határréteg leválásának helyét – így a test ellenállását is – a határréteg által közvetített csúsztatófeszültség határozza meg, mivel ez utóbbi jellemző érzékeny a fals diffúzió okozta viszkozitási hibára.  

A fals diffúzió kiküszöbölésére tehát az alábbi módszereket alkalmazhatjuk:

·         Second Order Upwindinget vagy más, nagy nagyobb formális pontosságú fluxus séma alkalmazása és a háló sűrítése;

·         A háló áramvonalasítása (pl. fali határrétegekben és szabad nyírórétegekben).

A numerikus fluxus sémák sajátosságaival kapcsolatban ajánlható alapvető szakirodalmi forrás: [Hirsch, 1988].

 

Fali inflációs réteg

2.7.ábra Fali hálósűrítés. Elvi elrendezés a cellaközéppontokkal (balra), ennek megvalósulása a gyakorlatban (jobbra)

·         Prizmatikus réteg (hexa vagy wedge elemek);

·         A fali cellasor igen lapos (legalább a logaritmikus fali rétegig lemegy, azaz y+<300);

·         Legalább 3 (inkább 5) réteg;

·         Fokozatos méretnövekedés (max. 50% cellasoronként);

·         Sima átmenet a belső tér hálójába.

Inflációs rétegek beillesztésére az ANSYS hálózójában speciális eszköz található.

Az inflációs réteg elemei – a fali háló típusától függően – prizma vagy hexa elemek lehetnek. A inflációs réteg méreteit megadhatjuk a rétegek számával, a növekedési rátával, továbbá vagy a fali réteg magasságával, vagy a teljes inflációs réteg magasságával, vagy a belső hálóba való sima átmenet előírásával.

Akkor áll fenn a sima átmenet a belső tér hálójába, ha az inflációs réteg legfelső celláinak középpontja ugyanakkora távolságra van az inflációs réteg határától, mint a belső (2.7.ábrán látható példában háromszög) háló legalsó elemeinek középpontja.

Nem mindig tudják a hálózó programok folytonossá tenni az inflációs réteget a határfelület éles sarkainál, továbbá problémát okozhat az inflációs réteg kialakításában, ha a háló túlságosan durva a geometriai modell tagoltságához képest. Ilyen esetben a háló sűrítésével vagy geometriai modell egyszerűsítésével segíthetjük a hálógeneráló program munkáját.

Hexa hálók

2.8.ábra Egy belső bordázattal rendelkező hőcserélő cső CFD elemzéséhez készítet, blokk-struktúrált háló [Hernádi, 2015]

·         Általános esetben (kb. 30-50%-al) kevesebb cellával érhető el azonos pontosság, jó minőségű háló esetén.

·         Anizotróp (egyirányú) sűrítés lehetséges: határrétegek, nyírórétegek, rések hálózhatók.

·         Áramvonalasított hálók készíthetők.

·         Többirányú sűrítéshez bonyolult háló topológiát kell felépíteni. Blokkstruktúrált hálót célszerű, készíteni: ICEM CFD hálózó javasolható. Ez különösen fontos pl. ventilátorok, szivattyúk, repülőgépek modellezése esetén.

Struktúrált háló alatt azt értjük, ha a háló tégla topológiával rendelkezik, azaz cellasorokra, oszlopokra és rétegekre bontható. Blokkstruktúrált háló esetében ez az egyes blokkokra igaz. A térfogati háló sűrűsége az egyes blokkok éleinek alkalmas felosztásával szabályozható.

Sajnos jó minőségű hexa háló ma még bonyolultabb geometriák esetében nem készíthető automatikus módszerekkel. A tartomány alkalmas feldarabolásáról a felhasználónak kell gondoskodnia, ami meglehetősen idő- és munkaigényes. Főként az okoz nehézséget, hogy a szomszédos blokkok határfelületein a felosztásnak egyezni kell. Megéri az erőfeszítést, ha nagy pontosság elérése a cél, ugyanis ezzel a hálózási módszerrel áramvonalasíthatjuk a hálót és az irányonkénti (anizotróp) hálósűrítéssel lényegesen nagyobb felbontást érhetünk el adott cellaszám mellett.

Blokkstruktúrált hálók esetében a felbontás O-blokkok beszúrásával szabályozható finoman. Az O-blokkok záródhatnak önmagukba a tartomány belsejében, vagy akár kifuthatnak a tartomány határára (nyitott O-blokk). Ezekre alkalmazására láthatunk példát a 2.8.ábrán, ahol a bordázat irányába sűrűsödő hálót egy zárt és két nyitott O-blokk beiktatásával sikerült kialakítani. Minden O-blokk tetszőleges számú rétegre felbontható, így egy újabb sűrítési paramétert vezet be.

Tetra és hibrid hálók

2.9.ábra Egy uszoda külső és belső légáramlásának vizsgálatára alkalmazott hibrid háló: tetra alapháló hexa magokkal

·         Háromirányú sűrítés egyszerűen végezhető, azonban egyirányú sűrítés nem lehetséges. (Háló torzulást okoz.)

·         Hexa hálóval jól összekapcsolható, kivéve, ha lapos hexa cellákat kell körbevenni vele.

·         Igen jó automatikus módszerek vannak rá, amelyek prizmatikus inflációs réteggekkel (pl. fali határréteg hálóval) is képesek kombinálni az alaphálót. 

·         Adaptáláskor kevésbé nő a cellaszám mint hexa cellák esetében.

 

A tetraéder elemekből álló (tetra) háló fő előnye, hogy elkészítése lényegesen kevesebb manuális munkát igényel, emellett a háló sűrűsége méretfüggvényekkel egyszerűen szabályozható.

A tetra háló jobb adaptálhatósága (azaz a megoldáshoz igazodó, utólagos lokális sűríthetősége) és cellatorzulásra való viszonylagos érzéketlensége deformálódó mozgóhálós modellek készítését is lehetővé teszi komplex mozgást végző testek körül.

A tetra hálók hátránya, hogy viszonylag izotróp felbontást igényelnek, ezért például nem tudjuk kihasználni, hogy egy csatornaáramlásban a mezőváltozók a csatorna tengelye irányában sokkal lassabban változnak, így ebben az irányban sokkal durvább felbontás is elegendő lenne. Ugyancsak korlátozást jelent, hogy a jelenleg elérhető hálózó eszközökkel a tetra hálókat nem tudjuk áramvonalasítani.   

Duális (poliéder) hálók

Tetra hálóból a 2.10.ábrán látható poliéder hálót alakíthatunk ki FLUENT-ben a Mesh/Polyhedra menüpont alkalmazásával. Az átalakítás eredményeként a cellaközéppontok és a csomópontok szerepet cserélnek.

2.10.ábra Tetraéder cellákból képezett duális (poliéder) háló

·         Patch Conforming hálózó használata célszerű, a Patch Independent módszerrel torzult poliéderek keletkeznek.

·         Zónálisan is alkalmazható.

·         A cellaszám kb. 60%-al csökken, a felületek száma azonban nem csökken, ezért közel azonos a műveletigény iterációs lépésenként, viszont sokkal kevesebb iterációs lépés kell (gyorsabb konvergencia).

·         Igen torz tetra hálóból is jó minőségű poliéder háló állítható elő: a torzulás kiküszöbölésének hatékony módszere.

·         Az effektív térbeli felbontás számottevően csökken, azonban a memória igény nem csökken jelentősen, ezért leginkább lokálisan célszerű alkalmazni.

Adaptált hálók

FLUENT rendszer Adapt menüpontjában elérhető parancs segítségével helyileg sűríthetjük a hálót, akár a megoldás sajátosságainak figyelemebe vételével is, ily módon a háló képes a megoldáshoz alkalmazkodni. A művelet során az eredeti cellák középpontjába elhelyezett új csomópont felhasználásával a cella több – hexa elemek esetében nyolc – részre bomlik. Az átalakítást jellemzően csak a tartomány kisebb részén alkalmazzuk. A 2.11.ábrán látható példában a sűrítés három egymást követő lépésben történt.

2.11.ábra Egy olajtranszformátor durva modellezése adaptált hálóval

ANSYS-FLUENT-ben az adaptálás visszafordítható, tehát az utólag beiktatott extra csomópontok eltávolításával vissza lehet állítani az eredeti hálót. Az adaptálás nem csak a felhasználó által megadott térrészekben végezhető, hanem megadhatunk összetett adaptálási kritériumokat is, például a mezőváltozók gradiensére vonatkozó határokat, ezért a háló időben változó megoldáshoz is alkalmazkodhat, például mozgó nyíróréteghez vagy lökéshullámhoz.

Hátrányai, hogy túl nagy (100%-os) ugrást hoz létre a cellaméretben, ezért egy jó minőségben elkészített hálóhoz képest növeli a diszkretizációs hibát, és annak becslését is megnehezíti.

Hálózási módszerek

2.12.ábra Hálózási módszerek ANSYS mesherben

1.       Sweep: A forrásfelület (Source Face) hálóját „extrudálja” a tartományba.

2.       Multi zone: Hasonló Sweep módszerhez, de annál összetettebb testeket is tud hálózni. Automatikus módon, blokkstruktúrált hálót próbál előállítani a megadott térrészben.

3.       Hex dominant: Kívülről befelé parkettáz. A térrész közepén beiktat néhány tertra cellát is, melyek gyakran rossz minőségűek.

4.       Patch Conforming tetra: A felületről kiindulva befelé növeszti a tetra hálót.

5.       Patch Independent tetra: Szabályos kocka hálóból készít tetra hálót, majd a felület közelében besűríti. Geometriai hibákra érzéketlen.

6.       Cut-cell: Adaptívan sűrített hexa hálót készít. A felület közelében a hexa elemeket elvágja. Az Assembly Meshing módszerek között található.

 

A hálózási módszert ANSYS Meshing programban Method Control beszúrásával szabályozhatjuk. Egyes térfogatokra (testekre) előírhatjuk a megfelelő egyedi hálózási módszert.

Viszonylag egyszerű alakú térrészek esetében blokkstruktúrált hálót a geometria alkalmas feldarabolását követően sweep módszerrel lehet létrehozni. Ilyen esetben az élek sűrítési módjának megadásával (Edge Sizing) alakíthatjuk ki a kívánt térfogati hálót. Ha egy összetett geometriai modell esetében blokkstruktúrált háló előállítására van szükség – például áramlástechnikai gépek vagy egy repülőgép vizsgálata esetén – az ICEM hálózó használata javasolható az ANSYS Meshing program helyett.

Több test hálózása esetén az ANSYS Meshing program először a Sweep és Multi Zone módszerrel hálózható testekkel kezdi. Az többi térrész esetében a Patch Conforming tetra háló az alapértelmezett. A testek különállóan is hálózhatók, így a hálózás sorrendjét is meghatározhatja a felhasználó.

A Patch Conforming és a Pach Independent tetra hálózók közül az utóbbi robusztusabban működik, annak köszönhetően, hogy az azonos peremfeltételi zónához (Named Selection) tartozó felületek közötti határvonalakra ez a módszer nincs tekintettel. E tulajdonságának köszönhetően a kontúrfelület élei általában nem határozott élként, hanem kissé töredezett formában jelennek meg a numerikus hálón.

A mérnöki gyakorlatban sok esetben szükséges rövid idő alatt összetett rendszerekre vonatkozó eredményt produkálni. Erre adnak megoldást az automatikus Assembly Meshing (öszeállítás hálózó) módszerek, melyek – figyelembe véve a felhasználó által előírt sűrítési sémákat és inflációs rétegeket – általában elfogadható minőségű hálót képesek létrehozni. Megjegyzendő, hogy az automatikus hálózó módszerek jellemzően jóval nagyobb cellaszám alkalmazásával képesek azonos pontosságot biztosítani, mint a szakszerűen optimalizált hálók, ezért a manuális munka csökkenése mellett a szimulációs kapacitásigény növekedésével is számolni kell. Assembly Meshing alkalmazása esetén nem használhatunk testenként eltérő hálózási módszert, csak két, globálisan alkalmazott módszer közül választhatunk:

A.     Cut-Cell

B.       Tetrahedrons (Patch Comforming)

Hálózás előkészítése

Jó minőségű háló előállítása érdekében egy CAD rendszerben elkészített, összetett geometriai modellen a következő előkészítő lépések végrehajtására lehet szükség:

a)       A geometriai modell alkalmas elmetszése a hálózás megkönnyítésére. Pl. Sweep módszerrel hálózható tartományok létrehozása.

b)      Rövid élek, keskeny felületek kiszűrése hálózáskor:

·         Model/Insert Virtual Topology

·         Mesh/Pinch (globális)

·         Mesh/Insert Pinch (lokális)

c)       Hálóméret szabályzása:

·         Mesh/Sizing (globális)

·         Mesh/Insert Sizing (lokális)

d)      Inflációs réteg (felületre merőleges sűrítés) létrehozása:

·         Mesh/Inflation (globális)

·         Mesh/Insert Inlaltion (lokális)

Sok esetben célszerű a geometriai modell előkészítésekor (Design Modeller-ben) olyan metszéseket alkalmazni, melyek a tartomány minél nagyobb részében hengerszerű, eltolási szimmetriával rendelkező testeket hoz létre. Ezek a térrészek Sweep vagy Multi Zone módszerrel hálózhatók, a tartomány többi részén pedig valamelyik tetra hálózó vagy Hex Dominant hálózó módszer használható.

Összetett geometriai modellek esetében előnyös lehet kialakítani egy egyszerű geometriával rendelkező (hexa) magot a tartomány belsejében, melyet – a 2.9.ábrán látható módon – kívülről tetra háló vesz közrül. A határfelület közelében jellemző geometriai részletekhez az automatikusan generált tetra háló könnyebben tud alkalmazkodni. Az így elkészített háló szilárd falak mentén a háló prizmatikus réteggel is kiegészíthető.

A rövid élek és keskeny felületek jelenléte torzult cellákat eredményez. A geometriai modellből célszerű eltávolítani minden olyan részletet, amely a cellák méreténél kisebb.

Virtuális topológia alkalmazása esetén a felhasználó manuálisan összeolvaszthat egyes határfelületeket és éleket. A Pinch eszköz az egyszerűsítéseket egy megadott minimális méret figyelembevételével a hálón végzi el.

A háló méretét (felbontását) előírhatunk térfogatokra, felületekre, élekre vagy pontokra. A méret lokális előírását az alapháló (globális) méretéhez képesti sűrítésre szokás alkalmazni. Pont körüli sűrítéshez megadható a befolyásolt tartomány (Sphere of Influence) sugara. Élek méretezése esetében a felbontás változhat az él mentén: sűrűsödhet bármelyik végpont, vagy az él középpontja felé.

Inflációs zóna létrehozásakor ki kell választanunk a testet és azon határfelületeit, ahol inflációs réteget kívánunk elhelyezni. Célszerű elkerülni, hogy az inflációs réteg a hálózott térrész belsejében érjen véget, tehát törekednünk arra, hogy az réteg vége a határfelületre essen, vagy a prizmatikus réteg önmagába zárójon. 

Szintén a hálózás előkészítéseként szokás (Named Selection-ként kijelölve) nevesíteni a geometriai modell egyes testeit, egyes határfelületeket és belső felületeket. A nevesített zónákon később térfogati forrásokat, peremfeltételeket vagy belső szakadási feltételeket alakíthatunk ki, továbbá felhasználhatjuk azokat az eredmények kiértékelésekor. A nevesített felületek kontúrélet a Patch Independent tetra hálózó és az Assembly hálózó módszerek is figyelembe veszik, azokra háló csomópontokat helyeznek, egyéb élek esetében ez nem szükségszerű.

Kellemetlen hibákhoz vezethet a felületek utólagos, hálótopológiát is érintő módosítása, ha már megkezdtük a modellezést FLUENT-ben, mivel a FLUENT a peremfeltételi beállításokat konkrét azonosítóval rendelkező peremfeltételi zónákhoz rendeli hozzá, mely azonosítók ilyenkor megváltoznak. E probléma elkerülésére a peremfeltételi zónák módosítását követően célszerű a Workbench-ben a Setup objektumra jobb gombbal, majd a Reset menüpontra kattintva törölhetjük a FLUENT műveleti adatait a modelladatbázisból. Ilyenkor a FLUENT-ben végzett beállításokat meg kell ismételni. Nem szükséges ezt elvégezni, amennyiben a módosítások csak geometriai paramétereket érintik, azonban a háló topológiája változatlan marad. Ilyen esetben, a megoldó a korábbi beállításoknak megfelelően, automatikusan képes ismételni az elemzést.

Deformálódó háló

A deformálódó hálók lehetővé teszik folyadék-szilárd kölcsönhatás (FSI) modellek felépítését az áramlástani megoldó és valamilyen mechanikai megoldó összekapcsolásával. Merevtest-szerű mozgás – például forgó szeleptest vagy egy szabadon zuhanó test – modellezéséhez az ANSYS-FLUENT beépített megoldóval (Six Degree of Freedom Solver, 6DOF) rendelkezik, összetettebb problémák – például szívbillentyű vagy lengő hídszerkezet [Szabó, ????] – modellezése a FLUENT és az ANSYS mechanikai megoldó Workbench-ben történő összekapcsolásával lehetséges. Kapcsolt modellekben az áramlástani megoldó a határfelületi erőket, a mechanikai megoldó pedig a felület elmozdulást szolgáltatja.

A háló mozgatása idő-függvényében megadott módon változó áramlási terek, például belsőégésű motor, vagy a 2.13.ábrán látható víztároló modellezését is tehetővé teszik.

2.13.ábra Egy ivóvíz tároló medence modellezése mozgó hálóval. A és B: numerikus háló; lentebb: a vízszint napszakos változása.

Deformálódó tartományt a perempontok sebességének előírásával adhatunk meg. A belső csomópontok mozgatását a FLUENT automatikusan elvégzi, erre háromféle módszer alkalmazható:

a)       Struktúrált hexa háló alkalmazásával rétegek eltávolítása/beszúrása (Layering). Ez a módszer alkalmazható például dugattyús gépek szimulációjára.

b)      A belső háló simítása rugó modell alapján.

c)       Lokális újrahálózás a torzult részeken.

b) és c) módszerek kombinálhatók is.