A Meshfree Isosurface Computation Method for Boundary Element Methods

A Meshfree Isosurface Computation Method for Boundary Element Methods

André Buchau Wolfgang M. Rucker

Universität Stuttgart, Institut für Theorie der Elektrotechnik, Germany

| |
| | Citation



Isosurfaces are an appropriate approach to visualize scalar fields or the absolute value of vector fields in three dimensions. The nodes of the corresponding isosurface mesh are determined using an efficient and accurate isovalue search method. Then, these nodes are typically connected by triangular elements, which are obtained with the help of an adapted advancing front algorithm. An important prerequisite of an isovalue search method is that volume data of the examined field is available in total space. That means, the field values are precomputed in the nodes of an auxiliary post-processing volume mesh or a novel meshfree method is developed that enables both efficient computations of field values in arbitrary points and fast determination of domains with a defined range of field values. If the first approach is applied, a classical isovalue search method is to use an octree scheme to find relevant volume elements, which are intersected by the isosurface. Finally, the surface elements of the isosurface are constructed based on the intersection points of the isosurface with the volume elements. In that case, the accuracy and the computational costs are mainly influenced by the density of the post-processing volume mesh. In contrast, an innovative coupling of established isovalue search methods, fast boundary element method (BEM) techniques, and advancing front meshing algorithms is here presented to compute isosurfaces with high accuracy only using the original BEM model. This novel meshfree method enables very accurate isovalue search methods along with nearly arbitrarily adjustable resolution of the computed isosurface. Furthermore, refinements of the isosurface are also possible, for instance in dependency of the current viewing position. The main idea to realize this meshfree method is to directly combine an octree-based isovalue search method with the octree-based fast multipole method (FMM). 


applied boundary element methods, fast multipole methods, isosurface computations, mesh- free post-processing, visualization


[1] Wilhelms, J. & van Gelder, A., Octrees for faster isosurface generation. ACM Transactions on Graphics, 11(3), pp. 201–227, 1992. http://dx.doi.org/10.1145/130881.130882

[2] Anderson, J.C., Bennett, J.C. & Joy, K.I., Marching diamonds for unstructured meshes. IEEE Visualization, pp. 423–429, 2005. http://dx.doi.org/10.1109/visual.2005.1532825

[3] Schreiner, J., Scheidegger, C.E. & Silva, C.T., High-quality extraction of isosurfaces from regular and irregular grids. IEEE Transactions on Visualization and Computer Graphics, 12(5), pp. 1205–1212, 2006. http://dx.doi.org/10.1109/TVCG.2006.149

[4] Kazhdan, M., Klein, A., Dalal, K. & Hoppe, H., Unconstrained isosurface extraction on arbitrary octrees. Proceedings of Eurographics Symposium on Geometry Processing, 2007.

[5] Dey, T.K. & Levine, J.A., Delaunay meshing of isosurfaces. The Visual Computer, 24(6), pp. 411–422, 2008. http://dx.doi.org/10.1007/s00371-008-0224-1

[6] Martin, S., Shen, H.W. & McCormick, P., Load-Balanced isosurfacing on multi-GPU clusters. Proceedings of the 10th Eurographics Conference on Parallel Graphics and Visualization, pp. 91–100, 2010.

[7] Greengard, L. & Rokhlin, V., A new version of the fast multipole method for the laplace equation in three dimensions. Acta Numerica, pp. 229–269, 1997. http://dx.doi.org/10.1017/S0962492900002725

[8] Buchau, A., Rieger, W. & Rucker, W.M., Fast field computations with the fast multipole method. COMPEL, 20(2), pp. 547–561, 2001. http://dx.doi.org/10.1108/03321640110383861

[9] Buchau, A. & Rucker, M., Meshfree computation of field lines across multiple domains using fast boundary element methods. IEEE Transactions on Magnetics, 51(3), 2015. http://dx.doi.org/10.1109/TMAG.2014.2359520

[10] Buchau, A. & Rucker, W.M., Feasibility of a meshfree post-processing for boundary element methods. Boundary Elements and Other Mesh Reduction Methods, XXXVIII, pp. 327–338, 2015. http://dx.doi.org/10.2495/BEM380261

[11] Buchau, A., Rieger, W. & Rucker, W.M., BEM computations using the fast multipole method in combination with higher order elements and the galerkin method. IEEE Transactions on Magnetics, 37(5), pp. 3181–3185, 2001. http://dx.doi.org/10.1109/20.952572

[12] Kobbelt, L.P. & Botsch, M., An interactive approach to point cloud triangulation. EUROGRAPHICS, 19(3), 2000.

[13] Hornung, A. & Kobbelt, L., Robust reconstruction of watertight 3D models from non-uniformly sampled point clouds without normal information. Eurographics Symposium on Geometry Processing, 2006.

[14] Löhner, R., Progress in grid generation via the advancing front technique. Engineering with Computers, 12, pp. 186–210, 1996. http://dx.doi.org/10.1007/BF01198734

[15] Lang, U. & Wössner, U., Virtual and augmented reality developments for engineering applications. Proceedings of ECCOMAS, Jyväskylä, 2004.

[16] Open Data Protocol (OData), OASIS, Advancing open standards for the information society, available at: www.oasis-open.org