Topology optimization of passively moving rigid bodies in unsteady flows

Cet article propose une méthode d'optimisation topologique pour la conception de corps rigides à mouvement passif dans des écoulements fluides instationnaires en couplant la dynamique des corps rigides avec des équations de fluide cinétique sur réseau et en appliquant une analyse de sensibilité basée sur l'adjoint pour optimiser les formes pour les mouvements de translation et de rotation.

L'essentiel

Imaginez que vous essayiez de concevoir la voile parfaite pour un bateau ou la pale parfaite pour une éolienne. Habituellement, les ingénieurs commencent par une forme, la placent dans le vent, observent comment elle se déplace, ajustent la forme, et recommencent. Ce document présente une méthode automatisée et ultra-intelligente pour y parvenir, mais avec une nuance : au lieu que la forme reste simplement là pendant que le vent souffle sur elle, la forme est autorisée à bouger parce que le vent la pousse.

Voici une décomposition de la façon dont cette « magie » fonctionne, en utilisant des analogies simples :

1. L'idée principale : Le danseur « passif »

La plupart des programmes informatiques qui conçoivent des formes supposent que l'objet est collé au sol (comme un pont ou un tuyau stationnaire). Si vous voulez concevoir une pièce mobile, comme une pale de ventilateur, vous dites généralement à l'ordinateur : « Faites tourner cette pale à 100 rotations par minute », et l'ordinateur calcule l'écoulement de l'air.

Ce papier inverse la logique. Il traite l'objet comme un danseur sur une piste de danse.

• L'ancienne méthode (Active) : Vous dites exactement au danseur comment bouger, et vous observez comment l'air se déplace autour de lui.
• La nouvelle méthode (Passive) : Vous ne dites pas au danseur comment bouger. Vous installez simplement la musique (le vent) et vous demandez à l'ordinateur de concevoir le corps du danseur de manière à ce que la musique le pousse naturellement à tournoyer ou à glisser le plus possible. Le mouvement du danseur est un résultat du vent, et non une commande.

2. L'astuce des deux grilles : La carte vs Le terrain

Pour que cela fonctionne, l'ordinateur utilise une astuce ingénieuse appelée « grilles séparées ». Imaginez que vous dessinez la carte d'une île mobile sur une feuille de papier millimétré.

• La grille de conception (La Carte) : C'est là que la forme est dessinée. C'est comme un carnet de croquis. L'ordinateur décide où se trouve le matériau « solide » (l'île) et l'espace « vide » (l'eau).
• La grille d'analyse (Le Terrain) : C'est là que la physique se produit. C'est une grille fixe d'eau et de vent.

À chaque infime fraction de seconde, la « Carte » (la forme) se déplace physémment et pivote. Ensuite, l'ordinateur projette cette carte mobile sur la grille du « Terrain » fixe pour calculer comment le vent pousse contre elle. Après la poussée du vent, l'ordinateur calcule comment l'objet doit bouger ensuite, met à jour la Carte, et répète le cycle. C'est comme prendre une photo d'un objet en mouvement, calculer la force du vent, déplacer l'objet, puis prendre la photo suivante instantanément.

3. La force « Fantôme » (Force de Brinkman)

Comment l'ordinateur sait-il où se trouve l'objet solide ? Il utilise un concept appelé force de Brinkman.
Considérez la zone de conception comme une pièce remplie de miel invisible et collant.

• Là où il y a du matériau solide, le miel est très épais et collant. Le vent ne peut pas passer à travers ; il pousse simplement contre la surface.
• Là où il y a de l'espace vide, le miel est fin ou inexistant, donc le vent circule librement.
  L'ordinateur n'a pas besoin de tracer une ligne dure ; il ajuste simplement la « viscosité » du miel en chaque point. Si la viscosité est élevée, c'est un mur ; si elle est faible, c'est de l'air. Cela permet à la forme de se transformer de manière fluide d'une forme à une autre.

4. Les mathématiques du « voyage dans le temps » (Méthode Adjointe)

Pour trouver la forme parfaite, l'ordinateur doit savoir : « Si je change ce minuscule point de matière ici, de combien la forme de l'objet sera-t-elle améliorée ? »
Calculer cela pour chaque point prendrait une éternité. C'est pourquoi les auteurs utilisent une méthode appelée Méthode de la Variable Adjointe.

