[논문 리뷰] High-Order Quadrature on Multi-Component Domains Implicitly Defined by Multivariate Polynomials

Robert I. Saye|arXiv (Cornell University)|2021. 05. 19.

Advanced Numerical Analysis Techniques참고 문헌 81인용 수 28

한 줄 요약

이 논문은 다변수 다항식으로 암묵적으로 정의된 복잡한 다성분 도메인에서의 적분을 위한 고차수 구적법을 제시한다. 기하 구조를 다가치 높이 함수로 재구성하고, 일차원 가우스-레지앙드르 또는 탄제인-신함 구적법을 적용함으로써, h-재분할에 대해 최대 22차 정확도를 달성하고 q-재분할에 대해 약 지수 수렴을 보이며, 특이점과 위상적 복잡성을 자동으로 처리한다.

ABSTRACT

A high-order quadrature algorithm is presented for computing integrals over curved surfaces and volumes whose geometry is implicitly defined by the level sets of (one or more) multivariate polynomials. The algorithm recasts the implicitly defined geometry as the graph of an implicitly defined, multi-valued height function, and applies a dimension reduction approach needing only one-dimensional quadrature. In particular, we explore the use of Gauss-Legendre and tanh-sinh methods and demonstrate that the quadrature algorithm inherits their high-order convergence rates. Under the action of $h$-refinement with $q$ fixed, the quadrature schemes yield an order of accuracy of $2q$, where $q$ is the one-dimensional node count; numerical experiments demonstrate up to 22nd order. Under the action of $q$-refinement with the geometry fixed, the convergence is approximately exponential, i.e., doubling $q$ approximately doubles the number of accurate digits of the computed integral. Complex geometry is automatically handled by the algorithm, including, e.g., multi-component domains, tunnels, and junctions arising from multiple polynomial level sets, as well as self-intersections, cusps, and other kinds of singularities. A variety of numerical experiments demonstrates the quadrature algorithm on two- and three-dimensional problems, including: randomly generated geometry involving multiple high-curvature pieces; challenging examples involving high degree singularities such as cusps; adaptation to simplex constraint cells in addition to hyperrectangular constraint cells; and boolean operations to compute integrals on overlapping domains.

연구 동기 및 목표

• 다변수 다항식의 영수준집합으로 정의된 곡면 및 부피에서의 적분을 위한 고차수 구적법을 개발하는 것.
• 메esh 또는 수동 분할 없이도 다성분 도메인, 터널, 교차점, 콩기, 자가교차 등 복잡한 기하 구조를 처리하는 것.
• 고차수 수렴을 보이는 1차원 구적법 규칙을 활용하여 정확도와 효율성을 확보하는 것.
• 암묵적으로 정의된 도메인의 개별 연결 성분을 자동으로 탐지하고, 적분 함수가 표면을 가로질러 비연속일 경우에도 정확도를 유지하는 것.
• 최소한의 기하 전처리로 2차원과 3차원 모두에서 동일하게 작동하는 차원 독립적 프레임워크를 제공하는 것.

제안 방법

• 암묵적으로 정의된 도메인을 다가치 높이 함수의 그래프로 재구성함으로써, 적분 문제를 좌표 방향에 따른 1차원 구적법으로 환원한다.
• 높은 정확도를 확보하기 위해 높이 함수에 대해 고차수 1차원 구적법 규칙—특히 가우스-레지앙드르 및 탄제인-신함—을 적용하여 도메인을 통해의 적분을 계산한다.
• 각 적분 셀에 대해, 탈카스텔리오 분할, 뉴턴 방법과 바이섹션 보호 장치, 그리고 안정성을 확보하기 위한 고유값 해법을 조합한 하이브리드 근 구하기 전략을 사용하여 (0,1) 구간 내 다항식 시스템의 모든 실근을 식별한다.
• 구적법 가중치는 정수 함수를 정확히 적분하도록 정규화되어 있으며, 탄제인-신함 규칙의 근사 지수 수렴 특성을 유지한다.
• 이 방법은 차원에 독립적이며, 부피적 및 표면적 적분에 동일하게 적용되며, 각각 d 차원에서 O(q^d) 및 O(q^{d-1})개의 노드를 사용한다.
• 이 방법은 자연스럽게 복잡한 위상을 처리한다: 다중 성분, 교차, 특이점은 다항식 수준집합에서 유기적으로 발생하며 추가 사용자 간섭 없이도 작동한다.

실험 결과

연구 질문

• RQ1 복잡한 위상을 지닌 암묵적으로 정의된 도메인에 대해, 오직 1차원 구적법만을 사용하여 고차수 구적법을 달성할 수 있는가?
• RQ2 이러한 도메인에 대해 h-재분할(메쉬 크기)과 q-재분할(구적법 순서)에 따른 수렴 속도는 어떻게 변화하는가?
• RQ3 메쉬 또는 분할 없이도 다중 연결 성분, 교차점, 콩기와 같은 특이점을 자동으로 탐지하고 적분할 수 있는가?
• RQ4 가우스-레지앙드르 및 탄제인-신함 구적법 규칙의 사용이 곡면 인터페이스와 기하학적 특이점이 존재하는 상황에서도 고차수 정확도를 유지하는 데 얼마나 기여하는가?
• RQ5 하이브리드 근 구하기 전략은 고차수 다항식 시스템에서의 인터페이스 교차점을 어떻게 안정적이고 정확하게 계산하는가?

주요 결과

• 고정된 q에서 h-재분할에 대해 약 2q차 정확도를 달성하며, 수치 실험을 통해 최대 22차 수렴이 확인되었다.
• 고정된 기하 구조에서 q-재분할에 대해 수렴은 약 지수 수렴을 보이며, q를 두 배로 늘일 경우 계산된 적분의 정확한 자릿수 약 두 배로 증가한다.
• 특수 처리나 메쉬 재분할 없이도 고곡률 특징(콩기, 자가교차 포함)이 있는 도메인의 적분도 성공적으로 계산한다.
• 사용자 정의 분해 없이도 다중 다항식 수준집합에서 유래한 다중 연결 성분, 터널, 교차점 등을 자동으로 탐지하고 적분한다.
• 탈카스텔리오 분할, 뉴턴 방법, 고유값 해법을 조합한 하이브리드 근 구하기 전략은 높은 정확도와 안정성을 확보하였으며, 10회 이내의 뉴턴 반복 내에 95% 이상의 근이 정확히 기계 정밀도까지 계산되었다.
• 정규화된 탄제인-신함 구적법을 사용함으로써 정수 함수를 정확히 적분할 수 있으며, 매끄러운 적분 함수에 대해 메서드의 근사 지수 수렴 특성을 유지한다.