[논문 리뷰] Spacetime Meshing for Discontinuous Galerkin Methods

Thite, Shripad|ArXiv.org|2008. 04. 07.

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

한 줄 요약

이 논문은 비연속 갈레르킨 방법을 위한 새로운 전진 면 시공간 메쉬 알고리즘을 제안하며, (d+1)-차원 시공간 영역에서 비정규, 적응형 단체 메쉬를 구축한다. 원뿔 제약 조건과 두 개의 파라미터 최적화를 활용하여 시공간 사면체의 높이를 최대화함으로써, 유한한 시간적 비율을 확보하고, 동적 메쉬 세분화 및 굴착을 지원하며, 비선형성과 이방성 행동을 보이는 PDE의 병렬 해법을 가능하게 한다.

ABSTRACT

Spacetime discontinuous Galerkin (SDG) finite element methods are used to solve such PDEs involving space and time variables arising from wave propagation phenomena in important applications in science and engineering. To support an accurate and efficient solution procedure using SDG methods and to exploit the flexibility of these methods, we give a meshing algorithm to construct an unstructured simplicial spacetime mesh over an arbitrary simplicial space domain. Our algorithm is the first spacetime meshing algorithm suitable for efficient solution of nonlinear phenomena in anisotropic media using novel discontinuous Galerkin finite element methods for implicit solutions directly in spacetime. Given a triangulated d-dimensional Euclidean space domain M (a simplicial complex) and initial conditions of the underlying hyperbolic spacetime PDE, we construct an unstructured simplicial mesh of the (d+1)-dimensional spacetime domain M x [0,infinity). Our algorithm uses a near-optimal number of spacetime elements, each with bounded temporal aspect ratio for any finite prefix M x [0,T] of spacetime. Our algorithm is an advancing front procedure that constructs the spacetime mesh incrementally, an extension of the Tent Pitcher algorithm of Ungor and Sheffer (2000). In 2DxTime, our algorithm simultaneously adapts the size and shape of spacetime tetrahedra to a spacetime error indicator. We are able to incorporate more general front modification operations, such as edge flips and limited mesh smoothing. Our algorithm represents recent progress towards a meshing algorithm in 2DxTime to track moving domain boundaries and other singular surfaces such as shock fronts.

연구 동기 및 목표

• 불연속 갈레르킨(SDG) 방법을 위한 메쉬 알고리즘 개발로, 불연속성이 있는 초구형 PDE의 정확하고 효율적인 해법을 지원한다.
• 비선형성과 이방성 매질에 적응하는 비정규 시공간 메쉬를 구축하면서 기하학적 및 수치적 안정성을 유지하는 도전 과제를 해결한다.
• 새로 추가된 시공간 요소에서 즉각적인 계산을 가능하게 하여, PDE의 증분적, 병렬적 해법을 가능하게 한다.
• 동적 세분화, 굴착, 모서리 뒤집기, 제한된 스무딩을 포함한 동적 메쉬 적응성과 면 전진을 통합한다.
• 실용적 개선을 위해 적응형 파라미터 선택을 지원하면서도, 시간적 비율과 메쉬 품질에 대한 이론적 보장을 제공한다.

제안 방법

• 시공간 메쉬 세분화를 두 개의 파라미터 최적화 문제로 공식화하여, 시공간 사면체의 높이를 극대화함으로써 시간적 및 공간적 전진를 균형 잡는다.
• 메쉬 면의 기울기 조건을 강제하기 위해 원뿔 제약 조건을 사용하여 인과성 구조를 보장하고, 견고한 면 전파를 가능하게 한다.
• 텐트 피치어 알고리즘에 영감을 얻은 전진 면 절차를 구현하며, 모서리 뒤집기 및 제한된 메쉬 스무딩을 포함한 동적 면 수정 연산을 확장한다.
• 비탐욕적 면 전진 전략을 도입하여 인접한 면을 굴착하고, 공면성(공면성)을 촉진함으로써 향후 메쉬 품질을 향상시킨다.
• 모서리 연장 및 수축과 같은 국소 메쉬 최적화 기법을 적용하여 고품질 시공간 요소의 수를 증가시킨다.
• 모든 요소가 약한 단체 복합체를 이루고, 최소한 하나의 인과적 외부 유출 면을 가짐을 보장함으로써 메쉬 유효성을 유지한다. (조합적 변화 조건 하에서도 동일하게 유지됨)

실험 결과

연구 질문

• RQ1 비선형성과 이방성 PDE 반응에 적응하는 시공간 메쉬 세분화는 어떻게 비선형성과 이방성 PDE 반응에 적응하면서도 유한한 시간적 비율을 유지할 수 있는가?
• RQ2 비연속 갈레르킨 방법을 사용할 때, 초구형 PDE의 효율적이고 증분적, 병렬 처리 가능한 해법을 위한 알고리즘적 메커니즘은 무엇인가?
• RQ3 미래 기하학적 변화를 예측하는 비탐욕적 면 전진 전략을 통해 메쉬 품질을 향상시킬 수 있는가?
• RQ4 SDG 방법에 적합한 메쉬에서 시공간 요소의 수와 그 비율에 대한 이론적 한계는 무엇인가?
• RQ5 모서리 뒤집기 및 굴착과 같은 면 연산은 동적 도메인 경계와 충격 전면을 지원하기 위해 전진 면 프레임워크에 어떻게 통합될 수 있는가?

주요 결과

• 알고리즘은 유한한 시간 프리픽스 T에 대해 d차원 공간 도메인 위에 비정규 단체 시공간 메쉬를 구성하며, 이는 임의의 유한한 시간에 대해 거의 최적의 요소 수를 갖는다.
• 모든 시공간 요소에서 유한한 시간적 비율을 유지하여, SDG 방법에서 수치적 안정성을 보장한다.
• 모서리 뒤집기 및 제한된 메쉬 스무딩을 포함한 동적 세분화 및 굴착을 지원하여 실질적인 요소 품질 향상을 가능하게 한다.
• 시공간 요소의 수는 동일한 조건 하에서 크기 최적화 메쉬의 수에 대해 증명 가능하게 가까운 수준이다.
• 제한된 메쉬 스무딩이 실용적 품질 향상에 기여함에도 불구하고, 고차원 d ≥ 3에서도 일정 요소의 상수 배수 이내로 진전을 보장한다.
• 실험 결과는 면의 공간 투영 품질 향상이 더 우수한 시공간 요소 품질로 이어진다는 것을 확인한다.