Abstract:Adjoint methods efficiently compute gradients for systems with many inputs and have been widely used for large-scale gradient-based optimization in fluid mechanics. To ensure optimization's numerical robustness, we need to develop an adjoint solution algorithm that has a similar, if not the same, convergence rate as the primal flow solver at each optimization iteration. This consistent primal-adjoint convergence behavior is also called duality-preserving (DP). Existing DP adjoint methods are mostly for fully coupled Navier–Stokes (NS) solvers with compressible flows. However, few DP adjoints have been proposed for segregated NS solvers, e.g., the semi-implicit method for pressure-linked equations (SIMPLE) algorithm solvers widely used for incompressible flow simulations. This study will fill this gap by deriving a fixed-point segregated adjoint formulation that fully preserves the convergence behavior of primal solvers. We first rewrite the steady-state segregated NS primal solution process into a fully coupled left-preconditioned Richardson format. Then, we transpose the preconditioner matrix to obtain a DP iterative adjoint formulation. In addition to algebraic derivation, we create a new graph representation that significantly streamlines the DP adjoint development for new segregated solvers. We prove that the proposed graph representation is equivalent to primal and adjoint solutions and guarantees duality. We evaluate our proposed algorithm using three cases with increasing complexity: a miniature-sized problem, an airfoil, and a wing. To quantify duality, we compare the eigenvalues between the primal and adjoint solutions and observe excellent agreements. We also evaluate to what extent various numerical settings (e.g., inner iteration tolerances and non-dual preconditioner) impact the eigenvalue distributions. Our proposed adjoint achieves machine-precision accurate gradient computation with competitive speed. Finally, we incorporate our adjoint solver into a large-scale gradient-based optimization framework and demonstrate its capability for wing aerodynamic shape optimization. The proposed fixed-point adjoint approach has the potential to make the adjoint solution more robust and efficient for any segregated NS solvers.

Ordinary differential equations for the adjoint Euler equations

Adjoint and direct characteristic equations for two-dimensional compressible Euler flows

Numerical Simulation for the System of Ordinary Differential Equations

Analytic Adjoint Solutions for the 2D Incompressible Euler Equations Using the Green's Function Approach

Shock equations and jump conditions for the 2D Adjoint Euler equations

Geometric Methods for Adjoint Systems

An Adjoint Method for a High-Order Discretization of Deforming Domain Conservation Laws for Optimization of Flow Problems

Euler – Poincaré Equations in Geophysical Fluid Dynamics 3 1 Introduction

Fredholm determinants, continued fractions, Jost and Evans functions for a Jacobi matrix associated with the 2D-Euler equations

Understanding differential equations through diffusion point of view: non-symmetric discrete equations

A mathematical model for flow field with condensation discontinuity: One‐dimensional Euler equations with singular sources

Unsteady Dynamic Evolution Adjoint Design Method

Solution of singularly perturbed system of 2D parabolic convection diffusion differential-difference equations using alternating direction method

A Duality-Preserving Adjoint Method for Segregated Navier–Stokes Solvers

P-Euler Equations and P-Navier-stokes Equations

The adjoint Rayleigh and Orr-Sommerfeld equations: Green function and eigenmodes

J. KSIAM Vol.6, No.1, 1-15, 2002 RELATIONSHIPS AMONG CHARACTERISTIC FINITE ELEMENT METHODS FOR ADVECTION-DIFFUSION PROBLEMS

On a classification of steady solutions to two-dimensional Euler equations

Analysis for Implicit and Implicit-Explicit ADER and DeC Methods for Ordinary Differential Equations, Advection-Diffusion and Advection-Dispersion Equations

Stability analysis and numerical results for some schemes discretising 2D nonconstant coefficient advection–diffusion equations

Twin Brownian particle method for the study of Oberbeck-Boussinesq fluid flows