[논문 리뷰] The Augmented Fast Marching Method for Level Set Reinitialization

David Salač|arXiv (Cornell University)|2011. 11. 29.

Advanced Numerical Methods in Computational Mathematics참고 문헌 11인용 수 3

한 줄 요약

이 논문은 고정밀도로 부호 거리 함수와 그 일阶 및 이阶 도함수를 계산하는 새로운 레벨 세트 재초기화 기법인 보완된 빠른 마이닝 방법(Augmented Fast Marching Method, AFMM)을 소개한다. |∇ϕ|² = 1 뿐만 아니라 ∇(|∇ϕ|²) = 0 및 ∇∇(|∇ϕ|²) = 0 도 강제함으로써, 조밀한 격자에서도 매끄러운 곡률 필드를 보장하며, 레벨 세트 및 기울기 필드에 대해 2차 정확도를 달성하고 곡률에 대해 1차 정확도를 확보한다.

ABSTRACT

Including derivative information in the modelling of moving interfaces has been proposed as one method to increase the accuracy of numerical schemes with minimal additional cost. Here a new level set reinitialization technique using the fast marching method is presented. This augmented fast marching method will calculate the signed distance function and up to the second-order derivatives of the signed distance function for arbitrary interfaces. In addition to enforcing the condition $| ablaϕ|^2=1$, where $ϕ$ is the level set function, the method ensures that $ abla(| ablaϕ|)^2=0$ and $ abla abla(| ablaϕ|)^2=0$ are also satisfied. Results indicate that for both two- and three-dimensional interfaces the resulting level set and curvature field are smooth even for coarse grids. Convergence results show that using first-order upwind derivatives and the augmented fast marching method result in a second-order accurate level set and gradient field and a first-order accurate curvature field.

연구 동기 및 목표

• 곡률 및 그 도함수가 핵심적인 움직이는 표면 문제에서 정확하고 매끄러운 레벨 세트 재초기화의 필요성을 해결한다.
• 표준 빠른 마이닝 방법이 재초기화 과정에서 고차 도함수 정보를 유지하지 못하는 한계를 극복한다.
• 레벨 세트 함수와 그 도함수가 에이케오날 방정식 및 기하학적 제약 조건과 일치하도록 매끄럽고 일관되게 유지되도록 보장한다.
• 조금 더 넓은 격자에서도 레벨 세트 및 기울기 필드에 대해 2차 정확도를 확보하면서 곡률에 대해 1차 정확도를 유지하는 것을 가능하게 한다.
• 굽힘 힘이 고차 도함수에 의존하는 복잡한 표면, 예를 들어 비세이드의 수치 시뮬레이션을 지원한다.

제안 방법

• 에이케오날 방정식 |∇ϕ| = 1 을 풀면서 |∇ϕ|² 의 기울기 및 헤시안에 대한 추가 제약 조건을 강제하여 빠른 마이닝 방법(FMM)을 도함수 정보를 포함하도록 확장한다.
• 빠른 마이닝 전파 과정에서 부호 거리 함수의 일阶 및 이阶 도함수를 추적하는 보완 상태 변수를 도입한다.
• 공급원 원리와 안정성을 확보하기 위해 공간 도함수 ϕ, ∇ϕ, ∇∇ϕ 를 상향 유한 차분으로 근사한다.
• 해결 과정의 정확한 전파 순서와 수렴을 보장하기 위해 수용된, 시도 중인, 먼 세 개의 노드 집합을 알고리즘 전반에 걸쳐 유지한다.
• 각 격자 점에서 ∇(|∇ϕ|²) = 0 과 ∇∇(|∇ϕ|²) = 0 을 요구함으로써 기하학적 일관성을 확보하여 기울기 크기의 매끄러움을 강제한다.
• 2차원 및 3차원 영역 모두에 동일한 이산화 체계를 사용하여 에이케오날 방정식 및 도함수 제약 조건을 적용한다.

실험 결과

연구 질문

• RQ1 빠른 마이닝 방법을 확장하여 부호 거리 함수뿐만 아니라 그 일阶 및 이阶 도함수를 고정밀도로 계산할 수 있는가?
• RQ2 ∇(|∇ϕ|²) = 0 과 ∇∇(|∇ϕ|²) = 0 을 강제하면 레벨 세트 시뮬레이션에서 곡률 필드가 더 매끄럽고 수치적 안정성이 향상되는가?
• RQ3 구조적 직각 격자에서 보완된 방법의 레벨 세트 함수, 기울기, 곡률에 대한 수렴 속도는 어떠한가?
• RQ4 이 방법은 조밀한 메esh에서도 레벨 세트 및 기울기 필드에 대해 2차 정확도를 달성하면서 곡률에 대해 1차 정확도를 유지할 수 있는가?
• RQ5 곡률 및 그 도함수가 물리적 힘을 결정하는 복잡한 표면, 예를 들어 비세이드의 시뮬레이션에서 이 방법은 어떻게 성능을 발휘하는가?

주요 결과

• 보완된 빠른 마이닝 방법은 부호 거리 함수와 레벨 세트 함수의 최대 이阶 도함수를 고정밀도로 성공적으로 계산한다.
• 이 방법은 ∇(|∇ϕ|²) = 0 과 ∇∇(|∇ϕ|²) = 0 이 충족됨을 보장하여, 조밀한 격자에서도 매끄러운 곡률 필드를 생성한다.
• 수렴 연구 결과, 일阶 상향 유한 차분을 사용할 경우 레벨 세트 및 기울기 필드에 대해 2차 정확도를 보였다.
• 곡률 필드는 이러한 체계에 대해 이론적 기대와 일치하여 1차 정확도를 보였다.
• 2차원 및 3차원에서의 수치 결과는 이 방법이 매끄럽고 안정적이며 기하학적으로 일관된 레벨 세트 표현을 생성함을 확인한다.
• 이 방법은 곡률에 의존하는 힘에 의해 지배되는 표면, 예를 들어 비세이드 막을 정확하게 시뮬레이션할 수 있도록 고차 도함수 정보를 유지함으로써 가능하게 한다.