2025-11-18T14:22:20.884401

Version 2.0 -- cashocs: A Computational, Adjoint-Based Shape Optimization and Optimal Control Software

Blauth
In this paper, we present version 2.0 of cashocs. Our software automates the solution of PDE constrained optimization problems for shape optimization and optimal control. Since its inception, many new features and useful tools have been added to cashocs, making it even more flexible and efficient. The most significant additions are a framework for space mapping, the ability to solve topology optimization problems with a level-set approach, the support for parallelism via MPI, and the ability to handle additional (state) constraints. In this software update, we describe the key additions to cashocs, which is now even better-suited for solving complex PDE constrained optimization problems.
academic

Version 2.0 -- cashocs: A Computational, Adjoint-Based Shape Optimization and Optimal Control Software

基本信息

  • 论文ID: 2306.09828
  • 标题: Version 2.0 -- cashocs: A Computational, Adjoint-Based Shape Optimization and Optimal Control Software
  • 作者: Sebastian Blauth (Fraunhofer ITWM, Kaiserslautern, Germany)
  • 分类: math.OC (Optimization and Control)
  • 发表时间: 2023年11月13日提交至arXiv，发表于SoftwareX 24
  • 论文链接: https://arxiv.org/abs/2306.09828
  • 软件版本: v2.0.12
  • 许可证: GNU GPL v3.0

摘要

本文介绍了cashocs软件的2.0版本更新。cashocs是一个基于伴随方法的计算软件，专门用于自动化求解PDE约束优化问题，包括形状优化和最优控制。自首次发布以来，该软件新增了许多功能和实用工具，使其更加灵活和高效。最重要的新增功能包括：空间映射框架、基于水平集方法的拓扑优化求解能力、通过MPI实现的并行计算支持，以及处理额外（状态）约束的能力。

研究背景与动机

问题定义

PDE约束优化问题在工程和科学计算中广泛存在，主要包括：

  1. 形状优化问题: 通过改变几何形状来优化系统性能
  2. 最优控制问题: 通过控制输入来优化系统行为
  3. 拓扑优化问题: 通过改变材料分布来优化结构性能

研究重要性

  • 这类问题在航空航天、汽车工业、生物医学工程等领域具有重要应用价值
  • 传统求解方法计算复杂度高，需要专业的数值优化技术
  • 工业应用中常使用商业求解器，缺乏直接的优化能力

现有方法局限性

  • 缺乏统一的软件框架处理不同类型的PDE约束优化问题
  • 现有工具在并行计算、复杂约束处理方面存在不足
  • 工业级复杂问题的求解效率有待提升

核心贡献

  1. 空间映射框架: 首次在PDE约束优化软件中实现了完整的空间映射技术，支持形状优化和最优控制问题
  2. 拓扑优化支持: 集成了基于水平集方法的拓扑优化求解器，包括最新的准牛顿方法
  3. 并行计算能力: 通过MPI实现并行计算支持，使软件能够在高性能计算系统上运行
  4. 约束处理机制: 新增二次罚函数法和增广拉格朗日方法来处理额外的状态和控制约束
  5. 用户友好性改进: 包括自动缩放、自定义标量积、改进的重网格化工作流程等

方法详解

软件架构设计

cashocs 2.0基于模块化设计，支持三大类优化问题：

PDE Constrained Optimization Problems
├── cashocs.ShapeOptimizationProblem
├── cashocs.OptimalControlProblem  
├── cashocs.TopologyOptimizationProblem
├── Additional Constraints
│   ├── cashocs.ConstrainedShapeOptimizationProblem
│   └── cashocs.ConstrainedOptimalControlProblem
└── Space Mapping Problems
    ├── cashocs.space_mapping.shape_optimization.SpaceMappingProblem
    └── cashocs.space_mapping.optimal_control.SpaceMappingProblem

空间映射技术

空间映射方法利用模型层次结构，包含精细模型（详细、复杂）和粗糙模型（近似、计算便宜）：

  1. 模型层次: 通过连续优化和修正粗糙模型近似来高效优化精细模型
  2. 避免直接优化: 无需直接优化精细模型，特别适合工业应用中使用商业求解器的场景
  3. 快速收敛: 实验显示约5次迭代即可达到收敛

拓扑优化方法

采用水平集方法求解拓扑优化问题：

  1. 算法支持: 包括传统拓扑优化算法和新颖的准牛顿方法
  2. 自动伴随系统: 虽然拓扑导数需要用户提供，但伴随系统的自动推导已实现
  3. 扩展应用范围: 将cashocs的适用性扩展到具有拓扑敏感性信息的拓扑优化领域

