Nondegeneracy tests for hexahedral cells

Abstract The aims of the paper are to consider the nondegeneracy requirement for computational grids and to analyze eight tests used to check the nondegeneracy of hexahedral cells. The paper starts with consideration of nondegeneracy requirement and formulation of definitions and common theorems utilized for estimation of nondegeneracy of grids both structured and unstructured. Then hexahedral cells are introduced and sufficient nondegeneracy conditions (Ushakova, 2000) for them are given. Sufficient nondegeneracy conditions are 27 inequalities for 32 tetrahedral volumes. Besides sufficient nondegeneracy conditions other conditions are applied as nondegeneracy tests in grid generation theory and practice. Considered nondegeneracy tests are the checks for positivity of different values. Tests 1, 2, 3, 4, 5, 6 check the positivity of 8, 10, 24, 32, 58, 48 tetrahedral volumes, correspondingly. Test 7 verifies the positivity of the volume of a cell. Test 8 checks the positivity of the Jacobian of the mapping used for generation of a cell. The check is performed at the corners of a cell and hex center. Tests 1, 7, 8 are often used in commercial packages. For the most part, nondegeneracy tests are not sufficient nondegeneracy conditions, however they are used for the purpose of constructing nondegenerate grids and, some times, instead of sufficient nondegeneracy conditions. The effectiveness and reliability of such substitutions are investigated in special numerical experiments with random numbers. In the numerical experiment for each test, hexahedral cells are generated randomly. Results of such experiments are the following. Among eight tests, test 2 is considered the best since it verifies the volumes of only 10 tetrahedra for positiveness, guarantees the nondegeneracy in most of cases (68.7% randomly generated hexahedral cells satisfying test 2) and covers a wide class of cells (about 60% of nondegenerate cells). Tests 1, 3, 4, 5, 6, 7, 8 have success in 31.7%, 83.1%, 100%, 100%, 39.5%, 0.2%, 34% of cases and cover 100%, 7.9%, 7.9%, 4.2%, 59.5%, 100%, 100% of nondegenerate cells, correspondingly. Because of high rate of success, tests 3, 4, 5 also can be used for grid generation purpose. All tests are illustrated by the examples of structured grids.

[1]  S. A. Ivanenko,et al.  Piecewise smooth homeomorphisms of bounded domains and their applications to the theory of grids , 2003 .

[2]  Olga V. Ushakova Conditions of Nondegeneracy of Three-Dimensional Cells. A Formula of a Volume of Cells , 2001, SIAM J. Sci. Comput..

[3]  Boris N. Azarenok A variational hexahedral grid generator with control metric , 2006, J. Comput. Phys..

[4]  Stephen A. Vavasis,et al.  A Bernstein-Bezier Sufficient Condition for Invertibility of Polynomial Mapping Functions , 2003, ArXiv.

[5]  S. S. Sritharan 10. Mathematical Aspects of Harmonic Grid Generation , 1991, Mathematical Aspects of Numerical Grid Generation.

[6]  J. Castillo Mathematical Aspects of Numerical Grid Generation , 1991, Frontiers in Applied Mathematics.

[7]  Jeffrey Grandy Conservative Remapping and Region Overlays by Intersecting Arbitrary Polyhedra , 1999 .

[8]  O Ushakova,et al.  Variational Methods of Construction of Optimal Grids , 1998 .

[9]  Bharat K. Soni,et al.  Handbook of Grid Generation , 1998 .

[10]  F. B. Fuller,et al.  HARMONIC MAPPINGS. , 1954, Proceedings of the National Academy of Sciences of the United States of America.

[11]  T. N. Bronina An algorithm for constructing initial three-dimensional structured grids for domains of revolution , 2008 .

[12]  Vladimir D. Liseikin,et al.  A Computational Differential Geometry Approach to Grid Generation (Scientific Computation) , 2003 .

[13]  Granino A. Korn,et al.  Mathematical handbook for scientists and engineers , 1961 .

[14]  M. Prokhorova Problems of homeomorphism arising in the theory of grid generation , 2008 .

[15]  Leon O. Chua,et al.  Global homeomorphism of vector-valued functions , 1972 .

[16]  O. V. Ushakova,et al.  Optimization algorithms for three-dimensional grids in domains of rotation , 2008 .

[17]  P. M. Knupp On the ivertability of the isoparametric map , 1990 .