💻 computer science

Adaptation of the hybrid fictitious domain-immersed boundary method for Reynolds-averaged turbulence modeling

Dit artikel presenteert een open-source adaptatie van de hybride fictitious domain-immersed boundary methode binnen OpenFOAM om robuuste, remeshing-vrije Reynolds-averaged Navier-Stokes (RANS) simulaties voor vorm- en topologieoptimalisatie mogelijk te maken, waarbij een nauwkeurigheid wordt aangetoond die vergelijkbaar is met standaard body-fitted CFD over diverse benchmarks en stromingscondities.

Oorspronkelijke auteurs: Lucie Kubíčková, Martin Isoz

Gepubliceerd 2026-06-05

📖 5 min leestijd🧠 Diepgaand

CC BY 4.0

Oorspronkelijke auteurs: Lucie Kubíčková, Martin Isoz

Oorspronkelijk artikel gelicentieerd onder CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Dit is een AI-gegenereerde uitleg van het onderstaande artikel. Het is niet geschreven of goedgekeurd door de auteurs. Raadpleeg het oorspronkelijke artikel voor technische nauwkeurigheid. Lees de volledige disclaimer

Stel je voor dat je een ingenieur bent die probeert de perfecte vorm te ontwerpen voor een auto-onderdeel, een buis of een vleugel. Om dit te doen, moet je simuleren hoe lucht of water rondom het object stroomt. Traditioneel is dit als het proberen te tekenen van een perfecte kaart van een kronkelende rivier door een stuk papier uit te snijden dat precies de vorm van de rivier heeft. Elke keer dat je de vorm van de rivier verandert (zelfs een klein beetje), moet je een compleet nieuw stuk papier uitknippen en de hele kaart opnieuw tekenen. Dit is traag, duur en frustrerend.

Dit artikel introduceert een nieuwe, slimmere manier om deze simulatie uit te voeren, genaamd openHFDIBRANS. Zo werkt het, met behulp van eenvoudige analogieën:

1. De "Ghost Wall"-truc (Immersed Boundary)

In plaats van voor elke nieuwe vorm een nieuw stuk papier (mesh) uit te snijden, gebruiken de onderzoekers een gigantisch, eenvoudig raster (zoals een schaakbord) dat het hele gebied beslaat. Ze veranderen het raster niet. In plaats daarvan gebruiken ze een "magische stift" om de computer te vertellen: "Hé, deze specifie dich specifieke vierkantjes zijn vaste wanden, en deze zijn lege lucht."

Als je de vorm van het object wilt veranderen, verplaats je gewoon de "magische stift" op het raster. Je hoeft het hele raster niet opnieuw op te bouwen. Het is als spelen met een digitale zandbak waarbij je een kasteel direct kunt hervormen zonder de hele zandtafel opnieuw op te bouwen.

2. Het "Verkeersregelaar"-probleem (Turbulence)

Het lastige deel is dat vloeistoffen (zoals lucht of water) rommelig en turbulent worden, vooral nabij wanden. In de echte wereld beweegt de lucht vlak naast een wand heel langzaam, terwijl de lucht iets verder weg snel beweegt. Dit wordt de "grenslaag" (boundary layer) genoemd.

Normaal gesproken hebben computers een zeer gedetailleerde, hoog-resolutie kaart nodig vlak naast de wand om deze kleine details te kunnen zien. Maar omdat ons "magische stift"-raster eenvoudig is en de wand niet perfect volgt, kan de computer deze minuscule details niet zien.

De auteurs hebben dit opgelost door een "Verkeersregelaar" (Wall Functions) toe te voegen. In plaats van te proberen elke kleine werveling van lucht bij de wand te zien, gebruikt de Verkeersregelaar een reeks regels (wiskundige formules) om te raden wat de lucht doet op basis van hoe snel deze een stukje verderop beweegt. Het is als een politieagent bij een druk kruispunt die niet elke individuele auto hoeft te zien om te weten dat de verkeersstroom zwaar is; ze gebruiken gewoon een vuistregel om de doorstroming te beheren.

