CFD-validatie van turbulente stroming rond NACA 0012
Het doel was te toetsen of een CFD-simulatie van de luchtstroming rond een vleugelprofiel betrouwbaar is: de lift kwam binnen 1–2% overeen met windtunnelmetingen van NASA, en een overschatting van de weerstand met 40% komt hoogstwaarschijnlijk voort uit één corrigeerbare meshkeuze. Tweedimensionale RANS-simulatie van turbulente stroming rond een NACA 0012-vleugelprofiel bij een hoog Reynoldsgetal.
- Project
- Casestudy, Cornell University (onlinecursus)
- Mijn rol
- Geometrie, meshing, CFD en validatie
- Software
- ANSYS Meshing, ANSYS Fluent, CFD-Post
- Methode
- 2D stationaire RANS, Realizable k-ε
- 1,059liftcoëfficiënt, ≈1–2% onder de NASA-referentie
- +40%overschatting van de weerstand, waarschijnlijke oorzaak gevonden
- 6×106Reynoldsgetal op basis van de koorde, bij een invalshoek van 10°
- 3/4validatie- en verificatiecontroles vervuld
Aangetoonde vaardigheden
Aanpak
Nagaan of een stationaire RANS-CFD-methodiek, vanuit de grondbeginselen opgebouwd en toegepast op een goed gedocumenteerd referentieprofiel, lift en weerstand nauwkeurig genoeg voorspelt om haar te vertrouwen bij een nieuwe, minder goed bekende geometrie.
-
De berekening opzetten
ANSYS Meshing
Eén symmetrische 2D-profieldoorsnede in een groot verre-velddomein, opgezet om één stationair werkpunt te isoleren voordat aan verfijning van de mesh of het turbulentiemodel wordt begonnen. Het NACA 0012-profiel heeft een koorde van 1,0 m en het domein reikt in elke richting 12,5 koordelengtes vanaf het profiel, met de vrije-stroomsnelheid en turbulentiegrootheden op de inlaat in het verre veld, statische druk op de uitlaat en de no-slipconditie op het profieloppervlak.
- Reynoldsgetal
- 6×106
- Invalshoek
- 10°
- Koorde
- 1,0 m
- Omvang van het domein
- 12,5 koordelengtes
Randvoorwaarden op het domein: vrije-stroomsnelheid en turbulentiegrootheden op de inlaat in het verre veld, statische druk op de uitlaat en de no-slipconditie op het profieloppervlak. 
Schets van het verre-velddomein, dat in elke richting 12,5 koordelengtes vanaf het profiel reikt. 
Gemaatvoerd NACA 0012-profiel waarmee het stromingsdomein is opgebouwd, koordelengte 1,0 m. -
De mesh opbouwen
ANSYS Meshing
Een gelaagde meshstrategie lost de grenslaag en het zog vlak achter het profiel op zonder die resolutie overal in een domein van 12,5 koordelengtes te hoeven betalen: een grof verre veld dat overgaat in een verfijningszone met een invloedsbol van 3c, een randelementgrootte van 0,3% c op het profieloppervlak en 10 inflatielagen met een groeifactor van 1,2.
De orthogonale kwaliteit ligt voor het overgrote deel van de elementen boven 0,9; het kleine aantal cellen van mindere kwaliteit zit precies bij de achterrand, waar de ongestructureerde mesh plaatselijk scheef wordt.
- Cellen
- 26.962
- Inflatielagen
- 10, groeifactor 1,2
- Celgrootte aan het oppervlak
- 0,3% van de koorde
- Verfijnde zone
- 3 koordelengtes rond het profiel
Grove mesh in het verre veld die overgaat in de fijnere verfijningszone met invloedsbol rond het profiel. 
Close-up van de oppervlaktemesh van het profiel en de inflatielagen in de grenslaag. -
De stroming doorrekenen
ANSYS Fluent
De stroming wordt gemodelleerd als stationair, onsamendrukbaar en tweedimensionaal, gesloten met het turbulentiemodel Realizable k-ε en standaardwandfuncties. De discretisatie verloopt in twee stappen: eerst een robuuste berekening met eerste-orde-upwind, daarna, zodra het stromingsveld zich heeft gevormd, een nauwkeurige berekening met tweede-orde-upwind, waarbij alle residuen tot 10−6 convergeren.
Het drukveld toont het verwachte stuwpuntsgebied net onder de neus en vlak daarachter een sterke zuigpiek aan de bovenzijde; de snelheidsvectoren bevestigen dat de stroming tot aan de achterrand aan beide zijden aanliggend blijft.
- Turbulentiemodel
- Realizable k-ε
- Wandbehandeling
- Standaardwandfuncties
- Discretisatie (definitief)
- Tweede-orde-upwind
- Residuen
- 10−6
Druk rond de voorrand: stuwpuntsgebied (rood) en zuigpiek (blauw). 
De snelheidsvectoren bevestigen een aanliggende, voorwaarts gerichte stroming aan beide zijden, zonder loslating. -
Vergelijken met gegevens van NASA
CFD-Post
De berekende verdeling van de drukcoëfficiënt is over de volledige koorde rechtstreeks vergeleken met gepubliceerde windtunnelmetingen van NASA voor dit profiel, met een polylijn langs de koorde die is opgebouwd met de boundary-intersection-functie van CFD-Post.
De twee curves vallen vrijwel samen, inclusief de scherpe zuigpiek bij de voorrand. Dat is in deze studie het sterkste bewijs dat het drukveld dat de lift genereert goed wordt weergegeven, los van hoe goed de grenslaag zelf is opgelost.
- Liftcoëfficiënt
- 1,059 tegenover ≈1,07–1,08
- Weerstandscoëfficiënt
- 0,0167 tegenover ≈0,012
CFD-drukcoëfficiënt (lijn) tegenover gepubliceerde windtunnelmetingen van NASA (punten), weergegeven met de gebruikelijke omgekeerde verticale as.
Controle van de resultaten
Van de vier validatie- en verificatiecontroles die voor deze berekening het zwaarst wegen, zijn er drie zonder meer vervuld; de vierde wijst op één specifieke, te verhelpen beperking.
3/4
validatie- en verificatiecontroles vervuld
De liftcoëfficiënt en de drukverdeling over het oppervlak komen tot op enkele procenten overeen met gepubliceerde windtunnelmetingen van NASA, en de oplossing behoudt massa met een relatieve onbalans die ver onder elke technische relevantie ligt. De weerstandscoëfficiënt overschat de experimentele waarde met ongeveer 40%; de y+-diagnose hieronder herleidt dat hoogstwaarschijnlijk tot een specifieke, corrigeerbare beperking van de mesh en niet tot een fout in het turbulentiemodel.
CFD-resultaat tegenover gepubliceerde experimentele waarde
Krachtcoëfficiënten als percentage van de experimentele referentiewaarde; 100% is exacte overeenstemming.
| Validatie- / verificatiecontrole | Berekend resultaat | Referentie / streefwaarde | Status |
|---|---|---|---|
| Drukcoëfficiënt aan het oppervlak (Cp) | Goede overeenstemming over de volledige koorde | Komt overeen met de gepubliceerde Cp-verdeling van NASA | Vervuld |
| Liftcoëfficiënt (CL) | 1,059 (≈1–2% te laag) | ≈ 1,07–1,08 (NASA-referentie) | Vervuld |
| Massabehoud | ≈4,8×10−7 kg/s bij een doorstroming van ≈1490 kg/s | Netto fluxonbalans ≈ 0 | Vervuld |
| Weerstandscoëfficiënt (CD) | 0,0167 (≈40% te hoog) | ≈ 0,012 (NASA-referentie) | Aandachtspunt, waarschijnlijke oorzaak gevonden |
Resolutie nabij de wand en y+-diagnose
Omdat al is aangetoond dat het drukveld nauwkeurig is, moet de overschatting van de weerstand ergens anders vandaan komen: de mesh nabij de wand. De dimensieloze wandafstand y+ geeft aan waar het middelpunt van de eerste cel ligt ten opzichte van de lagen van de turbulente grenslaag, en de standaardwandfuncties in deze studie zijn alleen nauwkeurig als die eerste cel in de logaritmische laag valt, ruwweg tussen 30 en 300.
De berekende y+ piekt precies bij de voorrand tot ongeveer 90 en zakt vlak daarachter tot 3 à 5, maar ligt over het grootste deel van de koorde op slechts 20 tot 25, onder het vereiste bereik. Daardoor valt de eerste cel in de bufferlaag in plaats van in de logaritmische laag, waar de logaritmische wandwet die standaardwandfuncties opleggen niet geldt, zodat de berekende wandschuifspanning aan nauwkeurigheid inboet. Dat sluit aan bij de overschatte weerstandscoëfficiënt die in de hele oplossing doorwerkt.
y+ langs de koorde
Technisch oordeel en aanbevelingen
Naast het oordeel voldoet/voldoet niet wijst deze studie op één dominante hefboom om het resultaat te verbeteren, en op een korte lijst punten die het controleren waard zijn voordat de methodiek op een nieuw profiel wordt toegepast.
-
Aandachtspunt
Verfijn de mesh nabij de wand voordat de weerstand wordt vertrouwd
De afwijking in de weerstand komt hoogstwaarschijnlijk van een wandmesh die de eerste cel in de bufferlaag plaatst in plaats van in de logaritmische laag. Verfijning van de inflatielagen tot y+ ≈ 1, of ten minste y+ ≈ 5 in combinatie met Enhanced Wall Treatment, zal naar verwachting de weerstandsvoorspelling aanzienlijk verbeteren zonder de al nauwkeurige lift- en drukresultaten te verstoren.
-
De meshkwaliteit aan de achterrand vraagt plaatselijk aandacht
Een klein aantal scheve cellen met een lage orthogonale kwaliteit ligt geconcentreerd bij de achterrand. Ze vormen een klein deel van de totale mesh, maar zitten precies waar de grenslagen van beide zijden samenkomen; het loont daarom ze samen met de wandverfijning hierboven aan te pakken.
-
Meshonafhankelijkheid en omvang van het verre veld staan nog open
Deze studie gebruikte één mesh en een domeinomvang van 12,5 koordelengtes. Voor geen van beide is formeel aangetoond dat verdere verfijning het resultaat niet verandert; beide blijven dus op de lijst voordat het resultaat als volledig geverifieerd kan gelden.
-
Bevestig dat de methodiek algemeen toepasbaar is voordat ze wordt hergebruikt
Voordat deze aanpak voor meshing en turbulentiemodellering op een nieuwe, minder goed bekende geometrie wordt toegepast, is het zinvol deze berekening te herhalen met een tweede turbulentiemodel (k-ω SST of Spalart-Allmaras) en bij extra invalshoeken, om na te gaan dat de huidige goede overeenstemming niet specifiek is voor dit ene werkpunt.
Technische details
Methode, mesh- en solverinstellingen en de analyse van het stromingsveld, voor wie zich verder wil verdiepen.
Methode, mesh en solveropzet Modelaannames, meshinstellingen, solverconfiguratie en meshkwaliteit
Methode
De stroming wordt beschouwd als stationair, onsamendrukbaar en volledig turbulent, en opgelost met de Reynolds-gemiddelde Navier-Stokesvergelijkingen in plaats van elke turbulente wervel rechtstreeks op te lossen.
Sluiting met wervelviscositeit: τt = μt (∂ū/∂y), met μt uit het turbulentiemodel Realizable k-ε
- Het opbouwen van het domein en de oppervlaktemeshing zijn volledig in ANSYS Meshing uitgevoerd, uitgaande van één symmetrisch NACA 0012-profiel.
- De stroming wordt gemodelleerd als stationair, onsamendrukbaar en tweedimensionaal, gesloten met het turbulentiemodel Realizable k-ε en standaardwandfuncties.
- De discretisatie verloopt in twee stappen: eerst een robuuste berekening met eerste-orde-upwind, daarna, zodra het stromingsveld zich heeft gevormd, een nauwkeurige berekening met tweede-orde-upwind.
- De scope is beperkt tot één werkpunt, één turbulentiemodel en één meshtopologie; meshonafhankelijkheid en reeksen bij meerdere invalshoeken zijn aangemerkt als vervolgwerk en hier niet uitgevoerd.
Meshgeneratie en solveropzet
Een gelaagde meshstrategie lost de grenslaag en het zog vlak achter het profiel op zonder die resolutie overal in een domein van 12,5 koordelengtes te hoeven betalen.
| Gebied | Instelling | Waarde |
|---|---|---|
| Ver veld | Elementgrootte achtergrond | Grof, overgaand in de verfijningszone |
| Omgeving van het profiel | Straal van de invloedsbol / elementgrootte | 3,0 m (3c) / 5,0×10−2 m (5% c) |
| Profieloppervlak | Randelementgrootte / biasfactor | 3,0×10−3 m (0,3% c) / 10,0 |
| Grenslaag | Inflatielagen / groeifactor / max. dikte | 10 / 1,2 / 6,0×10−3 m (0,6% c) |
| Meshtotalen | Knopen / elementen | 27.343 / 26.962 |
De solver is vervolgens ingesteld volgens het wiskundige model hierboven:
| Waarde | Waarde |
|---|---|
| Solvertype / tijd / ruimte | Drukgebaseerd / stationair / 2D vlak |
| Turbulentiemodel | Realizable k-ε, standaardwandfuncties |
| Dichtheid / viscositeit van het fluïdum | 1,1767 kg/m³ / 1,009×10−5 kg/(m·s) |
| Discretisatie (definitief) | Druk: tweede orde; impuls, k, ε: tweede-orde-upwind |
| Convergentiecriterium (definitief) | 10−6, alle residuen |
Vóór de solveropzet is de meshkwaliteit op twee standaardmaatstaven gecontroleerd. De orthogonale kwaliteit ligt voor het overgrote deel van de elementen boven 0,9; het kleine aantal cellen van mindere kwaliteit zit precies bij de achterrand, waar de ongestructureerde mesh plaatselijk scheef wordt en dat zonder verdere lokale onderverdeling moeilijk te voorkomen is. Deze zwakte aan de achterrand is meegenomen in de aanbevelingen.
Analyse van het stromingsveld Snelheid, druk, stroomlijnen en turbulentie, vergeleken met het gedrag dat vóór de berekening werd verwacht
Elke controle hieronder isoleert één deel van het geconvergeerde stromingsveld en vergelijkt dat met het kwalitatieve en kwantitatieve gedrag dat werd verwacht voordat de solver ooit draaide.
A1. Snelheidsveld en ontwikkeling van de grenslaag
Het veld van de snelheidsgrootte versnelt rond de neus en langs de bovenzijde (zuigzijde), precies zoals te verwachten is bij een positieve invalshoek, en vertraagt geleidelijk naar de achterrand zonder enig teken van grootschalige loslating.
Door de mesh over het snelheidsveld te leggen, wordt zichtbaar dat over een groot deel van de koorde slechts één controlevolume de volledige grenslaagdikte overspant: een probleem van meshresolutie dat in de y+-controle kwantitatief terugkomt. De snelheidsvectoren (stap 3) bevestigen dat de stroming tot aan de achterrand aan beide zijden aanliggend blijft, waarmee overtrek of loslating als verklaring voor de overschatte weerstand vervalt.
A2. Drukveld en zuigpiek
Het drukveld toont het verwachte stuwpuntsgebied net onder de neus en vlak daarachter een sterke zuigpiek aan de bovenzijde, waar de stroming rond de kromming van de voorrand het sterkst versnelt.
Richting achterrand stijgt de druk weer, wat de ongunstige drukgradiënt oplevert die verantwoordelijk is voor de verdikking van de grenslaag in A1. Omdat de lift vrijwel volledig uit deze drukverdeling voortkomt en niet uit viskeuze schuifspanning, is een nauwkeurig drukveld de belangrijkste reden dat de liftcoëfficiënt zo dicht bij de experimentele waarde uitkomt. Het drukveld is weergegeven in stap 3.
A3. Stromingsstructuur: stroomlijnen en turbulentieproductie
Stroomlijnen die stroomopwaarts van het profiel beginnen en gekleurd zijn naar de drukcoëfficiënt, buigen scherp rond de zuigpiek en gaan dan over in een rustig zog, zonder zichtbare recirculatie in het domein.
De turbulente kinetische energie is in de vrije stroming verwaarloosbaar en vrijwel volledig geconcentreerd in de grenslaag en het nabije zog. Dat bevestigt dat de turbulentie in deze stroming plaatselijk door het profiel wordt opgewekt en niet vanaf de inlaat wordt meegevoerd, zodat de geschatte turbulentiewaarden aan de inlaat hier geen belangrijke bron van onzekerheid zijn.
