Version 2.0 -- cashocs : Un Logiciel de Optimisation de Forme et de Contrôle Optimal Basé sur la Méthode Adjointe

Informations Fondamentales

• ID de l'article : 2306.09828
• Titre : Version 2.0 -- cashocs: A Computational, Adjoint-Based Shape Optimization and Optimal Control Software
• Auteur : Sebastian Blauth (Fraunhofer ITWM, Kaiserslautern, Allemagne)
• Classification : math.OC (Optimisation et Contrôle)
• Date de soumission : Soumis à arXiv le 13 novembre 2023, publié dans SoftwareX 24
• Lien de l'article : https://arxiv.org/abs/2306.09828
• Version du logiciel : v2.0.12
• Licence : GNU GPL v3.0

Résumé

Cet article présente la mise à jour de la version 2.0 du logiciel cashocs. cashocs est un logiciel de calcul basé sur la méthode adjointe, spécialisé dans la résolution automatisée de problèmes d'optimisation contraints par des EDP, incluant l'optimisation de forme et le contrôle optimal. Depuis sa première publication, le logiciel a acquis de nombreuses fonctionnalités et outils pratiques, le rendant plus flexible et efficace. Les nouvelles fonctionnalités les plus importantes comprennent : le cadre de mappage spatial, la capacité de résolution de l'optimisation topologique basée sur la méthode des ensembles de niveaux, le support du calcul parallèle via MPI, et la capacité de traiter des contraintes supplémentaires (d'état).

Contexte de Recherche et Motivation

Définition du Problème

Les problèmes d'optimisation contraints par des EDP sont largement présents dans le calcul scientifique et l'ingénierie, comprenant principalement :

1. Problèmes d'optimisation de forme : optimisation de la performance du système en modifiant la forme géométrique
2. Problèmes de contrôle optimal : optimisation du comportement du système par des entrées de contrôle
3. Problèmes d'optimisation topologique : optimisation de la performance structurelle en modifiant la distribution de matière

Importance de la Recherche

• Ces problèmes ont une valeur d'application importante dans les domaines de l'aéronautique, de l'industrie automobile, de l'ingénierie biomédicale, etc.
• Les méthodes de résolution traditionnelles présentent une complexité de calcul élevée, nécessitant des techniques d'optimisation numérique spécialisées
• Les applications industrielles utilisent couramment des solveurs commerciaux, manquant de capacités d'optimisation directes

Limitations des Méthodes Existantes

• Absence d'un cadre logiciel unifié pour traiter différents types de problèmes d'optimisation contraints par des EDP
• Les outils existants présentent des insuffisances dans le calcul parallèle et le traitement des contraintes complexes
• L'efficacité de résolution des problèmes industriels complexes nécessite une amélioration

Contributions Principales

1. Cadre de Mappage Spatial : première implémentation complète de la technique de mappage spatial dans un logiciel d'optimisation contraint par des EDP, supportant les problèmes d'optimisation de forme et de contrôle optimal
2. Support de l'Optimisation Topologique : intégration d'un solveur d'optimisation topologique basé sur la méthode des ensembles de niveaux, incluant les dernières méthodes quasi-Newton
3. Capacité de Calcul Parallèle : support du calcul parallèle via MPI, permettant au logiciel de fonctionner sur des systèmes de calcul haute performance
4. Mécanisme de Traitement des Contraintes : ajout de la méthode de pénalité quadratique et de la méthode du lagrangien augmenté pour traiter les contraintes d'état et de contrôle supplémentaires
5. Améliorations de la Convivialité : incluant la mise à l'échelle automatique, les produits scalaires personnalisés, les flux de remaillage améliorés, etc.

Détails Méthodologiques

Conception de l'Architecture Logicielle

cashocs 2.0 est basé sur une conception modulaire, supportant trois grandes catégories de problèmes d'optimisation :

Problèmes d'Optimisation Contraints par des EDP
├── cashocs.ShapeOptimizationProblem
├── cashocs.OptimalControlProblem  
├── cashocs.TopologyOptimizationProblem
├── Contraintes Supplémentaires
│   ├── cashocs.ConstrainedShapeOptimizationProblem
│   └── cashocs.ConstrainedOptimalControlProblem
└── Problèmes de Mappage Spatial
    ├── cashocs.space_mapping.shape_optimization.SpaceMappingProblem
    └── cashocs.space_mapping.optimal_control.SpaceMappingProblem

