Topology optimization of passively moving rigid bodies in unsteady flows

Dit artikel stelt een topologieoptimalisatiemethode voor voor het ontwerpen van passief bewegende starre lichamen in onstabiele stromingen door starre lichaamsdynamica te koppelen aan lattice-kinetische vloeistofvergelijkingen en een adjoint-gebaseerde gevoeligheidsanalyse toe te passen om vormen te optimaliseren voor translationele en rotationele bewegingen.

De kern

Stel je voor dat je het perfecte zeil voor een boot of de perfecte bladen voor een windturbine probeert te ontwerpen. Meestal beginnen ingenieurs met een vorm, plaatsen die in de wind, kijken hoe het beweegt, passen de vorm aan en proberen het opnieuw. Dit artikel introduceert een super slimme, geautomatiseerde manier om dat te doen, maar met een twist: in plaats van dat de vorm daar gewoon stilzit terwijl de wind eroverheen blaast, krijgt de vorm de vrijheid om te bewegen omdat de wind ertegenaan duwt.

Hier is een uitsplitsing van hoe deze "magie" werkt, met behulp van eenvoudige analogieën:

1. Het Grote Idee: De "Passieve" Danser

De meeste computerprogramma's die vormen ontwerpen, gaan ervan uit dat het object aan de vloer vastgelijmd zit (zoals een brug of een stationaire pijp). Als je een bewegend onderdeel wilt ontwerpen, zoals een ventilatorblad, vertel je de computer meestal: "Laat dit blad 100 rotaties per minuut draaien," en de computer berekent hoe de luchtstroom verloopt.

Dit artikel draait het scenario om. Het behandelt het object als een danser op een dansvloer.

• De Oude Manier (Actief): Je vertelt de danser precies hoe hij moet bewegen, en je kijkt hoe de lucht om hem heen beweegt.
• De Nieuwe Manier (Passief): Je vertelt de danser niet hoe hij moet bewegen. Je stelt alleen de muziek (de wind) in die wordt afgespeeld, en je vraagt de computer om het lichaam van de danser zo te ontwerpen dat de muziek hem vanzelf doet draaien of glijden zo ver mogelijk. De beweging van de danser is een gevolg van de wind, geen opdracht.

2. De Twee-Grid-Truc: De Kaart versus het Terrein

Om dit werkend te krijgen, gebruikt de computer een slimme truc genaamd "gescheiden grids" (separated grids). Stel je voor dat je een kaart tekent van een bewegend eiland op een stuk grafisch papier.

• Het Ontwerp-Grid (De Kaart): Dit is waar de vorm wordt getekend. Het is als een schetsblok. De computer bepaalt hier waar het "vaste" materiaal (het eiland) en de "lege" ruimte (het water) zich bevinden.
• Het Analyse-Grid (Het Terrein): Dit is waar de fysica plaatsvindt. Het is een vast grid van water en wind.

Elke fractie van een seconde beweegt en roteert de "Kaart" (de vorm) fysiek. Vervolgens projecteert de computer die bewegende kaart op het vaste "Terrein"-grid om te berekenen hoe de wind tegen het object duwt. Nadat de wind heeft geduwd, berekent de computer hoe het object zich vervolgens moet bewegen, werkt de "Kaart" bij en herhaalt de cyclus. Het is alsof je een foto maakt van een bewegend object, de windkracht berekent, het object verplaatst en direct de volgende foto maakt.

3. De "Geestkracht" (Brinkman-kracht)

Hoe weet de computer waar het vaste object zich bevindt? De computer gebruikt een concept genaamd de Brinkman-kracht.
Denk aan het ontwerpgebied als een kamer gevuld met onzichtbare, plakkerige honing.

• Waar vast materiaal is, is de honing super dik en plakkerig. De wind kan er niet doorheen bewegen; hij duwt alleen tegen het oppervlak.
• Waar lege ruimte is, is de honing dun of zelfs afwezig, waardoor de wind vrij kan stromen.
  De computer hoeft geen harde lijn te tekenen; hij past simpelweg de "plakkerigheid" van de honing aan elk punt aan. Als de plakkerigheid hoog is, is het een muur; als deze laag is, is het lucht. Hierdoor kan de vorm soepel transformeren van de ene naar de andere vorm.

4. De "Tijdreis"-wiskunde (Adjoint-methode)

Om de perfecte vorm te vinden, moet de computer weten: "Als ik dit minuscule stipje materiaal hier verander, hoeveel beter zal het object dan bewegen?"
Het berekenen hiervan voor elk afzonderlijk stipje zou eeuwig duren. Daarom gebruiken de auteurs een methode genaamd de Adjoint Variable Method.

• De Analogie: Stel je voor dat je in het donker de beste route naar de top van een berg probeert te vinden. In plaats van elke mogelijke route vooruit te bewandelen om te zien welke de beste is, schijn je een zaklamp achteruit vanaf de top. Het licht laat je precies zien welke stappen het meest efficiënt de berg op leiden.
• In dit artikel loopt de "zaklamp" achteruit door de tijd, waarbij wordt berekend hoe de windkrachten en de beweging van het object zouden hebben gereageerd op elke kleine verandering in de vorm. Dit geeft de computer een "gevoeligheidskaart" die precies aangeeft waar materiaal moet worden toegevoegd of verwijderd om het beste resultaat te behalen.

