[논문 리뷰] Space-time Galerkin isogeometric method and efficient solver for parabolic problems

Gabriele Loli, Monica Montardini외 2명|arXiv (Cornell University)|2019. 09. 16.

Advanced Numerical Analysis Techniques인용 수 4

한 줄 요약

이 논문은 텐서 곱 구조를 기반으로 한 새로운 조건자(preconditioner)와 확장된 빠른 대각화 방법을 활용한 공간-시간 갈레르킨 이소지오메트릭 방법을 제안한다. 이는 편미분방정식의 효율적 해법을 가능하게 하며, 자유도 수에 대해 거의 선형적 스케일러빌리티를 보장하고, 다항식 차수와 복잡한 기하 구조, 변수 계수를 포함한 상황에서도 안정성을 유지한다.

ABSTRACT

In this work we focus on the preconditioning of a Galerkin space-time isogeometric discretization of the heat equation. Exploiting the tensor product structure of the basis functions in the parametric domain, we propose a preconditioner that is the sum of Kronecker products of matrices and that can be efficiently applied thanks to an extension of the classical Fast Diagonalization method. The preconditioner is robust w.r.t. polynomial degree and the time required for the application is almost proportional to the number of degrees-of-freedom, for a serial execution. By incorporating some information on the geometry parametrization and on the equation coefficients, we keep high efficiency with non-trivial domains and variable thermal conductivity and heat capacity coefficients.

연구 동기 및 목표

• 공간-시간 이소지오메트릭 이산화를 통한 편미분방정식 문제의 계산 비효율성 문제를 해결한다.
• 다양한 다항식 차수와 복잡한 기하 구조에서 효율성을 유지하는 조건자를 개발한다.
• 기하 구조의 매개변수화와 변수 계수(예: 열전도도)를 해법에 통합하면서 성능 저하 없이 유지한다.
• 직렬 실행 환경에서 해법 시간 복잡도를 거의 선형으로 확보한다.

제안 방법

• 텐서 곱 기반의 기저 함수를 사용한 열 방정식에 대한 갈레르킨 공간-시간 이소지오메트릭 이산화를 활용한다.
• 공간 및 시간 행렬의 텐서 곱의 합으로 구성된 조건자를 통해 텐서 구조를 활용한다.
• 기존의 빠른 대각화 방법을 확장하여 조건자를 거의 최적의 시간 내에 효율적으로 적용한다.
• 기하 구조의 매개변수화와 변수 계수(예: 열전도도 및 비열)를 직접 조건자 설계에 통합한다.
• 비트ivial한 도메인과 비정수 계수 조건에서도 조건자가 효과적이고 효율적으로 유지됨을 보장한다.

실험 결과

연구 질문

• RQ1 공간-시간 이소지오메트릭 방법이 텐서 곱 기반 조건자를 사용하여 편미분방정식 문제를 거의 선형 시간 복잡도로 해결할 수 있는가?
• RQ2 제안된 조건자가 이소지오메트릭 이산화에서 다양한 다항식 차수에 대해 어떻게 성능을 유지하는가?
• RQ3 복잡한 기하 구조와 변수 계수를 처리할 때 해법의 효율성은 어느 정도 유지되는가?
• RQ4 확장된 빠른 대각화 방법이 직렬 실행 환경에서 조건자를 효율적으로 적용할 수 있는가?

주요 결과

• 제안된 조건자는 직렬 실행에서 자유도 수에 대해 거의 선형적으로 해법 시간이 증가함을 보여준다.
• 이소지오메트릭 기저 함수의 다항식 차수에 대해 해법이 안정적으로 유지된다.
• 조건자에 기하 구조의 매개변수화와 변수 계수를 통합하더라도 해법 효율성이 저하되지 않는다.
• 기저 함수의 텐서 곱 구조를 활용하여 확장된 빠른 대각화 방법이 조건자를 효과적으로 적용한다.