OpenMDAO技术咨询：NLBGS优化阶段失效但分析阶段正常的原因与解决办法

我基于OpenMDAO开展流固耦合优化工作，流体求解器采用DaFoam，结构求解器采用TACS。单独进行耦合分析时，使用NLBGS（非线性块高斯-塞德尔）方法可快速完成收敛；但在采用IPOPT优化器驱动优化流程时，即使是相同设计点，耦合求解的难度骤增，经常无法满足收敛准则。我的案例类似DaFoam官方的翼型流固耦合优化教程。

可能的原因分析 #

• 优化器迭代的步长与非线性区域影响：IPOPT在优化过程中会持续调整设计变量，即便初始点是已收敛的分析点，优化器的迭代步长可能将耦合系统推入非预期的强非线性区域。NLBGS对初始猜测的敏感度极高，原分析模式下的收敛策略在这类区域会失效。
• 导数计算的数值噪声干扰：IPOPT需要计算目标函数、约束对设计变量的导数，OpenMDAO的自动微分（AD）或有限差分（FD）过程可能引入数值噪声，导致耦合系统的残差计算出现微小扰动，破坏NLBGS的收敛稳定性。
• 收敛准则的设置差异：分析模式下的耦合收敛准则通常较宽松，而优化模式下为保证导数精度，往往需要更严格的准则，但NLBGS在严格准则下的收敛速度会大幅下降，甚至出现发散。
• 设计变量的耦合效应放大：优化时设计变量同时影响流体网格变形和结构刚度，两者的耦合效应被放大，而NLBGS作为块迭代方法，对强耦合问题的鲁棒性本身不足。

解决建议 #

• 调整NLBGS的迭代策略：
  • 增大NLBGS的最大迭代次数，将松弛因子从默认的1.0调整至0.8~0.9范围，避免迭代振荡。
  • 启用NLBGS的重启（restart）机制，当迭代停滞时重置迭代历史，提升鲁棒性。
• 优化导数计算方案：
  • 优先采用伴随导数（Adjoint）替代有限差分，减少数值噪声同时提升计算效率，需确保DaFoam和TACS的伴随求解器与OpenMDAO接口兼容配置正确。
  • 若使用有限差分，适当增大差分步长（需平衡精度与噪声），或改用中心差分替代向前差分，降低数值误差。
• 设置分层收敛准则：
  • 优化迭代早期采用宽松的耦合收敛准则加快进度，接近最优解时逐步收紧准则以保证导数精度，可通过OpenMDAO的回调函数或run_once机制实现。
  • 调整IPOPT的核心参数，如tol（容忍度）、max_iter（最大迭代次数），避免IPOPT因耦合残差未收敛而过早终止。
• 增强耦合系统鲁棒性：
  • 耦合迭代前加入预条件处理，对流体或结构的残差进行缩放，使各块残差量级相当，提升NLBGS的收敛速度。
  • 尝试切换耦合迭代方法，如用Newton-Krylov替代NLBGS，Newton法对强耦合问题收敛性更优，可借助OpenMDAO的自动微分计算耦合系统Jacobian矩阵。
• 规范设计变量配置：
  • 检查设计变量的边界约束与缩放设置，避免IPOPT产生过大步长导致系统进入极端非线性区域，可通过scaler对设计变量做归一化处理。

内容的提问来源于stack exchange，提问作者ja97