Authors

Navamita Ray

Type

Text

Type

Dissertation

Advisor

Jiao, Xiangmin | Glimm, James | Li, Xiaolin | Tautges, Timothy.

Date

2013-12-01

Keywords

Applied mathematics | high-order accuracy, mesh data structures, surface integrals, surface reconstruction, surface remeshing

Department

Department of Applied Mathematics and Statistics.

Language

en_US

Source

This work is sponsored by the Stony Brook University Graduate School in compliance with the requirements for completion of degree.

Identifier

http://hdl.handle.net/11401/76305

Publisher

The Graduate School, Stony Brook University: Stony Brook, NY.

Format

application/pdf

Abstract

High-order surface reconstruction is a numerical technique to obtain high-order approximations of both geometry and its differential quantities such as normals, curvatures, etc. | over a discrete surface mesh. Its computational framework is based on local polynomial fittings using a weighted least squares approach. In this dissertation, we complete the scope of this framework to compute high-order approximations of surface integrals and demonstrate the application of the complete framework to various mesh-based numerical computations for high-order numerical methods. The computational framework relies on an efficient underlying mesh data structure for various traversal queries. For this purpose, an array-based mesh data structure was developed to represent the mesh for efficient mesh query and modification operations. Our methods are mainly developed for applications with high-order methods in mind. Surface integration is a fundamental operation in many scientific and engineering applications. The standard methods for numerical computation are generally limited to second-order of accuracy due to lower-order approximations to geometry and integrand. This limitation is overcome by extending the computational framework for high-order surface reconstruction to a function defined over the surface and coupling it with high-order quadrature rules. We theoretically analyze the accuracy of our method and prove that it can achieve high-order of accuracy and verify it with numerical experiments as well. A widely used operation by many applications is the modification of the surface mesh by vertex redistribution, edge flipping, refinement or coarsening, such that the resulting mesh improves certain properties such as mesh quality, error distribution, etc. It is vital to preserve the geometric accuracy of the mesh as it undergoes the modification operations. Our computational framework provides an efficient high-order point projection strategy that can be easily coupled with various mesh quality improving techniques. We develop remeshing strategies coupling existing mesh quality improving techniques with high-order surface reconstruction, to produce high-quality and high-order accurate surface meshes. The developed algorithms are made robust to allow untangling mildly folded triangles and also take into account the approximation issues related to high-order approximations in under-resolved regions. All of our algorithms are based on an array-based half-facet mesh data structure called AHF, for efficient mesh query and modification operations. It was developed for 2D/3D non-manifold meshes with mixed-dimensional submeshes for increased applicability. We present the theoretical framework of our methods, show experimental comparisons against other methods, and demonstrate their utilization to geometric PDE's, high-order finite elements, biomedical image-based surface meshes, and complex interface meshes in fluid simulations. | 111 pages

Share

COinS
 
 

To view the content in your browser, please download Adobe Reader or, alternately,
you may Download the file to your hard drive.

NOTE: The latest versions of Adobe Reader do not support viewing PDF files within Firefox on Mac OS and if you are using a modern (Intel) Mac, there is no official plugin for viewing PDF files within the browser window.