An adjoint view on flux consistency and strong wall boundary conditions to the Navier-Stokes equations

Inconsistent discrete expressions in the boundary treatment of Navier-Stokes solvers and in the definition of force objective functionals can lead to discrete-adjoint boundary treatments that are not a valid representation of the boundary conditions to the corresponding adjoint partial differential equations. The underlying problem is studied for an elementary 1D advection-diffusion problem first using a node-centred finite-volume discretisation. The defect of the boundary operators in the inconsistently defined discrete-adjoint problem leads to oscillations and becomes evident with the additional insight of the continuous-adjoint approach. A homogenisation of the discretisations for the primal boundary treatment and the force objective functional yields second-order functional accuracy and eliminates the defect in the discrete-adjoint boundary treatment. Subsequently, the issue is studied for aerodynamic Reynolds-averaged Navier-Stokes problems in conjunction with a standard finite-volume discretisation on median-dual grids and a strong implementation of noslip walls, found in many unstructured general-purpose flow solvers. Going out from a base-line discretisation of force objective functionals which is independent of the boundary treatment in the flow solver, two improved flux-consistent schemes are presented; based on either body wall-defined or farfield-defined control-volumes they resolve the dual inconsistency. The behaviour of the schemes is investigated on a sequence of grids in 2D and 3D.

[1]  Ralf Hartmann,et al.  Adjoint-based error estimation and adaptive mesh refinement for the RANS and k-ω turbulence model equations , 2011, J. Comput. Phys..

[2]  Jason E. Hicken,et al.  Dual consistency and functional accuracy: a finite-difference perspective , 2014, J. Comput. Phys..

[3]  Ralf Hartmann,et al.  Adjoint Consistency Analysis of Discontinuous Galerkin Discretizations , 2007, SIAM J. Numer. Anal..

[4]  M. Giles,et al.  Adjoint Code Developments Using the Exact Discrete Approach , 2001 .

[5]  Kyriakos C. Giannakoglou,et al.  Adjoint wall functions: A new concept for use in aerodynamic shape optimization , 2010, J. Comput. Phys..

[6]  D. Mavriplis.  Discrete Adjoint-Based Approach for Optimization Problems on Three-Dimensional Unstructured Meshes , 2007 .

[7]  R. Dwight,et al.  Efficient and robust algorithms for solution of the adjoint compressible Navier–Stokes equations with applications , 2009 .

[8]  Norbert Kroll,et al.  The DLR Flow Solver TAU - Status and Recent Algorithmic Developments , 2014 .

[9]  James Lu,et al.  An a posteriori Error Control Framework for Adaptive Precision Optimization using Discontinuous Galerkin Finite Element Method , 2005 .

[10]  M. Giles,et al.  Algorithm Developments for Discrete Adjoint Methods , 2003 .

[11]  Rainald Löhner,et al.  An adjoint‐based design methodology for CFD problems , 2004 .

[12]  David L. Darmofal,et al.  Analysis of Dual Consistency for Discontinuous Galerkin Discretizations of Source Terms , 2009, SIAM J. Numer. Anal..

[13]  Antony Jameson,et al.  Aerodynamic design via control theory , 1988, J. Sci. Comput..

[14]  Jason E. Hicken,et al.  Superconvergent Functional Estimates from Summation-By-Parts Finite-Difference Discretizations , 2011, SIAM J. Sci. Comput..

[15]  Michael B. Giles,et al.  On the use of Runge-Kutta time-marching and multigrid for the solution of steady adjoint equations , 2000 .

[16]  Thomas Kaminski,et al.  Recipes for adjoint code construction , 1998, TOMS.

[17]  Arthur Stück,et al.  Adjoint complement to viscous finite-volume pressure-correction methods , 2013, J. Comput. Phys..

[18]  Forrester T. Johnson,et al.  Modi cations and Clari cations for the Implementation of the Spalart-Allmaras Turbulence Model , 2011 .

[19]  Adrian Sandu,et al.  On the properties of discrete adjoints of numerical methods for the advection equation , 2008 .

[20]  Ralf Hartmann,et al.  Higher order and adaptive DG methods for compressible flows , 2013 .

[21]  R. Swanson,et al.  Multistage Schemes With Multigrid for Euler and Navier-Stokes Equations , 1997 .

[22]  A. Jameson.  Optimum aerodynamic design using CFD and control theory , 1995 .

[23]  Thomas H. Pulliam,et al.  Artificial Dissipation Models for the Euler Equations , 1985 .

[24]  Jan Nordström,et al.  Weak and strong wall boundary procedures and convergence to steady-state of the Navier-Stokes equations , 2012, J. Comput. Phys..

[25]  Ralf Heinrich,et al.  The DLR TAU-Code: Recent Applications in Research and Industry , 2006 .

[26]  D. Darmofal,et al.  An implicit, exact dual adjoint solution method for turbulent flows on unstructured grids , 2004 .

[27]  Jason E. Hicken,et al.  The Role of Dual Consistency in Functional Accuracy: Error Estimation and Superconvergence , 2011 .