3. De "Koorddans"-uitdaging (Stability)

Wanneer je deze "Ghost Wall"-truc combineert met de "Verkeersregelaar"-regels, kunnen de zaken instabiel worden. Het is alsof je probeert een bezem op je vinger te balanceren terwijl iemand er steeds tegenaan duwt. De computer kan in de war raken, de getallen kunnen krankzinnig worden en de simulatie kan crashen.

De auteurs hebben een speciale "Veiligheidsgordel" in hun software gebouwd. Ze hebben ervoor gezorgd dat de computer de regels controleert voordat hij de volgende stap probeert te berekenen, en ze hebben een "correctiestap" toegevoegd aan het einde van elke berekening om ervoor te zorgen dat de getallen op het koord blijven. Dit zorgt ervoor dat de simulatie niet van de rand af valt, zelfs niet wanneer de vormen vreemd zijn of de snelheden zeer hoog zijn.

4. Heeft het gewerkt? (De Resultaten)

Het team heeft hun nieuwe hulpmiddel getest op verschillende klassieke uitdagingen:

Water in een buis: Ze testten dit van zeer langzame stroming tot extreem snelle stroming. De resultaten kwamen bijna perfect overeen met de traditionele, trage "snijd-het-papier"-methode.
Een trede in een kanaal: Ze simuleerden water dat over een plotselinge daling stroomt. Het nieuwe hulpmiddel voorspelde waar het water zou gaan wervelen en bezinken, net zo goed als de oude methode.
Een auto-vorm (Ahmed body) en een vliegtuigvleugel (NACA profiel): Ze testten hoe lucht rond deze vormen stroomt. Het nieuwe hulpmiddel voorspelde de weerstand (drag) en lift (opwaartse kracht) zeer nauwkeurig.

De adder onder het gras:
In sommige gevallen heeft het nieuwe hulpmiddel de grootte van het zog (wake) iets "overgeschat" (de rommelige, wervelende lucht die achter een object achterblijft). Het is alsof het hulpmiddel dacht dat het object iets groter was dan het in werkelijkheid was, waardoor er een iets groter spoor van turbulentie achterbleef. De auteurs zeggen echter dat dit een voorspelbaar gedrag is en de resultaten niet verpest.

Waarom is dit belangrijk?

Het hoofddoel van dit werk is Optimalisatie. Omdat dit nieuwe hulpmiddel geen rooster vereist om telkens opnieuw opgebouwd te worden wanneer de vorm verandert, is het perfect voor Topologische Optimalisatie. Dit is een proces waarbij een computer automatisch duizenden verschillende vormen probeert uit te vinden om de absoluut beste vorm te vinden (bijvoorbeeld de lichtste brug of het meest aerodynamische auto-onderdeel).

De auteurs hebben hun code open-source gemaakt (gratis voor iedereen om te gebruiken) en gebouwd op een populair platform genaamd OpenFOAM. Ze beweren dat het robuust genoeg is om direct te worden gebruikt voor het ontwerpen van echte technische componenten, zoals sensoren of mallen voor het gieten van metaal, zonder dat je een wiskundig genie hoeft te zijn om het op te zetten.

Kortom: Ze hebben een snellere, flexibelere manier gebouwd om de stroming van vloeistoffen te simuleren die uitstekend werkt voor het ontwerpen van vormen, mits je een kleine, voorspelbare fout in de "staart" van de stroming accepteert.

Technische Samenvatting: Adaptatie van de Hybride Fictieve Domein-Immersed Boundary Methode voor RANS-modellering

