[논문 리뷰] Domain decomposition algorithms for two dimensional linear Schrödinger equation

Christophe Besse, Feng Xing|arXiv (Cornell University)|2015. 06. 18.

Electromagnetic Simulation and Numerical Methods참고 문헌 15인용 수 4

한 줄 요약

이 논문은 일정한 잠재력과 일반 잠재력에 대해 각각 2차원 선형 슈뢰딩거 방정식을 해결하기 위한 새로운 도메인 분할 알고리즘과 조건화 알고리즘을 제안한다. 고전적인 로빈 조건 대신 새로운 투명 흡수 경계 조건을 도입함으로써, 수렴 속도가 크게 향상되고 계산 시간이 감소한다. 수치 결과에 따르면, 대규모 하위영역 수에서 총 계산 시간이 최대 80% 감소하는 것으로 나타났다.

ABSTRACT

This paper deals with two domain decomposition methods for two dimensional linear Schr{ö}dinger equation, the Schwarz waveform relaxation method and the domain decomposition in space method. After presenting the classical algorithms, we propose a new algorithm for the free Schr{ö}dinger equation and a preconditioned algorithm for the general Schr{ö}dinger equation. These algorithms are studied numerically, which shows that the two new algorithms could accelerate the convergence and reduce the computation time. Besides the traditional Robin transmission condition, we also propose to use a newly constructed absorbing condition as the transmission condition.

연구 동기 및 목표

• 일般 잠재력이 있는 2차원 선형 슈뢰딩거 방정식에 대한 도메인 분할 방법의 수렴 속도와 계산 효율성을 향상시키기 위해.
• 하위영역 수가 많아질수록 수렴 속도가 느려지는 고전적 샤우르 웨이브폼 리커리션(SWR) 및 공간 기반 도메인 분할(DDS) 방법의 문제를 해결하기 위해.
• 일정한 잠재력의 경우 새로운 알고리즘과 일반 잠재력의 경우 새로운 흡수 경계 조건을 사용한 조건화 알고리즘을 개발하기 위해.
• 다양한 메쉬 크기와 하위영역 수에 걸쳐 제안된 알고리즘의 강건성과 확장성을 검증하기 위해.
• 기존의 로빈 조건보다 새로운 경계 조건이 반복 횟수와 계산 시간 측면에서 우수한 성능을 보임을 입증하기 위해.

제안 방법

• 고정 잠재력이 있는 슈뢰딩거 방정식에 대한 새로운 도메인 분할 알고리즘을 제안하며, 새로 구성한 투명 흡수 경계 조건을 사용한다.
• 변동 잠재력이 있는 일반 슈뢰딩거 방정식에 대해 동일한 흡수 조건을 경계 조건으로 사용하는 조건화 알고리즘을 개발한다.
• 샤우르 웨이브폼 리커리션(SWR) 방법과 공간 기반 도메인 분할(DDS) 방법을 적용하며, 시간-공간 도메인 분할과 하위영역 간 반복적 갱신을 수행한다.
• 경계 문제를 유도하고 분석하여 새로운 알고리즘의 수렴 성질을 정당화한다.
• DDS 방법에서는 반명시적 시간 이산화를 사용하며, 반복 연산자의 스펙트럼 특성을 바탕으로 조건자(preconditioner)를 구성한다.
• 최대 1024개의 하위영역를 사용한 병렬 계산을 통해 다양한 메쉬 해상도에서 확장성과 성능을 테스트한다.

실험 결과

연구 질문

• RQ1 새로운 흡수 경계 조건이 2차원 선형 슈뢰딩거 방정식에 대한 도메인 분할 방법의 수렴 속도를 향상시킬 수 있는가?
• RQ2 제안된 조건화 알고리즘이 고전적 SWR 및 DDS 방법에 비해 반복 횟수와 계산 시간 측면에서 어떻게 비교되는가?
• RQ3 새로운 알고리즘이 하위영역 수와 메쉬 해상도의 변화에 대해 강건성과 확장성을 유지하는가?
• RQ4 기존의 로빈 조건에 비해 새로운 알고리즘의 경계 조건 파라미터 선택에 대해 얼마나 민감한가?
• RQ5 제안된 방법들이 단일 프로세서에서 전체 문제를 푸는 데 필요한 기준 시간 이하의 계산 시간을 달성할 수 있는가?

주요 결과

• 32개의 하위영역(N=32), Δx=1/128, Δy=1/8 조건에서 조건화 알고리즘이 반복 횟수를 17회에서 5회로 줄여 수렴 속도를 크게 향상시켰다.
• 동일한 설정에서 총 계산 시간은 고전적 방법의 12.1초에서 조건화 방법의 7.3초로 감소하여 40% 감소했다.
• 더 세밀한 메쉬(Δx=1/2048, Δy=1/64) 조건에서, 조건화 알고리즘은 하위영역 수 N=512일 때도 반복 횟수(m=5–6)를 낮게 유지하고, 계산 시간은 270초 이내로 유지했다.
• 고전적 방법은 하위영역 수가 증가함에 따라 반복 횟수가 증가(하나의 하위영역에서 N=256일 때 8회에서 N=512일 때 12회로 증가)하여 확장성이 떨어지는 것으로 나타났으며, 반면 조건화 방법은 안정성을 유지했다.
• N≥16일 때 조건화 알고리즘이 기준 시간(Tref=16.1초) 이하의 계산 시간을 달성하여 강력한 확장성과 효율성 향상을 입증했다.
• 새로운 흡수 경계 조건은 로빈, 패드(Padé) 등 다양한 경계 조건 유형과 파라미터에 대해 강건하며, 튜닝에 대해 민감하지 않음을 보였다.