Error control and propagation in Adaptive Mesh Refinement applied to elliptic equations on quadtree/octree grids

L. Prouvost,A. Belme,D. Fuster
DOI: https://doi.org/10.23967/admos.2023.030
Abstract:In this work we propose a new adaptive mesh refinement (AMR) method applied on isotropic oc-tree/quadtree meshes. The new AMR approach uses a metric-based linear interpolation error estimation [2] extended to square/cubic elements. The analysis of various examples shows that the minimization of the total numerical error can lead to a suboptimal mesh in terms of pure interpolation error. The grids that minimize the error for different values of N (the number of elements imposed) is related to a fixed ratio between the minimal and mean cell size named the compression ratio. Above a certain value, a clear proportionality between the interpolation and the total error allows us to use the former as a criterion to adapt the grid. However, below a certain critical value of the compression ratio, no correlation between both errors is observed and the interpolation error is no longer representative of the total error contained in the solution. Based on these results, we propose to add a model to estimate the discrete minimum grid size and to impose it as an additional constrain to the error minimization problem. The proposed minimum grid size depends on (i) the structure of the solution, (ii) the number of grid points specified and (iii) a security coefficient defined such that it controls the distance between the optimal pure interpolation error and the targeted performance. By increasing this user defined parameter we show that we effectively restrict the range of the minimization problem to regions where we can safely use the local estimation of the interpolation error to drive the mesh adaptation and reduce the total numerical error. The method is implemented in our in-house open-source solver Basilisk [1, 3] and the performance of our new approach is validated on a Poisson-Helmholtz solver and an incompressible Euler
Engineering,Mathematics
What problem does this paper attempt to address?
The paper aims to address the issues of error control and propagation encountered in solving elliptic equations using the Adaptive Mesh Refinement (AMR) method. Specifically, the researchers propose a new AMR method that is applied to isotropic octree/quadtree grids and extends to square/cube elements using a metric-based linear interpolation error estimation method. Through various instance analyses, it was found that simply minimizing the total numerical error may lead to suboptimal grids, meaning that the best grid cannot be obtained when only considering interpolation error. The study also points out that when the compression ratio (the ratio of the smallest cell size to the average cell size) exceeds a certain threshold, there is a significant proportional relationship between interpolation error and total error. At this point, interpolation error can be used as a criterion for adjusting the grid. However, when the compression ratio falls below a critical value, the relationship between the two errors becomes irrelevant, and interpolation error cannot accurately represent the total error in the solution. Based on these results, the authors suggest adding a model to estimate the discrete minimum grid size and use it as an additional constraint in the error minimization problem. The proposed minimum grid size depends on three factors: (i) the structure of the solution, (ii) the specified number of grid points, and (iii) a defined safety factor to control the distance between the optimal pure interpolation error and the target performance. By adding this user-defined parameter, the scope of the minimization problem can be effectively limited, ensuring operation in regions where local interpolation error estimates can be safely used to drive grid adaptation, thereby reducing the total numerical error. This method has been implemented in the open-source solver Basilisk, developed by the research team, and the performance of the new method has been verified on the Poisson-Helmholtz solver and the incompressible Euler equations solver through academic problems with known solutions. These problems of varying complexity have proven to be effective test cases for evaluating the efficiency of the AMR method.