[논문 리뷰] A space-time Trefftz discontinuous Galerkin method for the linear Schr\"odinger equation

Sergio Gómez, Andrea Moiola|arXiv (Cornell University)|2021. 06. 08.

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

한 줄 요약

이 논문은 조각별로 일정한 포텐셜을 가진 선형 슈뢰딩거 방정식을 위한 공간-시간 트레프츠 불연속 갈레르킨(DG) 방법을 제안한다. 비다항 복소 지수 기저 함수를 사용하여 국소적으로 방정식을 만족시키며, 1차원 및 2차원 문제에서 최적의 h수렴율을 달성한다. 다항식 DG 방법에 비해 자유도를 크게 감소시키면서도 안정성과 잘 정의됨을 유지한다.

ABSTRACT

A space-time Trefftz discontinuous Galerkin method for the Schr\"odinger equation with piecewise-constant potential is proposed and analyzed. Following the spirit of Trefftz methods, trial and test spaces are spanned by non-polynomial complex wave functions that satisfy the Schro\"odinger equation locally on each element of the space-time mesh. This allows for a significant reduction in the number of degrees of freedom in comparison with full polynomial spaces. We prove well-posedness and stability of the method, and, for the one- and two- dimensional cases, optimal, high-order, h-convergence error estimates in a skeleton norm. Some numerical experiments validate the theoretical results presented.

연구 동기 및 목표

• 조각별 일정한 포텐셜을 가진 시간에 의존하는 선형 슈뢰딩거 방정식을 위한 새로운 공간-시간 트레프츠-DG 방법을 개발한다.
• 국소적으로 방정식을 만족하는 비다항 기저 함수를 사용하여 자유도 수를 감소시킨다.
• 1차원 및 2차원 문제에서 이론적 안정성, 잘 정의됨, 최적의 h수렴 오차 추정치를 확립한다.
• 정사각형 우물 및 가우시안 포텐셜 문제에 대한 수치 실험을 통해 이론적 결과를 검증한다.
• 기저 매개변수 조정을 통한 p수렴 및 효율성 향상 잠재력을 탐색한다.

제안 방법

• 공간-시간 메esh에 다각형 요소를 사용하며 포텐셜 불연속성에 따라 정렬하고, 시간에는 업윈드 수치 플럭스, 공간에는 중심 플럭스를 사용하는 불연속 갈레르킨 설정을 적용한다.
• 시험 함수 및 테스트 함수는 국소적으로 슈뢰딩거 방정식을 만족하는 복소 지수 함수로 구성되며, 트레프츠 공간을 형성한다.
• 이산 트레프츠 공간는 정확한 해의 타일러 다항식을 차수 p까지 재현할 수 있도록 설계되어 최적의 근사 성질을 보장한다.
• 오차 분석에 스켈레톤 노름을 사용하며, 트레프츠 함수의 구조를 활용하여 역추정이 필요 없도록 한다.
• 수치 플럭스는 평균과 점프를 사용하며, 안정성과 일致성을 확보하기 위해 복소수 벌점 항을 적용한다.
• 분석은 이산 트레프츠 공간가 국소적으로 정확한 해의 타일러 다항식을 근사하는 데 필요한 조건을 맺는 데 기초한다.

실험 결과

연구 질문

• RQ1 조각별 일정한 포텐셜을 가진 선형 슈뢰딩거 방정식에 대해 공간-시간 트레프츠-DG 방법이 최적의 h수렴율을 달성할 수 있는가?
• RQ2 비다항 트레프츠 기저 함수의 사용이 표준 다항식 DG 방법에 비해 자유도를 얼마나 감소시키는가?
• RQ3 이산 트레프츠 공간에 어떤 조건이 있어야 정확한 해의 최적 국소 근사가 보장되는가?
• RQ4 이 방법은 p수렴율을 달성할 수 있으며, 다항식 DG 스킴에 비해 효율성이 높은가?
• RQ5 정확도는 기저 함수 매개변수(예: 파장 수와 방향)의 선택에 얼마나 민감한가?

주요 결과

• 정리 4.12에 의해 증명된 linemesh-스켈레톤 노름에서 1차원 및 2차원 문제 모두에서 최적의 h수렴율 hp를 달성한다.
• (1+1)-차원 정사각형 우물 문제에서 p=1일 때 약 2.07, p=2일 때 약 4.09, p=3일 때 약 6.49의 수렴율을 보이며, 최적의 h수렴을 확인한다.
• (2+1)-차원 가우시안 문제에서는 DG 노름에서 O(h^p) 수렴, 최종 시간에서 L2 노름에서 O(h^{p+1}) 수렴을 달성하여 이론적 예측과 일치한다.
• p수렴 분석 결과 오차는 O(e^{-b\sqrt{\#DOF}})로 감소하며, 다항식 DG 방법의 일반적인 O(e^{-b\sqrt[3]{\#DOF}})보다 훨씬 빠르게 감소한다.
• 수치 실험 결과 기저의 파장 수를 {−k*, 0, k*}로 조정하면 최적의 수렴율을 회복하지만, 부적절한 매개변수 선택은 열등한 성능을 초래한다.
• 이 방법은 역추정에 의존하지 않는 독립적인 오차 경계를 갖는다. 이는 표준 DG 설정에 비해 핵심적인 이점이다.