Probleemstelling
Technische toepassingen vereisen vaak vorm- of topologieoptimalisatie (TO) van fluïdumcomponenten, waarbij Computational Fluid Dynamics (CFD) dient als basis voor de evaluatie van de kostenfunctie. Een belangrijke flessenhals in op CFD gebaseerde TO-frameworks is de noodzaak van frequente remeshing (hervermeshing) naarmate geometrieën veranderen. Hoewel Immersed Boundary Methods (IBM's) een oplossing bieden door eenvoudige, niet-conforme meshes te gebruiken en geometrie te representeren via scalaire indicatorvelden, blijft hun integratie met Reynolds-averaged Navier-Stokes (RANS) turbulentiemodellen uitdagend. De primaire moeilijkheden omvat het modelleren van de grenslaag zonder anisotrope meshverfijning, het voorkomen van spookoscillaties in druk en wandwandschaarspanning bij hoge Reynoldsgetallen, en het waarborgen van een nauwe koppeling tussen momentum- en turbulentievariabelen om instabiliteit te voorkomen. Bovendien zijn bestaande IBM-RANS solvers vaak niet open-source of ontbreken de standaard CFD-praktijken.

Methodologie
De auteurs stellen openHFDIBRANS voor, een open-source solver gebouwd op de OpenFOAM-bibliotheek die de Hybride Fictieve Domein-Immersed Boundary (HFDIB) methode aanpast voor stationaire, incompressibele RANS-simulaties. De methodologie integreert de volgende kerncomponenten:

Governing Equations: De solver behandelt de stationaire RANS-vergelijkingen gekoppeld aan vier twee-vergelijkingen turbulentiemodellen: $k$ - $\epsilon$ , $k$ - $\omega$ , $k$ - $\omega$ SST, en realizable $k$ - $\epsilon$ . De vergelijkingen worden gemodificeerd door een fase-indicatorveld $\lambda$ om oplossingen te beperken tot het fluïdumdomein en bevatten brontermen ( $f_{ib}$ en $S_{ib}$ ) om randvoorwaarden af te dwingen.
Masking en Celclassificatie: Het computationele domein wordt verdeeld in in-solid cellen ( $\Omega^h_{ic}$ ), boundary cellen ( $\Omega^h_{bc}$ ) en free-stream cellen ( $\Omega^h_{fc}$ ) op basis van het $\lambda$ -veld. Maskeringsvelden ( $\alpha_u, \alpha_k, \alpha_\omega$ ) beperken de toepassing van brontermen tot relevante regio's.
Wandmodellering en Randbehandeling: Om de grenslaag te behandelen zonder deze op te lossen, gebruikt de methode wandfuncties.
- Wandsnelheid (Friction Velocity): Om stabiliteit en een nauwe koppeling tussen momentum en turbulentie te garanderen, wordt de wandsnelheid ( $u_\tau$ ) berekend met behulp van turbulentievariabelen in plaats van lokale snelheidsprofielen.
- Imposed Fields: In boundary cellen waar de grenslaag niet is opgelost ( $y^+_{ib} > y^+_{lam}$ ), wordt de snelheid niet gereconstrueerd via interpolatie, maar gecorrigeerd uitsluitend door turbulente viscositeit. In opgeloste regio's ( $y^+_{ib} \leq y^+_{lam}$ ) wordt polynoominterpolatie gebruikt om snelheid en turbulentie kinetische energie ( $k$ ) op te leggen op basis van wandfuncties.
- Turbulentievariabelen: Er worden geen specifieke brontermen toegevoegd voor de specifieke dissipatiesnelheid ( $\omega$ ), omdat deze nabij wanden divergeert. In plaats daarvan worden $\omega$ en de productie van $k$ ( $P_k$ ) direct toegewezen in de getroffen cellen door manipulatie van de discretisatie-matrix, wat effectief Dirichlet-achtige randvoorwaarden afdwingt.
Oplossingsstrategie: De solver maakt gebruik van een uitgebreid SIMPLE-algoritme. Cruciaal is dat vóór de momentum-predictor de opgelegde velden worden berekend. Nadat de momentumvergelijking is opgelost, wordt een correctiestap ( $u = u + \alpha_u(u_{ib} - u)$ ) toegepast om de randvoorwaarden direct af te dwingen. Deze stap is essentieel voor convergentie in stationaire simulaties waar onderrelaxatie vereist is. De incompressibiliteitsvoorwaarde wordt alleen in het fluïdumdomein afgedwongen om conflicten met de IBM-forcing te vermijden.

Belangrijkste Bijdragen

Open-Source Implementatie: Het werk levert een volledig open-source implementatie van een IBM-RANS solver binnen het breed gebruikte OpenFOAM-framework, inclusief standaard CFD-praktijken.
Stabiliteitsgeoriënteerd Ontwerp: De auteurs pakken veelvoorkomende IBM-RANS instabiliteiten aan door de wandsnelheid te evalueren via turbulentievariabelen en een nauwe verbinding tussen momentum- en turbulentietransportvergelijkingen af te dwingen.
Robuustheid voor Optimalisatie: De methode is specifiek ontworate voor geometrieoptimalisatie-toepassingen en demonstreert robuustheid over diverse Reynoldsgetallen ( $10^1$ tot $10^6$ ) en complexe geometrieën zonder de noodzaak van remeshing.
Uitgebreide Benchmarking: De solver wordt gevalideerd tegen standaard body-fitted simulaties (simpleFoam) en experimentele data over zes scenario's: 2D/3D buisstroom, backward-facing step, stroming over een cilinder, NACA-0009 profiel en de Ahmed body.

Resultaten
De prestaties van openHFDIBRANS werden geëvalueerd tegen standaard body-fitted CFD (simpleFoam) en experimentele data:

Nauwkeurigheid: De solver produceerde resultaten die consistent waren met simpleFoam over een breed scala aan Reynoldsgetallen en verschillende turbulentiemodellen. $L_2$ -normfouten voor snelheid, turbulente viscositeit en Reynolds-spanning waren over het algemeen laag, hoewel iets hoger in regio's met abrupte veranderingen in turbulente viscositeit (bijv. $k$ - $\omega$ SST en realizable $k$ - $\epsilon$ modellen).
Gevoeligheid: De methode vertoonde een bevredigende robuustheid met betrekking tot de positie van de vaste-vloeistof interface ten opzichte van de mesh. Hoewel uitlijning met mesh-vlakken de nauwkeurigheid verbeterde, bleef de solver functioneel, zelfs met niet-uitgelijnde interfaces.
Weerstand en Lift: Voor de cilinder en het NACA-0009 profiel legde de solver stromingspatronen vast en voorspelde het lift- en weerstandscoëfficiënten met goede overeenstemming met referentiedata.
Beperkingen: In specifieke gevallen (bijv. backward-facing step en NACA-0009) had de solver de neiging om de grootte van het zog (wake) achter de vaste objecten licht te overschatten. Daarnaast kon het, net als andere RAS-modellen, de drag crisis voor de cilinder niet oplossen, wat geavanceerdere twee-laags modellen vereist.

Betekenis en Claims
Het artikel claimt dat openHFDIBRANS een levensvatbaar instrument is voor engineeringtoepassingen die geometrie- of topologieoptimalisatie vereisen. De auteurs benadrukken dat, hoewel de methode mogelijk een lichte overschatting van de grootte van het zog vertoont, het gedrag voorspelbaar en consistent is over verschillende geometrieën. Deze consistentie wordt geïdentificeerd als de meest kritieke factor voor optimalisatieframeworks. Het werk positioneert de solver als een direct inzetbaar alternatief voor standaard body-fitted methoden voor TO, waardoor de computationele kosten van frequente remeshing worden geëlimineerd terwijl voldoende nauwkeurigheid voor engineering-ontwerp behouden blijft. Toekomstige toepassingen die door de auteurs worden genoemd, zijn de topologieoptimalisatie van biosensoren en gietmatrijzen voor gietvormen.