Technique de Mappage Spatial

La méthode de mappage spatial exploite la structure hiérarchique des modèles, contenant un modèle fin (détaillé, complexe) et un modèle grossier (approximatif, peu coûteux en calcul) :

1. Hiérarchie des Modèles : optimisation efficace du modèle fin par optimisation continue et correction du modèle grossier approximatif
2. Éviter l'Optimisation Directe : pas besoin d'optimiser directement le modèle fin, particulièrement adapté aux scénarios industriels utilisant des solveurs commerciaux
3. Convergence Rapide : les expériences montrent une convergence en environ 5 itérations

Méthode d'Optimisation Topologique

Utilise la méthode des ensembles de niveaux pour résoudre les problèmes d'optimisation topologique :

1. Support Algorithmique : incluant les algorithmes d'optimisation topologique traditionnels et la nouvelle méthode quasi-Newton
2. Système Adjoint Automatique : bien que la dérivée topologique doive être fournie par l'utilisateur, la dérivation automatique du système adjoint a été implémentée
3. Élargissement du Champ d'Application : extension de l'applicabilité de cashocs au domaine de l'optimisation topologique avec informations de sensibilité topologique

Mécanisme de Traitement des Contraintes

Implémentation de deux principales méthodes de traitement des contraintes :

1. Méthode de Pénalité Quadratique : traitement des contraintes en ajoutant des termes de pénalité à la fonction objectif
2. Méthode du Lagrangien Augmenté : méthode hybride combinant multiplicateurs de Lagrange et fonction de pénalité
3. Applicabilité Générale : support de multiples types de contraintes incluant les contraintes d'état et de contrôle

Configuration Expérimentale

Pile Technologique

• Langage de Programmation : Python
• Dépendances Principales : FEniCS, NumPy, PETSc, MPI
• Traitement des Maillages : meshio, Gmsh
• Contrôle de Version : git
• Documentation : https://cashocs.readthedocs.io/

Configuration des Problèmes Exemples

Cas d'Optimisation de Forme par Mappage Spatial :

• Objectif : réaliser une distribution uniforme du débit à trois sorties de tuyauterie
• Modèle Fin : équations de Navier-Stokes à nombre de Reynolds 1000, résolu avec Ansys Fluent
• Modèle Grossier : système Stokes linéaire, résolu avec FEniCS sur un maillage grossier

Résultats Expérimentaux

Performance de la Méthode de Mappage Spatial

À partir des résultats de la Figure 1, on observe :

• Vitesse de Convergence : convergence en environ 5 itérations
• Efficacité de l'Optimisation : la géométrie finale réalise une distribution de débit quasi-parfaitement uniforme
• Efficacité Computationnelle : amélioration significative par rapport à l'optimisation directe du modèle fin

Performance du Calcul Parallèle

• Pas de Modification de Code : la plupart des problèmes d'optimisation contraints par des EDP en série peuvent être parallélisés directement
• Calcul Haute Performance : support de la résolution de problèmes d'optimisation à grande échelle sur des systèmes HPC
• Basé sur FEniCS : exploitation complète du support MPI de FEniCS, avec modifications minimales du code utilisateur

Améliorations du Remaillage

• Élargissement de l'Applicabilité : le nouveau flux de remaillage supporte les problèmes de mappage spatial et de traitement des contraintes
• Mise à Jour de la Syntaxe : les utilisateurs doivent suivre une syntaxe légèrement différente mais toujours intuitive
• Goulot d'Étranglement de Performance : le remaillage lui-même n'est pas parallélisé en raison des limitations de Gmsh et meshio

Travaux Connexes

Développement de la Technique de Mappage Spatial

• Des recherches antérieures existent sur la méthode de mappage spatial pour les problèmes de contrôle optimal
• La méthode de mappage spatial pour l'optimisation de forme est une technique récemment introduite
• cashocs est le premier logiciel à fournir un cadre de mappage spatial pour l'optimisation de forme et le contrôle optimal simultanément pour l'optimisation contraint par des EDP

Algorithmes d'Optimisation Topologique

