Version 2.0 -- cashocs: Eine rechnergestützte, adjungierte Formoptimierungs- und optimale Kontrollsoftware

Grundinformationen

• Papier-ID: 2306.09828
• Titel: Version 2.0 -- cashocs: A Computational, Adjoint-Based Shape Optimization and Optimal Control Software
• Autor: Sebastian Blauth (Fraunhofer ITWM, Kaiserslautern, Deutschland)
• Klassifizierung: math.OC (Optimierung und Kontrolle)
• Einreichungszeit: Eingereicht bei arXiv am 13. November 2023, veröffentlicht in SoftwareX 24
• Papierlink: https://arxiv.org/abs/2306.09828
• Softwareversion: v2.0.12
• Lizenz: GNU GPL v3.0

Zusammenfassung

Dieser Artikel stellt die Aktualisierung von Version 2.0 der cashocs-Software vor. cashocs ist eine rechnergestützte Software, die auf der Adjungiertenmethode basiert und speziell zur automatisierten Lösung von PDE-beschränkten Optimierungsproblemen entwickelt wurde, einschließlich Formoptimierung und optimaler Kontrolle. Seit der ersten Veröffentlichung hat die Software viele neue Funktionen und praktische Werkzeuge erhalten, die sie flexibler und effizienter machen. Die wichtigsten neuen Funktionen sind: ein Raumabbildungsrahmen, Topologieoptimierungslösungsfähigkeiten basierend auf der Level-Set-Methode, Unterstützung für parallele Berechnung durch MPI sowie die Fähigkeit, zusätzliche (Zustands-)Beschränkungen zu verarbeiten.

Forschungshintergrund und Motivation

Problemdefinition

PDE-beschränkte Optimierungsprobleme sind in der Ingenieur- und wissenschaftlichen Berechnung weit verbreitet und umfassen hauptsächlich:

1. Formoptimierungsprobleme: Optimierung der Systemleistung durch Änderung der geometrischen Form
2. Optimale Kontrollprobleme: Optimierung des Systemverhaltens durch Kontrolleingaben
3. Topologieoptimierungsprobleme: Optimierung der Strukturleistung durch Änderung der Materialverteilung

Forschungsbedeutung

• Diese Probleme haben wichtige Anwendungswerte in der Luft- und Raumfahrt, der Automobilindustrie, der biomedizinischen Technik und anderen Bereichen
• Traditionelle Lösungsmethoden haben hohe Rechenkomplexität und erfordern spezialisierte numerische Optimierungstechniken
• In industriellen Anwendungen werden häufig kommerzielle Solver verwendet, denen direkte Optimierungsfähigkeiten fehlen

Einschränkungen bestehender Methoden

• Mangel an einheitlichem Softwarerahmen zur Behandlung verschiedener Arten von PDE-beschränkten Optimierungsproblemen
• Bestehende Werkzeuge weisen Mängel bei paralleler Berechnung und komplexer Beschränkungsverarbeitung auf
• Die Lösungseffizienz für industrielle komplexe Probleme muss verbessert werden

Kernbeiträge

1. Raumabbildungsrahmen: Erstmalige Implementierung einer vollständigen Raumabbildungstechnik in PDE-beschränkter Optimierungssoftware mit Unterstützung für Formoptimierungs- und optimale Kontrollprobleme
2. Topologieoptimierungsunterstützung: Integration eines Topologieoptimierungslösers basierend auf der Level-Set-Methode, einschließlich neuester Quasi-Newton-Methoden
3. Parallelberechnungsfähigkeit: Implementierung von Parallelberechnungsunterstützung durch MPI, die es der Software ermöglicht, auf Hochleistungsrechensystemen zu laufen
4. Beschränkungsverarbeitungsmechanismus: Hinzufügung von quadratischer Strafmethode und erweiterter Lagrange-Methode zur Verarbeitung zusätzlicher Zustands- und Kontrollbeschränkungen
5. Verbesserungen der Benutzerfreundlichkeit: Einschließlich automatischer Skalierung, benutzerdefinierter Skalarprodukte, verbesserter Remeshing-Workflows usw.

Methodische Details

Softwarearchitektur-Design

cashocs 2.0 basiert auf modularem Design und unterstützt drei Hauptklassen von Optimierungsproblemen:

PDE-beschränkte Optimierungsprobleme
├── cashocs.ShapeOptimizationProblem
├── cashocs.OptimalControlProblem  
├── cashocs.TopologyOptimizationProblem
├── Zusätzliche Beschränkungen
│   ├── cashocs.ConstrainedShapeOptimizationProblem
│   └── cashocs.ConstrainedOptimalControlProblem
└── Raumabbildungsprobleme
    ├── cashocs.space_mapping.shape_optimization.SpaceMappingProblem
    └── cashocs.space_mapping.optimal_control.SpaceMappingProblem

