A method is introduced to decrease the computational labor of the standard level set method for propagating interfaces. The fast approach uses only points close to the curve at every time step. We describe this new algorithm and compare its efficiency and accuracy with the standard level set approac
A PDE-Based Fast Local Level Set Method
β Scribed by Danping Peng; Barry Merriman; Stanley Osher; Hongkai Zhao; Myungjoo Kang
- Publisher
- Elsevier Science
- Year
- 1999
- Tongue
- English
- Weight
- 431 KB
- Volume
- 155
- Category
- Article
- ISSN
- 0021-9991
No coin nor oath required. For personal study only.
β¦ Synopsis
We develop a fast method to localize the level set method of Osher and Sethian (1988, J. Comput. Phys. 79, 12) and address two important issues that are intrinsic to the level set method: (a) how to extend a quantity that is given only on the interface to a neighborhood of the interface; (b) how to reset the level set function to be a signed distance function to the interface efficiently without appreciably moving the interface. This fast local level set method reduces the computational effort by one order of magnitude, works in as much generality as the original one, and is conceptually simple and easy to implement. Our approach differs from previous related works in that we extract all the information needed from the level set function (or functions in multiphase flow) and do not need to find explicitly the location of the interface in the space domain. The complexity of our method to do tasks such as extension and distance reinitialization is O(N ), where N is the number of points in space, not O(N log N ) as in works by Sethian (1996, Proc. Nat. Acad. Sci. 93, 1591) and Helmsen and co-workers (1996, SPIE Microlithography IX, p. 253). This complexity estimation is also valid for quite general geometrically based front motion for our localized method.
π SIMILAR VOLUMES
Level set techniques are numerical techniques for tracking the evolution of interfaces. They rely on two central embeddings; first, the embedding of the interface as the zero level set of a higher dimensional function, and second, the embedding (or extension) of the interface's velocity to this high
Modeling of microstructural evolution during thin-film deposition requires a knowledge of several key activation energies (surface diffusion, island edge atom diffusion, adatom migration over descending step edges, etc.). These and other parameters must be known as a function of crystal orientation.
We present a level set based numerical algorithm for simulating a model of epitaxial growth. The island dynamics model is a continuum model for the growth of thin films. In this paper, we emphasize the details of the numerical method used to simulate the island dynamics model.