[논문 리뷰] Superconvergence of Discontinuous Galerkin method for linear hyperbolic equations

Waixiang Cao, Zhimin Zhang외 1명|arXiv (Cornell University)|2013. 11. 27.

Advanced Numerical Methods in Computational Mathematics참고 문헌 11인용 수 8

한 줄 요약

이 논문은 선형 초구역 방정식에 대한 불연속 갈레르킨(DG) 방법에서 하류 점 및 도메인 평균에서 오랫동안 추측되어 온 $2k+1$ 초수렴 속도를 증명하며, 새로운 보정 함수 접근법을 사용한다. 또한, 내부 왼쪽 라다우 점에서 도함수 근사에 대해 처음으로 $k+1$ 초수렴 속도를 확립하며, 적절한 초기 이산화 조건 하에서 수치적 검증을 통해 모든 이론적 결과를 확인한다.

ABSTRACT

In this paper, we study superconvergence properties of the discontinuous Galerkin (DG) method for one-dimensional linear hyperbolic equation when upwind fluxes are used. We prove, for any polynomial degree $k$, the $2k+1$th (or $2k+1/2$th) superconvergence rate of the DG approximation at the downwind points and for the domain average under quasi-uniform meshes and some suitable initial discretization. Moreover, we prove that the derivative approximation of the DG solution is superconvergent with a rate $k+1$ at all interior left Radau points. All theoretical finding are confirmed by numerical experiments.

연구 동기 및 목표

• 선형 초구역 방정식에 대한 DG 방법에서 하류 점 및 도메인 평균에서 오랫동안 열려 있던 $2k+1$ 초수렴 속도를 증명하는 데 목적이 있다.
• 내부 왼쪽 라다우 점에서 도함수 근사에 대해 처음으로 $k+1$ 초수렴 속도를 확립하는 데 목적이 있다.
• 초밀착성(초밀착성)을 달성하기 위해 보정 함수를 사용한 엄밀한 수학적 프레임워크를 제공하는 데 목적이 있다.
• 특히 $k > 3$일 경우, $2k+1$ 초수렴 속도를 달성하기 위해 보정 절차에 기반한 적절한 초기 이산화 조건이 필수적임을 보여주는 데 목적이 있다.
• 기존의 초수렴 결과를 통합하고 확장하여 일차원 선형 초구역 시스템에서의 초수렴 성질에 대한 완전한 그림을 제시하는 데 목적이 있다.

제안 방법

• 모든 하류 점에서 값이 0이 되는 보정 함수 $w$를 구성하여, DG 해와 $P_h^-u - w$ 사이의 초밀착성을 확보한다. 여기서 $P_h^-u$는 정확한 해를 하류 점에서 보간하는 절단된 라다우 투영이다.
• 분석은 에너지 추정과 하류 점에서 정확한 해를 보간하는 특수 투영 $P_h^-$에 기반하며, 오차 구조에서 고차 근사의 가능성을 보장한다.
• 시간 적분을 정확하게 가정하고, 이전 연구들과는 다름없는 새로운 분석 경로를 사용하며, 보정 함수가 초구역 문제에서 오차 전파를 제어하는 데서의 역할을 강조한다.
• 이 방법은 DG 해와 보정된 라다우 투영 사이의 초밀착성을 순서 $2k+1$로 확립하여, 하류 점 및 도메인 평균에서의 초수렴을 이끈다.
• 내부 왼쪽 라다우 점에서의 도함수 초수렴은 동일한 초밀착성 프레임워크에서 유도되며, DG 형식의 수직성과 안정성 성질을 활용한다.
• 다양한 초기 이산화 방법을 사용한 수치 실험을 통해 이론적 수렴 속도를 검증하며, $k=4$일 때 유일하게 제안된 보정 기반 방법만이 $2k+1$ 속도를 달성하는 것을 보여주는 비교 분석도 포함된다.

실험 결과

연구 질문

• RQ1 선형 초구역 방정식에 대한 DG 방법에서 하류 점에서의 $2k+1$ 초수렴 속도를 엄밀히 증명할 수 있는가?
• RQ2 DG 해의 도함수 근사는 내부 왼쪽 라다우 점에서 초수렴하는가? 만약 그렇다면, 어떤 속도로 초수렴하는가?
• RQ3 특히 $k > 3$일 경우, 높은 차수의 초수렴을 달성하기 위해 초기 이산화 조건이 어떤 역할을 하는가?
• RQ4 동일한 프레임워크 하에서 도메인 평균에 대한 $2k+1$ 초수렴 속도를 증명할 수 있는가?
• RQ5 제안된 보정 함수 접근법은 기존 방법과 비교해 최적의 초수렴 속도를 달성하는 데 어떻게 다른가?

주요 결과

• 논문은 DG 방법에서 하류 점에서 평균적으로 $2k+1$ 초수렴 속도를 증명하며, 선형 초구역 방정식에 대한 오랫동안 추측되어 온 결과를 확인한다.
• 도메인 평균에 대해서도 $2k+1$ 초수렴 속도를 확립하였으며, 이는 이전에 알려진 결과보다 $1/2$ 차수 높은 결과이다.
• 내부 왼쪽 라다우 점에서 DG 해의 도함수 근사에 대해 처음으로 $k+1$ 초수렴 속도를 증명한다.
• 수치 실험을 통해 이론적 수렴 속도를 확인하였으며, $k=4$일 때는 오직 제안된 보정 기반 초기 이산화 방법만이 $2k+1$ 초수렴 속도를 달성한다.
• 하류 점에서의 $2k+1/2$ 초수렴 속도는 증명 상에서 최적화되지 않아, 수치적 속도는 $2k+1$이지만 이론적 추정치와의 격차가 존재함을 시사한다.
• 셀 평균에 대한 $2k+1$ 초수렴 속도의 증명은 강력한 수치적 증거가 있음에도 불구하고 여전히 열린 문제로 남아 있다.