• L'analogie : Imaginez que vous essayiez de trouver le meilleur chemin pour monter une montagne dans l'obscité. Au lieu de tester chaque chemin possible vers l'avant pour voir lequel est le meilleur, vous éclairez vers l'arrière depuis le sommet avec une lampe de poche. La lumière vous montre exactement quels pas mènent le plus efficacement vers le haut.
• Dans ce papier, la « lampe de poche » remonte le temps, calculant comment les forces du vent et le mouvement de l'objet auraient réagi à chaque infime changement de la forme. Cela donne à l'ordinateur une « carte de sensibilité » lui indiquant exactement où ajouter ou retirer du matériau pour obtenir le meilleur résultat.

5. Les résultats : Qu'ont-ils construit ?

L'équipe a testé cela sur trois scénarios :

1. La voile 2D : Ils ont conçu une forme qui part de l'immobilité et qui est poussée par le vent pour glisser horizontalement. Le résultat ressemblait à une aile d'avion incurvée (profil aérodynamique). Le vent poussait plus fort sur le dessus que sur le dessous, créant une portance qui tirait l'objet vers l'avant.
2. La turbine 2D : Ils ont conçu une forme qui tourne. Le résultat ressemblait à un hélice à quatre pales. Le vent frappait les pales courbes, créant une torsion qui la faisait pivoter.
3. La turbine 3D : Ils ont fait la même chose en 3D. Le résultat ressemblait à une éolienne du monde réel.

6. Le problème de la « nuance de gris »

Dans ces conceptions informatiques, les bords des formes ne sont pas toujours des lignes parfaitement nettes en noir et blanc. Parfois, ils sont en « nuances de gris » — un mélange de solide et d'air.

• Dans les exemples 2D, les auteurs ont constaté que même si l'on rendait la forme parfaitement nette (noir et blanc), la performance était presque identique. Les bords « flous » n'ont pas nui au résultat.
• Dans l'exemple 3D, les bords « flous » comptaient davantage. Parce que la grille informatique était un peu « grossière » (basse résolution), les bords flous modifiaient la façon dont le vent frappait les pales. Cela suggère que pour des formes 3D complexes, nous avons besoin d'une « carte » plus fine pour obtenir un résultat parfait.

Résumé

Ce papier présente une nouvelle façon de concevoir des machines mobiles (comme des voiles ou des turbines) où l'ordinateur détermine simultanément la forme et le mouvement. Il traite l'objet comme un danseur passif poussé par le vent, utilise l'astuce du « miel collant » pour définir la forme, et exécute une simulation mathématique remontant le temps pour trouver la forme la plus efficace. Le résultat sont des formes qui ressemblent naturellement à des ailes et des hélices, optimisées pour se déplacer le plus loin possible ou tourner le plus vite possible sous l'effet des forces fluides.

Résumé technique

Résumé technique : Optimisation topologique de corps rigides en mouvement passif dans des écoulements instationnaires

Énoncé du problème
Cette étude traite de l'optimisation topologique de corps rigides qui sont entraînés par des forces induites par l'écoulement du fluide, tels que des voiles et des turbines. Contrairement à l'optimisation topologique traditionnelle en dynamique des fluides où le mouvement du corps rigide est prescrit (actif) pour piloter l'écoulement, ou lorsque des particules sont suivies dans un couplage unidirectionnel où elles n'affectent pas le fluide, ce travail se concentre sur un système à couplage bidirectionnel. Ici, le mouvement du corps rigide est passif, ce qui signifie que sa trajectoire et sa rotation sont entièrement déterminées par les forces fluides agissant sur lui. Le défi majeur consiste à résoudre simultanément les équations de l'écoulement fluide instationnaire avec les équations du mouvement du corps rigide, où les forces fluides servent d'entrées aux équations de mouvement, et la position du corps sert de condition limite pour le fluide.

Méthodologie
Le cadre proposé intègre plusieurs techniques numériques pour gérer l'interaction complexe entre le solide en mouvement et le fluide instationnaire :

