[논문 리뷰] 3D Volume Calculation For the Marching Cubes Algorithm in Cartesian Coordinates

Shuqiang Wang|arXiv (Cornell University)|2013. 08. 02.

Computational Geometry and Mesh Generation참고 문헌 5인용 수 6

한 줄 요약

이 논문은 메르팅 큐브 알고리즘에서 전통적인 테트라하드랄 분할 대신 다이버전스 정리(Divergence Theorem)를 사용하여 3차원 카르테시안 격자에서 국소 부피를 계산하는 새로운 방법을 제안한다. 일정한 발산을 가지는 벡터장(예: f = (x,y,z))를 선택함으로써, 등면적 표면에서의 적분을 통해 부피를 계산하게 되어, 테트라하드랄 분할에 비해 더 견고하고 효율적이며 코드 구현에 유리한 접근 방식이 된다.

ABSTRACT

From a scalar field defined at the corner of a cube, an isosurface can be extracted using the Marching Cube algorithm. The isosurface separates the cell into two or more partial cells. A similar situation arises when an material interface in the Front Tracking method cuts through the computational cells. A popular method to calculate the volumes of the partial cells is to first partition the cells into tetrahedra and then sum together the volumes of the tetrahedra for the corresponding partial cells. In this paper, the divergence theorem is used to calculate the volumes of the partial cells generated by the Marching Cubes algorithm. This method is both more robust and efficient compared with the tetrahedralization approach.

연구 동기 및 목표

• 메르팅 큐브 알고리즘에서 부분 세포 부피를 계산하는 데 있어 효율적이고 견고한 방법의 부족을 해결하기 위해.
• 복잡한 인터페이스 구성을 가진 경우 테트라하드랄 분할 기반의 부피 계산의 복잡성과 코드 작성 난이도를 극복하기 위해.
• 부피 적분을 위한 다이버전스 정리를 활용하여 수치적으로 안정적이고 계산 효율성이 높은 대안을 제공하기 위해.
• 다상 및 인터페이스 문제를 위한 임베디드 바운더리 및 프론트 트래킹 방법에서 정확한 부피 계산을 가능하게 하기 위해.
• 모든 23개의 고유한 큐브 구성에 대해 일관되고 라벨링된 템플릿 시스템을 사용하여 구현을 단순화하기 위해.

제안 방법

• 부피는 표면 적분으로 다이버전스 정리를 적용하여 계산한다: 부피 = (1/3) ∫∂Ω n·f dS, 여기서 f = (x,y,z)이면 ∇·f = 3이 된다.
• 벡터장 f = (x,y,z)를 사용함으로써 ∇·f = 3이 되어, 닫힌 인터페이스를 따라 표면 적분을 통해 정확한 부피 계산이 가능해진다.
• 인터페이스는 평면 다각형 면으로 삼각분할되며, 정확한 결과를 얻기 위해 이차 적분 기법을 사용하여 표면 적분을 평가한다.
• 코드 구현을 위해 정점, 모서리, 면을 체계적으로 식별하기 위해 두 글자 코드(예: v0, e0, W, I)를 사용하는 라벨링 체계를 도입한다.
• 대칭성을 활용하여 256개의 가능한 큐브 구성이 23개의 고유한 케이스로 압축되며, 각 케이스에 대해 사전 계산된 삼각분할과 표면 적분 기여도가 제공된다.
• 각 케이스에 대해 한 구성 요소의 부피는 표면 적분을 통해 계산되고, 나머지 부분은 전체 큐브 부피에서 빼내어 계산된다.

실험 결과

연구 질문

• RQ1 다이버전스 정리는 메르팅 큐브 알고리즘에서 3차원 카르테시안 격자에서 국소 부피를 계산하는 데 효과적으로 적용될 수 있는가?
• RQ2 다이버전스 정리 기반 방법은 테트라하드랄 분할에 비해 부피 계산에서 견고성과 효율성 측면에서 어떻게 비교되는가?
• RQ3 모든 23개의 고유한 큐브 구성에 대해 일관되고 코드 구현에 유리한 라벨링 및 삼각분할 체계를 개발할 수 있는가?
• RQ4 이 방법은 정확한 인터페이스 부피 계산이 요구되는 임베디드 바운더리 및 프론트 트래킹 방법에 적합한가?
• RQ5 이 맥락에서 표면 적분을 통한 정확한 부피 계산을 보장하기 위해 어떤 벡터장 선택이 필요한가?

주요 결과

• 다이버전스 정리 방법은 삼각형 메esh에서 이차 적분 기법을 사용하여 정확한 부피 계산을 가능하게 하며, 테트라하드랄 분할에서 발생하는 수치 오차를 피한다.
• 이 방법은 특히 다수의 모서리 및 면 교차를 포함한 복잡한 인터페이스 구조를 처리할 때 테트라하드랄 분할보다 더 견고하고 효율적이다.
• 정점, 모서리, 면/인터페이스 삼각형에 대해 일관된 라벨링 체계(예: v0, e0, I, W)를 사용함으로써 코드 작성 과정이 단순화된다.
• 23개의 고유한 큐브 구성은 삼각분할과 표면 적분 기여도를 사전 계산하여 완전히 열거되어 있다.
• 이 방법은 타원형/포물형 인터페이스 문제, 이상성 불가약성 유체 흐름, 자기유체역학적 흐름 등의 응용 분야에서 정확한 부피 계산을 가능하게 한다.
• 큐브의 총 부피를 사용하여 보완되는 구성 요소의 부피를 계산함으로써 계산의 중복을 줄인다.