A Cost/Benefit Analysis of Simplicial Mesh Improvement Techniques as Measured by Solution Efficiency

The quality of unstructured meshes has long been known to affect both the efficiency and the accuracy of the numerical solution of application problems. Mesh quality can often be improved through the use of algorithms based on local reconnection schemes, node smoothing, and adaptive refinement or coarsening. These methods typically incur a significant cost, and in this paper, we provide an analysis of the tradeoffs associated with the cost of mesh improvement in terms of solution efficiency. We first consider simple finite element applications and show the effect of increasing the number of poor quality elements in the mesh and decreasing their quality on the solution time of a number of different solvers. These simple application problems are theoretically well-understood, and we show the relationship between the quality of the mesh and the eigenvalue spectrum of the resulting linear system. We then consider realistic finite element and finite volume application problems, and show that the cost of mesh improvement is significantly less than the cost of solving the problem on a poorer quality mesh.

[1]  E. Amezua,et al.  A method for the improvement of 3D solid finite-element meshes , 1995 .

[2]  V. Parthasarathy,et al.  A constrained optimization approach to finite element mesh smoothing , 1991 .

[3]  R. K. Smith,et al.  Mesh Smoothing Using A Posteriori Error Estimates , 1997 .

[4]  O. Axelsson,et al.  On the rate of convergence of the preconditioned conjugate gradient method , 1986 .

[5]  C. Ollivier-Gooch Multigrid acceleration of an upwind Euler solver on unstructured meshes , 1995 .

[6]  Lloyd N. Trefethen,et al.  How Fast are Nonsymmetric Matrix Iterations? , 1992, SIAM J. Matrix Anal. Appl..

[7]  A. Jennings Influence of the Eigenvalue Spectrum on the Convergence Rate of the Conjugate Gradient Method , 1977 .

[8]  Charles L. Lawson,et al.  Properties of n-dimensional triangulations , 1986, Comput. Aided Geom. Des..

[9]  S. Canann,et al.  Optismoothing: an optimization-driven approach to mesh smoothing , 1993 .

[10]  B. Joe Three-dimensional triangulations from local transformations , 1989 .

[11]  Barry Joe,et al.  Construction of Three-Dimensional Improved-Quality Triangulations Using Local Transformations , 1995, SIAM J. Sci. Comput..

[12]  I. Fried Condition of finite element matrices generated from nonuniform meshes. , 1972 .

[13]  O. Axelsson A class of iterative methods for finite element equations , 1976 .

[14]  I. Babuska,et al.  ON THE ANGLE CONDITION IN THE FINITE ELEMENT METHOD , 1976 .

[15]  M. Rivara Mesh Refinement Processes Based on the Generalized Bisection of Simplices , 1984 .