Elliptic Relaxation Strategies to Support Numerical Stability of Segregated Continuous Adjoint Flow Solvers

TL;DR

This paper presents an elliptic relaxation strategy to enhance the numerical stability of segregated continuous adjoint flow solvers, especially on unstructured grids, by smoothing and redistributing key coupling terms.

Contribution

The authors introduce a PDE-based elliptic relaxation method that stabilizes adjoint flow computations, addressing numerical instabilities caused by cross-coupling terms on complex grids.

Findings

Effective stabilization of adjoint simulations on unstructured grids.

Improved optimization results in shape design problems.

Method reduces divergence issues in practical CFD applications.

Abstract

This paper introduces a novel method for numerically stabilizing sequential continuous adjoint flow solvers utilizing an elliptic relaxation strategy. The proposed approach is formulated as a Partial Differential Equation (PDE) containing a single user-defined parameter, which analytical investigations reveal to represent the filter width of a probabilistic density function or Gaussian kernel. Key properties of the approach include (a) smoothing features with redistribution capabilities while (b) preserving integral properties. The technique targets explicit adjoint cross-coupling terms, such as the Adjoint Transpose Convection (ATC) term, which frequently causes numerical instabilities, especially on unstructured grids common in industrial applications. A trade-off is made by sacrificing sensitivity consistency to achieve enhanced numerical robustness. The method is validated on a two-phase, laminar, two-dimensional cylinder flow test case at Re=20 and Fn=0.75, focusing on minimizing resistance or maximizing lift. A range of homogeneous and inhomogeneous filter widths is evaluated. Subsequently, the relaxation method is employed to stabilize adjoint simulations during shape optimizations that aim at drag reduction of ship hulls. Two case studies are considered: A model-scale bulk carrier traveling at Re=7.246E+06 and Fn=0.142 as well as a harbor ferry cruising at Re=2.43E+08 and Fn=0.4 in full-scale conditions. Both cases, characterized by unstructured grids prone to adjoint divergence, demonstrate the effectiveness of the proposed method in overcoming stability challenges. The resulting optimizations achieve superior outcomes compared to approaches that omit problematic coupling terms.