5. De Resultaten: Wat hebben ze gebouwd?

Het team heeft dit getest op drie scenario's:

1. De 2D-zeil: Ze ontwierpen een vorm die stilstaat en door de wind wordt weggeduwd om horizontaal te glijden. Het resultaat leek op een gebogen vliegtuigvleugel (airfoil). De wind duwde harder op de bovenkant dan op de onderkant, wat een liftkracht creëerde die het object naar voren trok.
2. De 2D-turbine: Ze ontwierpen een vorm die draait. Het resultaat leek op een vierbladige propeller. De wind raakte de gebogen bladen, wat een draaiing veroorzaakte die het object liet roteren.
3. De 3D-turbine: Ze deden hetzelfde in 3D. Het resultaat leek op een echte windturbine uit de praktijk.

6. Het "Grijswaarden"-probleem

In deze computerontwerpen zijn de randen van de vormen niet altijd perfecte zwart-wit lijnen. Soms zijn ze "grijswaarden" — een beetje van beide, vast en lucht.

• In de 2D-voorbeelden ontdekten de auteurs dat zelfs als ze de vorm perfect scherp maakten (zwart en wit), de prestaties bijna hetzelfde waren. De "vage" randen beïnvloedden het resultaat nauwelijks.
• In het 3D-voorbeeld waren de "vage" randen belangrijker. Omdat het computernetwerk een beetje "korrelig" was (lage resolutie), veranderde de manier waarop de wind de bladen raakte. Dit suggereert dat we voor complexe 3D-vormen een fijnere "kaart" nodig hebben om een perfect resultaat te krijgen.

Samenvatting

Dit artikel presenteert een nieuwe manier om bewegende machines (zoals zeilen of turbines) te ontwerpen, waarbij de computer zowel de vorm als de beweging gelijktijdig bepaalt. Het behandelt het object als een passieve danser die door de wind wordt voortgestuwd, gebruikt een "plakkerige honing"-truc om de vorm te definiëren, en voert een wiskundige simulatie achteruit in de tijd uit om de meest efficiënte vorm te vinden. Het resultaat zijn vormen die er natuurlijk uitzien als vleugels en propellers, geoptimaliseerd om zo ver mogelijk te glijden of zo snel mogelijk te draaien onder invloed van vloeistofkrachten.

Technische samenvatting

Technische Samenvatting: Topologieoptimalisatie van passief bewegende starre lichamen in onstabiele stromingen

Probleemstelling
Deze studie behandelt de topologieoptimalisatie van starre lichamen die worden aangedreven door door vloeistofstroming geïnduceerde krachten, zoals zeilen en turbines. In tegenstelling tot traditionele topologieoptimalisatie in de vloeistofdynamica, waarbij de beweging van het starre lichaam wordt voorgeschreven (actief), of waarbij deeltjes worden gevolgd in een eenrichtingskoppeling waarbij zij de stroming niet beïnvloeden, richt dit werk zich op een tweeweggekoppeld systeem. Hierbij is de beweging van het starre lichaam passief, wat betekent dat het traject en de rotatie volledig worden bepaald door de vloeistofkrachten die erop inwerken. De kernuitdaging bestaat uit het gelijktijdig oplossen van de onstabiele vloeistofstromingsvergelijkingen met de bewegingsvergelijkingen van het starre lichaam, waarbij de vloeistofkrachten dienen als inputs voor de bewegingsvergelijkingen en de positie van het lichaam dient als een randvoorwaarde voor de vloeistof.

Methodologie
Het voorgestelde kader integreert verschillende numerieke technieken om de complexe interactie tussen de bewegende vaste stof en de onstabiele vloeistof te beheersen:

