期刊
JOURNAL OF COMPUTATIONAL PHYSICS
卷 259, 期 -, 页码 617-635出版社
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2013.11.027
关键词
Advection; Two-way wave equation; Level set; Interface tracking
资金
- Stifterverband fur die Deutsche Wissenchaft
- Centre de Cooperation Universitaire Franco-Bavarois
- Institut National des Sciences de l'Univers (CNRS) via the Observatoire Midi-Pyrenees (project ADVECT)
We present a new approach to perform Eulerian level set interface tracking. It consists in advecting the level set function using a second-order two-way wave equation instead of the standard one-way wave advection equation. The resulting numerical schemes are simple to implement, more stable, and significantly less prone to dissipation errors than popular (e.g., WENO) discretizations of the one-way advection equation of higher order, in particular for long advection times. While both the two-way wave advection and the associated numerical schemes were derived previously, these approaches have never been combined with level set advection. Since the level set function to advect is smooth by construction, this ensures the stability of the solution when using the two-way advection equation discretized using centered finite difference schemes. Our numerical tests show that the two-way wave equation approach yields more accurate results than the standard one-way wave equation for the level set advection, at a considerably smaller computational cost. (C) 2013 Elsevier Inc. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据