Grid Adaptation Using Adjoint-Based Error Minimization

We have recently introduced 1 a new adjoint-based grid adaptation methodology that directly solves an error minimization problem. The methodology has been previously applied to the problems for which the true error is available. In this paper, the methodology is extended to general formulations in which the true error is not employed. Two sets of adjoint variables are introduced. The rst set is used to estimate the error in a functional output of interest (e.g., lift, drag, thrust, etc.). The output error is estimated as an inner product of the residual and the adjoint variables associated with the functional output. The second set is used to compute sensitivity of the error estimate with respect to the mesh coordinates. For the practical error functional considered in this paper, only one set of adjoint equations needs to be solved to nd the both sets of adjoint variables. The methodology of grid adaptation combines grid renemen t with grid redistribution. One advantage of this automated grid adaptation methodology is that the number of adjoint solutions used at each optimization iteration is independent of the number of grid nodes. Another advantage is that grid redistribution provides a rigorous metric-free mechanism to relate the estimated error in the functional output to the size, shape, and orientation of grid elements. Thus this approach is very well suited for anisotropic grid adaptation. In this paper, the adjoint-based grid adaptation strategy is tested on several problems governed by the 2-D steady Poisson and Euler equations.

[1]  J. N. Lyness,et al.  Numerical Differentiation of Analytic Functions , 1967 .

[2]  Christopher L. Rumsey,et al.  Prediction of high lift: review of present CFD capability , 2002 .

[3]  Boris Diskin,et al.  Towards Verification of Unstructured-Grid Solvers , 2008 .

[4]  Nail K. Yamaleev,et al.  Error Minimization via Adjoint-Based Anisotropic Grid Adaptation , 2010 .

[5]  Boris Diskin,et al.  Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Inviscid Fluxes , 2011 .

[6]  D. Darmofal,et al.  Review of Output-Based Error Estimation and Mesh Adaptation in Computational Fluid Dynamics , 2011 .

[7]  John C. Vassberg,et al.  Summary of Data from the First AIAA CFD Drag Prediction Workshop , 2002 .

[8]  Timothy J. Barth,et al.  Numerical aspects of computing high Reynolds number flows on unstructured meshes , 1991 .

[9]  Michael A. Park,et al.  Three-Dimensional Turbulent RANS Adjoint-Based Error Correction , 2003 .

[10]  B. Diskin,et al.  Adjoint-based Methodology for Time-Dependent Optimization , 2008 .

[11]  Michael B. Giles,et al.  Adjoint Recovery of Superconvergent Functionals from PDE Approximations , 2000, SIAM Rev..

[12]  O. C. Zienkiewicz,et al.  Adaptive remeshing for compressible flow computations , 1987 .

[13]  D. Venditti,et al.  Anisotropic grid adaptation for functional outputs: application to two-dimensional viscous flows , 2003 .

[14]  W. K. Anderson,et al.  An implicit upwind algorithm for computing turbulent flows on unstructured grids , 1994 .

[15]  W. K. Anderson,et al.  Recent improvements in aerodynamic design optimization on unstructured meshes , 2001 .

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

[17]  M. Fortin,et al.  Anisotropic mesh adaptation: towards user‐independent, mesh‐independent and solver‐independent CFD. Part I: general principles , 2000 .