Dynamic tubular grid: An efficient data structure and algorithms for high resolution level sets
2006 (English)In: Journal of Scientific Computing, ISSN 0885-7474, E-ISSN 1573-7691, Vol. 26, no 3, 261-299 p.Article in journal (Refereed) Published
Level set methods [Osher and Sethian. Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations. J. Comput. Phys. 79 (1988) 12] have proved very successful for interface tracking in many different areas of computational science. However, current level set methods are limited by a poor balance between computational efficiency and storage requirements. Tree-based methods have relatively slow access times, whereas narrow band schemes lead to very large memory footprints for high resolution interfaces. In this paper we present a level set scheme for which both computational complexity and storage requirements scale with the size of the interface. Our novel level set data structure and algorithms are fast, cache efficient and allow for a very low memory footprint when representing high resolution level sets. We use a time-dependent and interface adapting grid dubbed the "Dynamic Tubular Grid" or DT-Grid. Additionally, it has been optimized for advanced finite difference schemes currently employed in accurate level set computations. As a key feature of the DT-Grid, the associated interface propagations are not limited to any computational box and can expand freely. We present several numerical evaluations, including a level set simulation on a grid with an effective resolution of 10243. © 2006 Springer Science+Business Media, Inc.
Place, publisher, year, edition, pages
2006. Vol. 26, no 3, 261-299 p.
Compact model representations, Deformable surfaces, Interface tracking, Level sets, Partial differential equations
Engineering and Technology
IdentifiersURN: urn:nbn:se:liu:diva-50280DOI: 10.1007/s10915-005-9062-8OAI: oai:DiVA.org:liu-50280DiVA: diva2:271176