Raumabbildungstechnik

Die Raumabbildungsmethode nutzt Modellhierarchien, die ein feines Modell (detailliert, komplex) und ein grobes Modell (approximativ, rechengünstig) enthalten:

1. Modellhierarchie: Effiziente Optimierung des feinen Modells durch kontinuierliche Optimierung und Korrektur der groben Modellapproximation
2. Vermeidung direkter Optimierung: Keine direkte Optimierung des feinen Modells erforderlich, besonders geeignet für industrielle Anwendungen mit kommerziellen Solvern
3. Schnelle Konvergenz: Experimente zeigen Konvergenz in etwa 5 Iterationen

Topologieoptimierungsmethode

Verwendung der Level-Set-Methode zur Lösung von Topologieoptimierungsproblemen:

1. Algorithmusunterstützung: Einschließlich traditioneller Topologieoptimierungsalgorithmen und neuartiger Quasi-Newton-Methoden
2. Automatisches Adjungierungssystem: Obwohl topologische Ableitungen vom Benutzer bereitgestellt werden müssen, wurde die automatische Ableitung des Adjungierungssystems implementiert
3. Erweiterte Anwendungsbereiche: Erweiterung der Anwendbarkeit von cashocs auf Topologieoptimierung mit topologischen Sensitivitätsinformationen

Beschränkungsverarbeitungsmechanismus

Implementierung von zwei Hauptmethoden zur Beschränkungsverarbeitung:

1. Quadratische Strafmethode: Verarbeitung von Beschränkungen durch Hinzufügung von Straftermen zur Zielfunktion
2. Erweiterte Lagrange-Methode: Hybridmethode, die Lagrange-Multiplikatoren und Strafmethoden kombiniert
3. Breite Anwendbarkeit: Unterstützung verschiedener Beschränkungstypen wie Zustandsbeschränkungen und Kontrollbeschränkungen

Experimentelle Einrichtung

Technologie-Stack

• Programmiersprache: Python
• Kernabhängigkeiten: FEniCS, NumPy, PETSc, MPI
• Netzverarbeitung: meshio, Gmsh
• Versionskontrolle: git
• Dokumentation: https://cashocs.readthedocs.io/

Beispielprobleme

Raumabbildungs-Formoptimierungsfall:

• Ziel: Erreichung einer gleichmäßigen Durchsatzverteilung an drei Auslassrohren
• Feines Modell: Navier-Stokes-Gleichungen bei Reynolds-Zahl 1000, gelöst mit Ansys Fluent
• Grobes Modell: Lineares Stokes-System auf grobem Netz, gelöst mit FEniCS

Experimentelle Ergebnisse

Leistung der Raumabbildungsmethode

Aus den Ergebnissen in Abbildung 1 ist ersichtlich:

• Konvergenzgeschwindigkeit: Konvergenz in etwa 5 Iterationen
• Optimierungseffekt: Die endgültige Geometrie erreicht eine nahezu perfekte gleichmäßige Durchsatzverteilung
• Recheneffizienz: Signifikante Verbesserung im Vergleich zur direkten Optimierung des feinen Modells

Parallelberechnungsleistung

• Keine Codeänderungen erforderlich: Die meisten seriellen PDE-beschränkten Optimierungsprobleme können direkt parallelisiert werden
• Hochleistungsrechnen: Unterstützung für die Lösung großer Optimierungsprobleme auf HPC-Systemen
• Basierend auf FEniCS: Vollständige Nutzung der MPI-Unterstützung von FEniCS mit minimalen Codeänderungen auf Benutzerseite

Remeshing-Verbesserungen

• Erweiterte Anwendbarkeit: Der neue Remeshing-Workflow unterstützt Raumabbildungs- und Beschränkungsverarbeitungsprobleme
• Syntaxaktualisierung: Benutzer müssen eine etwas andere, aber immer noch intuitive Syntax befolgen
• Leistungsengpass: Aufgrund von Einschränkungen in Gmsh und meshio ist das Remeshing selbst nicht parallelisiert

Verwandte Arbeiten

Entwicklung der Raumabbildungstechnik