• Intégration des derniers algorithmes d'optimisation topologique proposés par Amstutz et Andrä
• Implémentation de la méthode d'optimisation topologique quasi-Newton proposée par Blauth et Sturm
• Des améliorations restent possibles dans le calcul automatique des informations de sensibilité topologique

Conclusions et Discussion

Conclusions Principales

1. Complétude Fonctionnelle : cashocs 2.0 élargit significativement les fonctionnalités du logiciel, devenant une solution complète pour l'optimisation contraint par des EDP
2. Valeur d'Application Industrielle : le cadre de mappage spatial est particulièrement adapté aux applications industrielles, permettant une utilisation efficace des solveurs commerciaux existants
3. Amélioration des Capacités Computationnelles : le support parallèle MPI permet au logiciel de traiter des problèmes réels à grande échelle
4. Convivialité : plusieurs améliorations augmentent la facilité d'utilisation et la flexibilité du logiciel

Limitations

1. Calcul de la Dérivée Topologique : la dérivée topologique en optimisation topologique doit toujours être fournie manuellement par l'utilisateur
2. Goulot d'Étranglement du Remaillage : le processus de remaillage n'est pas parallélisé, pouvant devenir un goulot d'étranglement de performance pour certaines applications
3. Dépendance aux Outils : dépendance aux limitations des outils externes tels que Gmsh et meshio

Directions Futures

1. Dérivée Topologique Automatique : développement de méthodes pour calculer automatiquement la dérivée topologique
2. Remaillage Parallèle : amélioration de la performance parallèle du flux de remaillage
3. Algorithmes d'Optimisation Supplémentaires : intégration de plus d'algorithmes d'optimisation avancés et de méthodes de traitement des contraintes

Évaluation Approfondie

Avantages

1. Innovation Forte : première intégration dans un seul logiciel du mappage spatial, de l'optimisation topologique et du calcul parallèle
2. Valeur Pratique Élevée : particulièrement adapté aux applications industrielles, capable de s'intégrer aux logiciels commerciaux existants
3. Technologie Avancée : adoption des dernières techniques d'optimisation numérique et algorithmes
4. Documentation Complète : fourniture de tutoriels détaillés et de documentation complète
5. Disponibilité Open Source : la licence GNU GPL assure l'accessibilité du logiciel

Insuffisances

1. Dépendance Plateforme : dépendance élevée de l'écosystème FEniCS
2. Courbe d'Apprentissage : pour les utilisateurs non spécialisés, la maîtrise des concepts d'optimisation contraint par des EDP présente une certaine difficulté
3. Goulot d'Étranglement de Performance : la sérialisation du remaillage limite l'extensibilité de certaines applications

Impact

1. Contribution Académique : fourniture d'un outil important pour la recherche en optimisation contraint par des EDP
2. Valeur Industrielle : amélioration significative possible de l'efficacité de résolution des problèmes d'optimisation industriels
3. Construction Communautaire : la nature open source favorise la formation d'une communauté active d'utilisateurs et de développeurs

Scénarios d'Application

• Conception d'optimisation de forme dans l'industrie aéronautique
• Optimisation de dynamique des fluides dans l'industrie automobile
• Optimisation d'équipements en ingénierie biomédicale
• Optimisation topologique en ingénierie structurelle
• Projets d'optimisation nécessitant l'intégration avec des logiciels CFD commerciaux

Références Bibliographiques

Les références clés incluent :

1. Blauth, S. (2021). cashocs: A Computational, Adjoint-Based Shape Optimization and Optimal Control Software. SoftwareX, 13, 100646.
2. Blauth, S. (2023). Space Mapping for PDE Constrained Shape Optimization. SIAM J. Optim., 33(3), 1707–1733.
3. Amstutz, S., & Andrä, H. (2006). A new algorithm for topology optimization using a level-set method. J. Comput. Phys., 216(2), 573–588.

Résumé : cashocs 2.0 représente une mise à jour logicielle importante dans le domaine de l'optimisation contraint par des EDP. En intégrant le mappage spatial, l'optimisation topologique et le calcul parallèle et autres technologies avancées, le logiciel améliore significativement sa fonctionnalité et son utilité pratique. Ce logiciel est particulièrement adapté aux applications industrielles et académiques nécessitant une résolution efficace de problèmes d'optimisation complexes, constituant un outil important pour les chercheurs et ingénieurs du domaine.