A stitching method for the generation of unstructured meshes for use with co-volume solution techniques

Abstract The successful implementation of co-volume time-domain solution techniques requires the use of high quality, smooth dual meshes. The generation of Delaunay–Voronoi diagrams with these properties, for two-dimensional domains of complicated geometrical shape, is considered. In the adopted approach, near-boundary regions are discretised according to prescribed criteria and the remainder of the computational domain is discretised using an ideal mesh. The two meshes are stitched together to provide a consistent mesh for the complete domain. After smoothing, the resulting mesh is found to be of a quality which exceeds that of meshes produced by standard, automatic, unstructured mesh generation methods. Examples, involving electromagnetic wave scattering, are included to demonstrate the computational performance that can be achieved with a co-volume time domain solution algorithm on these meshes.

[1]  Qiang Du,et al.  Grid generation and optimization based on centroidal Voronoi tessellations , 2002, Appl. Math. Comput..

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

[3]  William H. Frey,et al.  Mesh relaxation: A new technique for improving triangulations , 1991 .

[4]  Alexander G. Churbanov,et al.  A Unified Algorithm to Predict Both Compressible and Incompressible Flows , 2002, Numerical Methods and Application.

[5]  N. Madsen Divergence preserving discrete surface integral methods for Maxwell's curl equations using non-orthogonal unstructured grids , 1995 .

[6]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[7]  S.,et al.  Numerical Solution of Initial Boundary Value Problems Involving Maxwell’s Equations in Isotropic Media , 1966 .

[8]  Barry Hilary Valentine Topping,et al.  A new point creation scheme for uniform Delaunay triangulation , 2002 .

[9]  Stephen D. Gedney,et al.  Full wave analysis of printed microstrip devices using a generalized Yee-algorithm , 1993, Proceedings of IEEE Antennas and Propagation Society International Symposium.

[10]  Roy A. Nicolaides,et al.  Convergence analysis of a covolume scheme for Maxwell's equations in three dimensions , 1998, Math. Comput..

[11]  Nigel P. Weatherill,et al.  Arbitrary order edge elements for electromagnetic scattering simulations using hybrid meshes and a PML , 2002 .

[12]  Jaime Peraire,et al.  A time domain unstructured grid approach to the simulation of electromagnetic scattering in piecewise homogeneous media , 1996 .

[13]  Alexander Fuchs Automatic Grid Generation with Almost Regular Delaunay Tetrahedra , 1998, IMR.

[14]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .