Implementation and Assessment of a Residual-Based r-Adaptation Technique on Structured Meshes

Abstract

In this study, an r-adaptation technique for mesh adaptation is employed for reducing the solution discretization error, which is the error introduced due to spatial and temporal discretization of the continuous governing equations in numerical simulations. In r-adaptation, mesh modification is achieved by relocating the mesh nodes from one region to another without introducing additional nodes. Truncation error (TE) or the discrete residual is the difference between the continuous and discrete form of the governing equations. Based upon the knowledge that the discrete residual acts as the source of the discretization error in the domain, this study uses discrete residual as the adaptation driver. The r-adaptation technique employed here uses structured meshes and is verified using a series of one-dimensional (1D) and two-dimensional (2D) benchmark problems for which exact solutions are readily available. These benchmark problems include 1D Burgers equation, quasi-1D nozzle flow, 2D compression/expansion turns, and 2D incompressible flow past a Karman–Trefftz airfoil. The effectiveness of the proposed technique is evident for these problems where approximately an order of magnitude reduction in discretization error (when compared with uniform mesh results) is achieved. For all problems, mesh modification is compared using different schemes from literature including an adaptive Poisson grid generator (APGG), a variational grid generator (VGG), a scheme based on a center of mass (COM) analogy, and a scheme based on deforming maps. In addition, several challenges in applying the proposed technique to real-world problems are outlined.