[논문 리뷰] An analysis of the practical DPG method

Jay Gopalakrishnan, Weifeng Qiu|arXiv (Cornell University)|2011. 07. 21.

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

한 줄 요약

이 논문은 라플라스 방정식과 선형 탄성에 대한 실용적 비연속 갈레르킨(DPG) 방법의 완전한 오차 분석을 제공하며, 시험 함수 공간의 다항식 차수 $ r \geq p + N $일 때 시험 함수의 계산에 사용되는 시도-시험 연산자 $ T $ 의 근사가 최적 수렴 속도를 유지함을 보여준다. 여기서 $ p $ 는 시도 다항식의 차수이고 $ N $ 은 공간 차원이다. 분석은 메시의 이산화, 시험 함수의 이산화, 강성 행렬의 조건수 등 메서드의 모든 실용적 근사에 대해 고려한다.

ABSTRACT

In this work we give a complete error analysis of the Discontinuous Petrov Galerkin (DPG) method, accounting for all the approximations made in its practical implementation. Specifically, we consider the DPG method that uses a trial space consisting of polynomials of degree $p$ on each mesh element. Earlier works showed that there is a "trial-to-test" operator $T$, which when applied to the trial space, defines a test space that guarantees stability. In DPG formulations, this operator $T$ is local: it can be applied element-by-element. However, an infinite dimensional problem on each mesh element needed to be solved to apply $T$. In practical computations, $T$ is approximated using polynomials of some degree $r > p$ on each mesh element. We show that this approximation maintains optimal convergence rates, provided that $r\ge p+N$, where $N$ is the space dimension (two or more), for the Laplace equation. We also prove a similar result for the DPG method for linear elasticity. Remarks on the conditioning of the stiffness matrix in DPG methods are also included.

연구 동기 및 목표

• 실용적 DPG 방법에 대한 완전한 오차 분석을 제공하여, 구현 과정에서 발생하는 모든 근사를 고려한다.
• 시도-시험 연산자 $ T $ 가 고차수 다항식을 사용해 근사될 경우 DPG 방법의 최적 수렴 속도를 확립한다.
• 실용적 DPG 설정에서 강성 행렬의 조건수를 분석한다.
• 이deal DPG 방법의 이론적 프레임워크를 완전히 구현 가능한 실용적 변형으로 확장한다.
• 두 번째 차수 탄성 문제, 특히 선형 탄성 문제에 대해 메서드의 강건성과 수렴 성질을 검증한다.

제안 방법

• 실용적 DPG 방법은 무한차원 시도-시험 연산자 $ T $ 를 유한차원 시험 공간의 다항식 차수 $ r $ 에 대한 갈레르킨 투영을 통해 정의된 이산 연산자 $ T^r $ 으로 근사함으로써 구성된다.
• 국소적으로 시험 함수 계산이 가능하도록, 혼합된 미지수(유량, 경계값, 수치적 경계값)를 사용하는 초우월한 변분 형식을 사용한다.
• 이산 시험 공간 $ V_h^r $ 는 $ T^r(U_h) $ 으로 정의되며, 여기서 $ U_h $ 는 차수 $ p $ 의 다항식으로 이루어진 시도 공간이고, $ T^r $ 은 각 요소에서 국소 해를 통해 계산된다.
• 핵심 기술 도구로는 변분 형식 내의 미분 연산자와 가환하는 발산을 유지하는 프로젝션 연산자 $ \varPi_{p+2}^{(\mathrm{div},\mathbb{S})} $ 이 사용된다.
• 분석은 새로운 보조정리들을 활용하여 안정성과 일致성 조건을 검증하며, 특히 잔여항이 이산 시험 공간과 수직임을 보장한다.
• 시험 공간의 비연속 갈레르킨 구조로 인해 강성 행렬이 블록 대각형임을 보여주며, 이는 국소 역행렬 계산을 통한 효율적 해법을 가능하게 한다.

실험 결과

연구 질문

• RQ1 시도-시험 연산자 $ T $ 가 다항식 차수 $ r $ 로 근사될 경우 실용적 DPG 방법의 수렴 속도는 어떻게 되는가?
• RQ2 시도-시험 연산자 $ T $ 가 $ T^r $ 으로 근사될 때, 이deal DPG 방법의 최적 수렴 성질이 유지되는가?
• RQ3 라플라스 방정식에서 최적 수렴 속도를 유지하기 위해 필요한 최소 다항식 차수 $ r $ 는 얼마인가?
• RQ4 특히 $ T $ 가 근사될 경우 실용적 DPG 방법에서 강성 행렬의 조건수는 어떻게 행동하는가?
• RQ5 이deal DPG 방법의 이론적 프레임워크는 선형 탄성 문제에 대해 완전히 이산화된 실용적 변형으로 확장될 수 있는가?

주요 결과

• 라플라스 방정식의 경우, $ r \geq p + N $ 를 만족할 경우 실용적 DPG 방법에서 최적 수렴 속도가 유지된다. 여기서 $ N $ 은 공간 차원이다.
• 실용적 DPG 방법의 오차 추정은 준최적이며, 오차는 시도 공간 내 최선의 근사 오차로 유계이다.
• 분석 결과, $ r \geq p + N $ 를 만족할 경우, 시도-시험 연산자 $ T $ 의 근사 $ T^r $ 이 수렴 속도를 떨어뜨리지 않음을 확인하였다.
• 초우월한 변분 형식과 하이브리드화 덕분에, 각 요소에서 국소 해를 통해 시험 함수를 계산하더라도 메서드는 안정성과 최적 수렴 성질을 유지한다.
• 실용적 DPG 방법의 강성 행렬은 블록 대각형이므로, 효율적인 해법 절차와 유리한 조건수를 보장한다.
• 선형 탄성 문제에 대해서도 동일한 조건 $ r \geq p + N $ 하에서 유사한 준최적 오차 추정이 확립되며, 이는 탄성 계수 텐서가 요소별로 일정하다는 가정 하에 성립한다.