• Raumabbildungsmethoden für optimale Kontrollprobleme sind bereits erforscht
• Raumabbildungsmethoden für Formoptimierung sind eine kürzlich eingeführte neue Technologie
• cashocs ist die erste Software, die sowohl einen Formoptimierungs- als auch einen optimalen Kontrollrahmen für Raumabbildung für PDE-beschränkte Optimierung bietet

Topologieoptimierungsalgorithmen

• Integration neuester Topologieoptimierungsalgorithmen von Amstutz und Andrä
• Implementierung der Quasi-Newton-Topologieoptimierungsmethode von Blauth und Sturm
• Verbesserungspotenzial bei der automatischen Berechnung topologischer Sensitivitätsinformationen

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

1. Funktionale Vollständigkeit: cashocs 2.0 erweitert die Softwarefunktionalität erheblich und wird zur umfassenden Lösung für PDE-beschränkte Optimierung
2. Industrieller Anwendungswert: Der Raumabbildungsrahmen ist besonders für industrielle Anwendungen geeignet und kann bestehende kommerzielle Solver effektiv nutzen
3. Verbesserte Rechenkapazität: MPI-Parallelunterstützung ermöglicht es der Software, große praktische Probleme zu bewältigen
4. Verbesserte Benutzerfreundlichkeit: Mehrere Verbesserungen erhöhen die Benutzerfreundlichkeit und Flexibilität der Software

Einschränkungen

1. Topologische Ableitungsberechnung: Topologische Ableitungen in der Topologieoptimierung müssen vom Benutzer manuell bereitgestellt werden
2. Remeshing-Engpass: Der Remeshing-Prozess ist nicht parallelisiert und kann ein Leistungsengpass für einige Anwendungen sein
3. Werkzeugabhängigkeiten: Abhängigkeit von externen Werkzeugen wie Gmsh und meshio

Zukünftige Richtungen

1. Automatische topologische Ableitungen: Entwicklung von Methoden zur automatischen Berechnung topologischer Ableitungen
2. Paralleles Remeshing: Verbesserung der Parallelleistung des Remeshing-Workflows
3. Weitere Optimierungsalgorithmen: Integration weiterer fortschrittlicher Optimierungsalgorithmen und Beschränkungsverarbeitungsmethoden

Tiefgreifende Bewertung

Stärken

1. Hohe Innovativität: Erstmalige Integration von Raumabbildung, Topologieoptimierung und paralleler Berechnung in einer einzigen Software
2. Hoher praktischer Wert: Besonders geeignet für industrielle Anwendungen und Integration mit bestehender kommerzieller Software
3. Fortgeschrittene Technologie: Verwendung neuester numerischer Optimierungstechniken und Algorithmen
4. Umfassende Dokumentation: Bereitstellung detaillierter Tutorials und Dokumentationsunterstützung
5. Open Source verfügbar: GNU GPL-Lizenz gewährleistet die Zugänglichkeit der Software

Schwächen

1. Plattformabhängigkeit: Hohe Abhängigkeit vom FEniCS-Ökosystem
2. Lernkurve: Für Nicht-Spezialisten ist die Beherrschung von PDE-beschränkten Optimierungskonzepten immer noch schwierig
3. Leistungsengpässe: Die Serialisierung des Remeshing begrenzt die Skalierbarkeit einiger Anwendungen

Einfluss

1. Akademischer Beitrag: Bereitstellung wichtiger Werkzeuge für die Forschung zu PDE-beschränkter Optimierung
2. Industrieller Wert: Kann die Lösungseffizienz industrieller Optimierungsprobleme erheblich verbessern
3. Gemeinschaftsaufbau: Die Open-Source-Natur fördert die Bildung einer aktiven Benutzer- und Entwicklergemeinde

Anwendungsszenarien

• Formoptimierungsdesign in der Luft- und Raumfahrtindustrie
• Strömungsdynamische Optimierung in der Automobilindustrie
• Geräteoptimierung in der biomedizinischen Technik
• Topologieoptimierung in der Strukturtechnik
• Optimierungsprojekte, die Integration mit kommerzieller CFD-Software erfordern

Referenzen

Wichtige Referenzen umfassen:

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.

Zusammenfassung: cashocs 2.0 ist eine wichtige Softwareaktualisierung im Bereich der PDE-beschränkten Optimierung. Durch die Integration fortschrittlicher Techniken wie Raumabbildung, Topologieoptimierung und parallele Berechnung verbessert die Software ihre Funktionalität und Praktikabilität erheblich. Diese Software ist besonders für industrielle und akademische Anwendungen geeignet, die komplexe Optimierungsprobleme effizient lösen müssen, und ist ein wichtiges Werkzeug für Forscher und Ingenieure in diesem Bereich.