An adjoint method for a high-order discretization of deforming domain conservation laws for optimization of flow problems

JOURNAL OF COMPUTATIONAL PHYSICS

卷:326

Article

Zahr, M. J.¹  Persson, P. -O.² 

[1] Stanford Univ, Inst Computat & Math Engn, Stanford, CA 94035 USA

[2] Univ Calif Berkeley, Dept Math, Berkeley, CA 94720 USA

关键词: Fully discrete adjoint method;    Discontinuous Galerkin methods;    Deforming domain conservation law;    Optimal control;    High-order methods;    PDE-constrained optimization;   

DOI  :  10.1016/j.jcp.2016.09.012

来源: Elsevier

【 摘 要 】

The fully discrete adjoint equations and the corresponding adjoint method are derived for a globally high-order accurate discretization of conservation laws on parametrized, deforming domains. The conservation law on the deforming domain is transformed into one on a fixed reference domain by the introduction of a time-dependent mapping that encapsulates the domain deformation and parametrization, resulting in an Arbitrary Lagrangian-Eulerian form of the governing equations. A high-order discontinuous Galerkin method is used to discretize the transformed equation in space and a high-order diagonally implicit Runge-Kutta scheme is used for the temporal discretization. Quantities of interest that take the form of space-time integrals are discretized in a solver-consistent manner. The corresponding fully discrete adjoint method is used to compute exact gradients of quantities of interest along the manifold of solutions of the fully discrete conservation law. These quantities of interest and their gradients are used in the context of gradient-based PDE-constrained optimization. The adjoint method is used to solve two optimal shape and control problems governed by the isentropic, compressible Navier-Stokes equations. The first optimization problem seeks the energetically optimal trajectory of a 2D airfoil given a required initial and final spatial position. The optimization solver, driven by gradients computed via the adjoint method, reduced the total energy required to complete the specified mission nearly an order of magnitude. The second optimization problem seeks the energetically optimal flapping motion and time-morphed geometry of a 2D airfoil given an equality constraint on the x-directed impulse generated on the airfoil. The optimization solver satisfied the impulse constraint to greater than 8 digits of accuracy and reduced the required energy between a factor of 2 and 10, depending on the value of the impulse constraint, as compared to the nominal configuration. (C) 2016 Elsevier Inc. All rights reserved.

【 授权许可】

Free