• Governing Equations (Sturende vergelijkingen): De vloeistofstroom wordt gemodelleerd met behulp van de incompressibele Navier-Stokes-vergelijkingen, opgelost via het Lattice Kinetic Scheme (LKS), een uitgebreide versie van de Lattice Boltzmann Method (LBM) die geschikt is voor onstabiele stromingen. De dynamica van het starre lichaam wordt beheerst door de Newton-Euler-vergelijkingen, waarbij de krachten en momenten worden berekend als reactie van de Brinkman-kracht binnen het vloeistofdomein.
• Gescheiden ontwerp- en analysegrids: Een kernonderdeel van de methodologie is de scheiding van het ontwerpgrid (waar de vorm wordt gedefinieerd) en het analysegrid (waar de vloeistofvariabelen worden berekend). Het ontwerpgrid ondergaat in elke tijdstap een starre lichaamsbeweging en wordt op het analysegrid geprojecteerd met behulp van een wegingsfunctie. Dit maakt de directe representatie van de beweging van het starre lichaam mogelijk zonder het herberekenen van het mesh (remeshing) van het vloeistofdomein.
• Adjoint Sensitivity Analysis (Adjoint gevoeligheidsanalyse): Om de ontwerpgevoeligheden efficiënt te berekenen, maakt de studie gebruik van het Adjoint Lattice Kinetic Scheme (ALKS). De auteurs leiden continue adjoint-vergelijkingen af die de vloeistofstromingsvergelijkingen koppelen aan de bewegingsvergelijkingen van het starre lichaam. Deze formulering behandelt de kinematische variabelen van het starre lichaam (positie, snelheid, rotatiehoek en hoeksnelheid) als aanvullende toestandsvelden binnen het adjoint-systeem.
• Optimalisatiealgoritme: De vorm wordt bijgewerkt met behulp van de Method of Moving Asymptotes (MMA). Om produceerbare, binair-achtige ontwerpen te garanderen, wordt een dichtheidsfilter met Heaviside-projectie gebruikt. Een continuatieschema wordt toegepast op de convexiteitsparameter en de projectiessteilheid om het optimalisatieproces van een diffuse beginstaat naar een duidelijke vaste-stof-vloeistofverdeling te leiden.
• Periodieke Benadering: Voor problemen met cyclische beweging (bijv. turbines) wordt een periodieke probleembenadering gehanteerd. De beginwaarden voor de toestand en de adjoint-velden aan het begin van een optimalisatie-iteratie worden ingesteld op de eindwaarden van de vorige iteratie, wat de computationele kosten verlaagt door het simuleren van transiënte opstartfasen voor elke iteratie te vermijden.

Belangrijkste Bijdragen

1. Passieve Bewegingsformulering: Het artikel presenteert een nieuw topologieoptimalisatiekader specifief voor starre lichamen die worden aangedreven door vloeistofkrachten, wat het onderscheidt van actieve beweging of eenrichtingskoppeling bij deeltjes-tracking.
2. Gekoppelde Adjoint Afleiding: Het biedt een rigoureuze afleiding van de adjoint-vergelijkingen en ontwerpgevoeligheden die de vloeistofdynamica volledig koppelen aan de bewegingsvergelijkingen van het starre lichaam, waarbij de kinematische variabelen van het lichaam als toestandsvelden worden behandeld.
3. Implementatie van Grid-scheiding: De studie demonstreert de effectiviteit van de grid-scheidingsaanpak voor onstabiele, bewegende starre lichaamsproblemen, wat de overdracht van grootheden zoals de Brinkman-coëfficiënt en de snelheid van de vaste stof tussen bewegende en stationaire grids vergemakkelijkt.
4. Numerieke Validatie: De methode is geverifieerd tegen eindige-verschillenbenaderingen voor zowel translationele als rotationele gevallen, wat de nauwkeurigheid van de afgeleide gevoeligheden bevestigt.

Resultaten
De methode werd toegepast op drie numerieke voorbeelden:

• 2D Zeil: De optimalisatie maximaliseerde de horizontale verplaatsing van een star lichaam. De resulterende vorm leek op een profiel (airfoil) die liftkrachten genereert om het lichaam vooruit te stuwen. De geoptimaliseerde vorm presteerde aanzienlijk beter dan een referentievorm met een plat bovenoppervlak.
• 2D Turbine: Het doel was om de rotatiehoek te maximaliseren. Door een cyclische periodieke conditie toe te passen, produceerde de optimalisatie een vorm met vier profielachtige bladen die effectief roteren onder invloed van vloeistofkrachten. De resultaten gaven aan dat de aanwezigheid van grijswaardenregio's in de geoptimaliseerde vorm een verwaarloosbaar effect had op de prestaties wanneer deze gebinariseerd werd.
• 3D Turbine: Een driedimensionale roterende vorm werd geoptimaliseerd, wat een turbine-achtige structuur opleverde. Vanwege computationele beperkingen die leidden tot een grovere mesh-resolutie, vertoonde de gebinariseerde versie echter niet-verwaarloosbare verschillen in fase en periode vergeleken met de geoptimaliseerde vorm, wat de gevoeligheid van 3D dunne structuren voor mesh-resolutie benadrukt.

Betekenis en Claims
De auteurs stellen dat de voorgestelde methode effectief de vormen kan optimaliseren van apparaten zoals zeilen en turbines waarbij de prestaties worden bepaald door de interactie tussen vloeistofkrachten en de beweging van het starre lichaam. De studie toont aan dat een onstabiele analyse noodzakelijk is voor dergelijke problemen, aange omdat de effectieve invalshoek en de krachten dynamisch veranderen met de beweging van het lichaam.

Het artikel merkt bescheiden op dat hoewel de methode effectief is, er beperkingen bestaan, met name in driedimensionale gevallen waar mesh-resolutiebeperkingen kunnen leiden tot discrepanties tussen de geoptimaliseerde en de gebinariseerde vormen. De auteurs suggereren dat toekomstig werk zich moet richten op het verbeteren van de mesh-resolutie (bijv. via adaptieve methoden) en het uitbreiden van het kader naar turbulente stromingen. De studie concludeert dat de voorgestelde gekoppelde topologieoptimalisatie een robuust hulpmiddel biedt voor het ontwerpen van vloeistofgestuurde bewegende componenten, en inzichten biedt in de fysieke mechanismen van lift- en weerstandgeneratie in onstabiele regimes.