[논문 리뷰] Efficient entropy stable Gauss collocation methods

Jesse Chan, David C. Del Rey Fernández외 1명|arXiv (Cornell University)|2018. 09. 04.

Advanced Numerical Methods in Computational Mathematics참고 문헌 28인용 수 7

한 줄 요약

이 논문은 다각형 및 육면체 요소에 대해 고차수 불연속 갈레르킨 스킴을 위한 효율적인 엔트로피 안정적 가우스 콜로케이션 방법을 제안한다. 이는 모든 자유도 간의 전반적 결합을 제거하기 위해 분리된 합성곱-통합(SBP) 연산자를 도입함으로써 이루어진다. 이 방법은 가우스 체적점(Gauss quadrature points)을 사용하여 곡선 메쉬에서 더 높은 정확도를 달성하면서도, 유량 차분 형식과 면 기반 요소 간 결합을 통해 엔트로피 안정성을 유지한다. 이로 인해 계산 비용과 통신 오버헤드가 크게 감소한다.

ABSTRACT

The construction of high order entropy stable collocation schemes on quadrilateral and hexahedral elements has relied on the use of Gauss-Legendre-Lobatto collocation points and their equivalence with summation-by-parts (SBP) finite difference operators. In this work, we show how to efficiently generalize the construction of semi-discretely entropy stable schemes on tensor product elements to Gauss points and generalized SBP operators. Numerical experiments suggest that the use of Gauss points significantly improves accuracy on curved meshes.

연구 동기 및 목표

• 기존의 가우스 점에서의 엔트로피 안정적 콜로케이션 스킴이 자유도 간 전반적 결합으로 인해 높은 계산 비용과 통신 오버헤드를 야기하는 문제를 해결한다.
• 가우스 체적점을 사용하여 비정렬 다각형 및 육면체 메쉬에서 고차수 불연속 갈레르킨 방법으로 엔트로피 안정성을 확장한다.
• 전체 텐서 곱 결합을 면 기반 결합으로 대체하여 계산 복잡도를 감소시켜 확장성과 효율성을 향상시킨다.
• 가우스 점과 일반화된 SBP 연산자를 활용하여 곡선 메쉬에서 고차수 정확도와 안정성을 유지한다.
• 이전에 단체 요소에서 사용된 분리된 SBP 프레임워크를 텐서 곱 요소에 일반화하여 가우스 노드를 사용한다.

제안 방법

• 체적 및 표면 노드에 모두 작용하는 분리된 SBP 연산자 $\bm{Q}_N$를 제안하여 안정적이고 정확한 도함수 근사치를 가능하게 한다.
• 변분 형식을 사용하여 도함수 근사치를 $\bm{u} = \bm{W}^{-1} \begin{bmatrix}\bm{I} & \bm{V}_f\end{bmatrix}^T \mathrm{diag}(\bm{f}_N) \bm{Q}_N \bm{g}_N$ 로 표현하며, 이는 비대칭 방식으로 경계 보정을 통합한다.
• 유량 차분 형식을 적용하여 체적점에서의 해 상태 간 이중점 엔트로피 보존 유량을 계산함으로써 이산 엔트로피 안정성을 확보한다.
• 이웃 요소 간에 면 값만을 통신하는 면 기반 요소 간 결합 전략을 도입하여 전반적 노드 간 결합을 피한다.
• 텐서 곱 요소에서 가우스 체적점을 사용하여 $d$ 차원에서 두 점 유량 평가 횟수를 $O(N^{2d})$ 에서 $O(N^{d+1})$ 으로 감소시킨다.
• SBP-SAT 프레임워크 내에서 이론을 유도하여, 연속 엔트로피 부등식의 이산 형태를 통해 엔트로피 안정성을 보장한다.

실험 결과

연구 질문

• RQ1 텐서 곱 요소에서 가우스 점에 대해 전반적 결합 없이도 효율적인 엔트로피 안정적 고차수 콜로케이션 스킴을 구성할 수 있는가?
• RQ2 분리된 SBP 연산자를 다각형 및 육면체 요소의 가우스 노드로 일반화하여 통신 및 계산 비용을 줄일 수 있는가?
• RQ3 곡선 메쉬 문제에서 가우스 체적점과 가우스-레지오드-로바토 점(Gauss-Legendre-Lobatto points)을 사용할 경우 정확도와 안정성에 어떤 영향을 미치는가?
• RQ4 유량 차분 형식을 가우스 콜로케이션 스킴에 적응시켜 엔트로피 안정성을 유지하면서 요소 간 결합을 최소화할 수 있는가?
• RQ5 고차수 엔트로피 안정 스킴에서 면 기반 결합이 전체 노드 기반 결합보다 계산 효율에 얼마나 높은 성능 향상을 가져오는가?

주요 결과

• 제안된 방법은 곡선 메쉬에서 고차수 정확도를 달성하며, 수치 실험 결과 GLL 기반 스킴보다 뚜렷한 향상을 보였다.
• 가우스 점의 사용은 특히 부족한 해상도와 난류 흐름 시뮬레이션에서 GLL 점보다 더 높은 정확도를 제공한다.
• 분리된 SBP 연산자는 요소 간 결합을 오직 면 값으로만 제한하여 기존 SBP 형식에서 나타나는 전반적 결합 패tern을 제거한다.
• $d$ 차원에서 요소당 두 점 유량 평가 횟수를 $O(N^{2d})$ 에서 $O(N^{d+1})$ 으로 감소시켜 확장성에 크게 기여한다.
• 유량 차분 형식과 SBP 성질을 통한 이산 엔트로피 부등식을 만족시킴으로써 엄격한 엔트로피 안정성을 유지한다.
• 수치 결과는 2D 및 3D에서 압축성 올레르 방정식의 매끄럽고, 불연속적이며, 난류 해에 대해 안정성과 고차수 수렴성을 확인한다.