• Équations de base : L'écoulement du fluide est modélisé par les équations de Navier-Stokes incompressibles, résolues via le Schéma Cinétique de Réseau (LKS), une version étendue de la méthode Lattice Boltzmann (LBM) adaptée aux écoulements instationnaires. La dynamique du corps rigide est régie par les équations de Newton-Euler, où les forces et les couples sont calculés comme la réaction de la force de Brinkman au sein du domaine fluide.
• Grilles de conception et d'analyse séparées : Un composant central de la méthodologie est la séparation de la grille de conception (où la forme est définie) et de la grille d'analyse (où les variables d'état du fluide sont calculées). La grille de conception subit un mouvement de corps rigide à chaque pas de temps et est projetée sur la grille d'analyse à l'aide d'une fonction de pondération. Cela permet la représentation directe du mouvement du corps rigide sans nécessiter de remeshage du domaine fluide.
• Analyse de sensibilité adjointe : Pour calculer efficacement les sensibilités de conception, l'étude emploie le Schéma Cinétique de Réseau Adjoint (ALKS). Les auteurs dérivent des équations adjointes continues qui couplent les équations de l'écoulement fluide avec les équations du mouvement du corps rigide. Cette formulation traite les variables de mouvement du corps rigide (position, vitesse, angle de rotation et vitesse angulaire) comme des champs d'état supplémentaires au sein du système adjoint.
• Algorithme d'optimisation : La forme est mise à jour à l'aide de la Méthode des Asymptotes Mobiles (MMA). Pour garantir des conceptions binaires et fabricables, un filtre de densité avec projection de Heaviside est utilisé. Un schéma de continuation est appliqué au paramètre de convexité et au paramètre de raideur de la projection pour guider l'optimisation d'un état initial diffus vers une distribution solide-fluide nette.
• Approximation périodique : Pour les problèmes impliquant un mouvement cyclique (ex: turbines), une approximation de problème périodique est employée. Les valeurs initiales des champs d'état et adjoints au début d'une itération d'optimisation sont fixées selon les valeurs terminales de l'itération précédente, réduisant ainsi le coût de calcul en évitant la simulation des phases de démarrage transitoires pour chaque itération.

Contributions clés

1. Formulation du mouvement passif : L'article présente un nouveau cadre d'optimisation topologique spécifiquement conçu pour les corps rigides entraînés par les forces fluides, se distinguant du mouvement actif ou du suivi de particules en couplage unidirectionnel.
2. Dérivation de l'adjoint couplé : Il fournit une dérivation rigoureuse des équations adjointes et des sensibilités de conception qui couplent pleinement la dynamique des fluides avec les équations du mouvement du corps rigide, en traitant les variables cinématiques du corps comme des champs d'état.
3. Mise en œuvre de la séparation de grille : L'étude démontre l'efficacité de l'approche de séparation de grille pour les problèmes de corps rigides en mouvement instationnaire, facilitant le transfert de quantités telles que le coefficient de Brinkman et la vitesse du solide entre les grilles mobiles et stationnaires.
4. Validation numérique : La méthode est vérifiée par rapport à des approximations par différences finies pour les cas de translation et de rotation, confirmant l'exactitude des sensibilités dérivées.

Résultats
La méthode a été appliquée à trois exemples numériques :

• Voile 2D : L'optimisation visait à maximiser le déplacement horizontal d'un corps rigide. La forme résultante ressemble à un profil aérodynamique, générant des forces de portance qui propulsent le corps vers l'avant. La forme optimisée surpasse significativement une forme de référence présentant une surface supérieure plane.
• Turbine 2D : L'objectif était de maximiser l'angle de rotation. En employant une condition périodique cyclique, l'optimisation a produit une forme dotée de quatre pales de type profil aérodynamique qui tournent efficacement sous l'effet des forces fluides. Les résultats indiquent que la présence de régions de niveaux de gris dans la forme optimisée a un effet négligeable sur la performance après binarisation.
• Turbine 3D : Une forme rotative tridimensionnelle a été optimisée, produisant une structure de type turbine. Cependant, en raison des limitations de calcul résultant d'une résolution de maillage plus grossière, la version binarisée présentait des différences non négligeables de phase et de période par rapport à la forme optimisée, soulignant la sensibilité des structures fines 3D à la résolution du maillage.

Signification et affirmations
Les auteurs affirment que la méthode proposée optimise efficacement les formes de dispositifs tels que des voiles et des turbines dont la performance est dictée par l'interaction entre les forces fluides et le mouvement du corps rigide. L'étude démontre qu'une analyse instationnaire est nécessaire pour de tels problèmes, car l'angle d'attaque effectif et les forces changent dynamiquement avec le mouvement du corps.

L'article note modestement que bien que la méthode soit efficace, des limitations existent, particulièrement dans les cas tridimensionnels où les contraintes de résolution de maillage peuvent conduire à des écarts entre les formes optimisées et binarisées. Les auteurs suggèrent que les travaux futurs devraient se concentrer sur l'amélioration de la résolution du maillage (par exemple, via des méthodes adaptatives) et sur l'extension du cadre aux écoulements turbulents. L'étude conclut que l'approche d'optimisation topologique couplée proposée constitue un outil robuste pour la conception de composants mobiles pilotés par les fluides, offrant des perspectives sur les mécanismes physiques de génération de portance et de traînée dans les régimes instationnaires.