约束处理机制

实现了两种主要的约束处理方法：

  1. 二次罚函数法: 通过在目标函数中添加罚项来处理约束
  2. 增广拉格朗日方法: 结合拉格朗日乘数和罚函数的混合方法
  3. 广泛适用性: 支持状态约束、控制约束等多种约束类型

实验设置

技术栈

示例问题设置

空间映射形状优化案例:

  • 目标: 实现三个出口管道的均匀流量分布
  • 精细模型: 雷诺数为1000的Navier-Stokes方程，使用Ansys Fluent求解
  • 粗糙模型: 线性Stokes系统，在粗网格上使用FEniCS求解

实验结果

空间映射方法性能

从Figure 1的结果可以看出：

  • 收敛速度: 约5次迭代达到收敛
  • 优化效果: 最终几何形状实现了近乎完美的均匀流量分布
  • 计算效率: 相比直接优化精细模型显著提升效率

并行计算性能

  • 无代码修改: 大多数串行PDE约束优化问题可直接并行化
  • 高性能计算: 支持在HPC系统上求解大规模优化问题
  • 基于FEniCS: 充分利用FEniCS的MPI支持，用户端代码修改最小

重网格化改进

  • 扩展适用性: 新的重网格化工作流程支持空间映射和约束处理问题
  • 语法更新: 用户需要遵循略有不同但仍然直观的语法
  • 性能瓶颈: 由于Gmsh和meshio的限制，重网格化本身未并行化

相关工作

空间映射技术发展

  • 最优控制问题的空间映射方法已有相关研究
  • 形状优化的空间映射方法是最近才引入的新技术
  • cashocs是首个为PDE约束优化同时提供形状优化和最优控制空间映射框架的软件

拓扑优化算法

  • 集成了Amstutz和Andrä提出的最新拓扑优化算法
  • 实现了Blauth和Sturm提出的准牛顿拓扑优化方法
  • 在拓扑敏感性信息的自动计算方面仍有改进空间

结论与讨论

主要结论

  1. 功能完备性: cashocs 2.0显著扩展了软件功能，成为PDE约束优化的综合解决方案
  2. 工业应用价值: 空间映射框架特别适合工业应用，能够有效利用现有商业求解器
  3. 计算能力提升: MPI并行支持使软件能够处理大规模实际问题
  4. 用户友好性: 多项改进提升了软件的易用性和灵活性

局限性

  1. 拓扑导数计算: 拓扑优化中的拓扑导数仍需用户手动提供
  2. 重网格化瓶颈: 重网格化过程未并行化，可能成为某些应用的性能瓶颈
  3. 工具依赖: 依赖于Gmsh和meshio等外部工具的限制

未来方向

  1. 自动拓扑导数: 开发自动计算拓扑导数的方法
  2. 并行重网格化: 改进重网格化工作流程的并行性能
  3. 更多优化算法: 集成更多先进的优化算法和约束处理方法

深度评价

优点

  1. 创新性强: 首次在单一软件中集成空间映射、拓扑优化和并行计算
  2. 实用价值高: 特别适合工业应用，能够与现有商业软件集成
  3. 技术先进: 采用最新的数值优化技术和算法
  4. 文档完善: 提供详细的教程和文档支持
  5. 开源可用: GNU GPL许可证确保了软件的可访问性

不足

  1. 平台依赖: 高度依赖FEniCS生态系统
  2. 学习曲线: 对于非专业用户，掌握PDE约束优化概念仍有一定难度
  3. 性能瓶颈: 重网格化的串行化限制了某些应用的扩展性

影响力

  1. 学术贡献: 为PDE约束优化研究提供了重要工具
  2. 工业价值: 能够显著提升工业优化问题的求解效率
  3. 社区建设: 开源性质有利于形成活跃的用户和开发者社区

适用场景

  • 航空航天工业的形状优化设计
  • 汽车工业的流体动力学优化
  • 生物医学工程的设备优化
  • 结构工程的拓扑优化
  • 需要与商业CFD软件集成的优化项目

参考文献

关键参考文献包括：

  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.

总结: cashocs 2.0是PDE约束优化领域的一个重要软件更新，通过集成空间映射、拓扑优化和并行计算等先进技术，显著提升了软件的功能性和实用性。该软件特别适合需要高效求解复杂优化问题的工业和学术应用，是该领域研究者和工程师的重要工具。