[논문 리뷰] Compressive Space-Time Galerkin Discretizations of Parabolic Partial Differential Equations

Stig Larsson, Christoph Schwab|arXiv (Cornell University)|2015. 01. 19.

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

한 줄 요약

이 논문은 시간 도함수의 순서가 1/2인 분수적 미분을 바탕으로 한 변분 형식을 사용하여 선형 포아송 PDE에 대해 압축 가능한 공간-시간 갈레르킨 이산화를 제안한다. 희소 텐서곱 공간에서 이산 인프-서프 조건을 확립하여, 최적 수렴 속도(로그 항을 제외한)를 입증하며, 이는 타원 문제의 최적 N-항 근사와 동일하며, 조건수는 일관되게 유계이고, 음의 순서의 소볼레프 안정성 요구 없이 적응형 기법을 통한 안정성을 확보한다.

ABSTRACT

We study linear parabolic initial-value problems in a space-time variational formulation based on fractional calculus. This formulation uses "time derivatives of order one half" on the bi-infinite time axis. We show that for linear, parabolic initial-boundary value problems on $(0,\infty)$, the corresponding bilinear form admits an inf-sup condition with sparse tensor product trial and test function spaces. We deduce optimality of compressive, space-time Galerkin discretizations, where stability of Galerkin approximations is implied by the well-posedness of the parabolic operator equation. The variational setting adopted here admits more general Riesz bases than previous work; in particular, no stability in negative order Sobolev spaces on the spatial or temporal domains is required of the Riesz bases accommodated by the present formulation. The trial and test spaces are based on Sobolev spaces of equal order $1/2$ with respect to the temporal variable. Sparse tensor products of multi-level decompositions of the spatial and temporal spaces in Galerkin discretizations lead to large, non-symmetric linear systems of equations. We prove that their condition numbers are uniformly bounded with respect to the discretization level. In terms of the total number of degrees of freedom, the convergence orders equal, up to logarithmic terms, those of best $N$-term approximations of solutions of the corresponding elliptic problems.

연구 동기 및 목표

• 선형 포아송 초기경계값 문제에 대해 안정적이고 압축 가능한 공간-시간 갈레르킨 이산화를 개발하기 위해.
• 이전 방법들이 공간 또는 시간 Riesz 기저의 음의 순서 소볼레프 안정성 요구를 필요로 하는 한계를 극복하기 위해.
• 공간-시간 영역에서 적응형 압축 갈레르킨 근사의 최적 수렴 속도를 확립하여, 해당 타원 문제의 최적 N-항 근사 속도와 일치시키기 위해.
• 희소 텐서곱 이산화에서 조건수의 균일 유계성을 입증하여 V*에서 공간 Riesz 기저 안정성에 의존하지 않기 위해.
• 유도된 큰 비대칭 선형 연립방정식의 적응형 웨이블릿 솔버의 s*-계산 가능성과 최적성 증명하기 위해.

제안 방법

• 시간 축이 양의 무한과 음의 무한으로 이어지는 이중 무한 시간 축에서 시간 도함수의 순서가 1/2인 분수적 미분을 기반으로 한 공간-시간 변분 형식을 사용하여 포아송 문제를 수립한다.
• 보흐너 공간과 보간 이론을 사용하여 고전적 형식 사이의 중간 해공간을 정의하여 문제의 잘 정의됨과 유일한 해 존재를 보장한다.
• 공간과 시간에서 다중 수준 웨이블릿 분해의 희소 텐서곱 공간을 사용하여 압축 가능한 갈레르킨 이산화를 구성한다.
• 공간 및 시간 변수 모두에 대해 H^{1/2} 공간에서 Riesz 기저를 사용한 페트로프-갈레르킨 형식의 이산 인프-서프 조건을 확립한다.
• 논문 [6]의 적응형 웨이블릿 기법을 적용하여 유도된 무한대 행렬 문제를 최적 복잡도로 해결하며, '적응형에 의한 안정성'을 활용한다.
• 웨이블릿 기저에 대한 온건한 가정 하에 시스템 행렬과 그 구성 요소(분수적 시간 도함수의 질량 및 강성 행렬)의 s*-계산 가능성을 증명한다.

실험 결과

연구 질문

• RQ1 동일한 순서의 시간 소볼레프 공간을 사용하는 대칭적 공간-시간 변분 형식이, Riesz 기저의 음의 순서 소볼레프 안정성 요구 없이 포아송 PDE의 안정적이고 최적의 갈레르킨 이산화를 보장할 수 있는가?
• RQ2 압축 가능한 공간-시간 갈레르킨 이산화의 수렴 속도는 해당 타원 문제의 최적 N-항 근사와 어떻게 비교되는가?
• RQ3 유도된 큰 비대칭 선형 연립방정식의 조건수는 희소 텐서곱 이산화 수준에 따라 균일하게 유계로 유지되는가?
• RQ4 어떻게 적응형 웨이블릿 기법을 사용하여 시스템 행렬과 그 구성 요소를 효율적이고 안정적으로 계산할 수 있는가?
• RQ5 적응형 웨이블릿 솔버로 이산 문제를 해결할 때의 계산 복잡도는 무엇이며, 부드러운 해에 대해 최적 수렴을 달성하는가?

주요 결과

• 제안된 공간-시간 변분 형식은 고전적 형식 사이의 중간 해공간인 보흐너 공간에서 유일한 해를 가지며, 문제의 잘 정의됨을 보장한다.
• 희소 텐서곱 갈레르킨 이산화에 대해 이산 인프-서프 조건이 성립하여, V*에서 Riesz 기저 안정성 요구 없이도 안정성을 보장한다.
• 압축 가능한 갈레르킨 이산화의 수렴 속도는 타원 문제의 최적 N-항 근사 속도와 로그 항을 제외한 범위에서 일치한다.
• 유도된 큰 비대칭 선형 연립방정식의 조건수는 이산화 수준에 관계없이 균일하게 유계이다.
• 적응형 웨이블릿 솔버는 최적 복잡도를 달성한다: 허용 오차 ε에 대해 필요한 연산 수와 저장 용량이 O(ε^(-1/s))로 스케일링되며, 이 상수는 공간 웨이블릿 기저의 Riesz 및 적합성 상수에만 의존한다.
• 웨이블릿 기저에 대한 온건한 가정 하에 시스템 행렬과 구성 요소(분수적 시간 도함수의 질량 및 강성 행렬)는 s*-계산 가능하며, 이는 효율적인 적응형 알고리즘에 의한 해를 가능하게 한다.