[논문 리뷰] Cross-Points in Domain Decomposition Methods with a Finite Element Discretization

Martin J. Gander, Kévin Santugini-Repiquet|arXiv (Cornell University)|2014. 04. 18.

Advanced Numerical Methods in Computational MathematicsEngineering참고 문헌 15인용 수 16

한 줄 요약

이 논문은 유한요소 기반 비중첩 도메인 분할 방법에서 두 개 이상의 하위영역이 만날 때 발생하는 교차점에서 노이만 조건을 일관되게 이산화하는 문제를 다룬다. 두 가지 방법을 제안한다: 추가 변수 방법은 안정적인 수렴을 위해 추가적인 경계 변수를 도입하고, 완전한 통신 방법은 교차점에 공유되는 모든 하위영역이 참여함으로써 일관된 노이만 교환을 보장한다. 이 두 방법 모두 수렴 증명과 최적의 매개변수 선택을 가능하게 하며, 질량 럼핑이 대수적으로 고차원의 벤츠셀 유사 전달 조건을 생성함을 보여준다.

ABSTRACT

Non-overlapping domain decomposition methods necessarily have to exchange Dirichlet and Neumann traces at interfaces in order to be able to converge to the underlying mono-domain solution. Well known such non-overlapping methods are the Dirichlet-Neumann method, the FETI and Neumann-Neumann methods, and optimized Schwarz methods. For all these methods, cross-points in the domain decomposition configuration where more than two subdomains meet do not pose any problem at the continuous level, but care must be taken when the methods are discretized. We show in this paper two possible approaches for the consistent discretization of Neumann conditions at cross-points in a Finite Element setting.

연구 동기 및 목표

• 유한요소 도메인 분할 방법에서 교차점에서의 노이만 조건 이산화의 비일관성을 해결하기 위해.
• 여러 하위영역이 만날 때 노이만 트레이스에 대해 수학적으로 일관되고 수렴 가능한 이산화 전략을 제공하기 위해.
• 교차점 특이성이 존재하더라도 비중첩 최적화 샤우츠 방법의 수렴성을 입증하기 위해.
• 질량 럼핑을 통해 대수적으로 고차원의 벤츠셀 전달 조건을 생성할 수 있는지 보여주기 위해, 복잡한 탄성 도함수 이산화를 피하기 위해.
• 교차점이 존재할 경우 수렴 진동을 최소화하는 최적의 루빈 매개변수와 오버럼핑 인자 식별하기 위해.

제안 방법

• 보조 변수 방법을 도입하여 하위영역이 노이만 자료를 표현하기 위해 추가적인 경계 변수를 저장함으로써 교차점에서 일관된 노이만 트레이스 이산화를 가능하게 한다.
• 완전한 통신 방법을 적용하여 교차점에 연결된 모든 하위영역(포인트로만 연결된 경우도 포함)과의 정보 교환을 요구함으로써 일관성을 유지한다.
• 에너지 추정을 사용하여 보조 변수를 사용할 경우 도메인 분할 알고리즘의 수렴성을 증명하며, 이는 기존의 유한요소 이산화 방식으로는 교차점에서 달성할 수 없는 결과이다.
• 교차점에서 여러 가능한 노이만 트레이스 분할 방식 중에서 반복 과정에서의 진동을 최소화하는 방법을 제안한다.
• 루빈 조건의 질량 럼핑이 대수적으로 고차원의 벤츠셀 전달 조건을 동일하게 생성함을 보여주며, 이는 교차점에서도 탄성 도함수의 직접 이산화가 필요 없음을 의미한다.
• 2×2 하위영역과 $Q_1$ 유한요소를 사용하여 라플라스 방정식에 대한 수치 실험을 수행하여 다양한 루빈 매개변수와 오버럼핑 수준에서 수렴 인자를 평가한다.

실험 결과

연구 질문

• RQ1 유한요소 도메인 분할 방법에서 교차점에서의 노이만 조건을 어떻게 일관되게 이산화할 수 있는가?
• RQ2 교차점이 존재할 경우 비중첩 최적화 샤우츠 방법의 수렴성을 어떻게 증명할 수 있는가?
• RQ3 반복 과정에서의 진동을 최소화하기 위해 교차점에서의 노이만 트레이스를 어떻게 분할하는 것이 최적인가?
• RQ4 유한요소 맥락에서 질량 럼핑을 통해 루빈 조건으로부터 대수적으로 고차원의 벤츠셀 전달 조건을 얻을 수 있는가?
• RQ5 교차점이 존재할 경우 최적의 루빈 매개변수와 오버럼핑 인자는 메쉬 정밀도 향상에 따라 어떻게 척도화되는가?

주요 결과

• 보조 변수 방법은 에너지 추정을 통한 엄밀한 수렴 증명을 가능하게 하며, 이는 기존의 유한요소 이산화 방식으로는 교차점에서 달성할 수 없는 결과이다.
• 완전한 통신 방법은 교차점에 공유되는 모든 하위영역(포인트로만 연결된 경우도 포함)을 포함함으로써 일관된 노이만 조건 이산화를 보장한다.
• $100\times100$ 메쉬(하위영역당)에서 완전한 통신 방법을 사용할 경우 $p=2.0$ 및 $\omega=122.5$ 조건에서 최고의 수렴 인자 $\kappa=0.6013464$ 를 달성하였다.
• 보조 변수 방법을 사용할 경우 $100\times100$ 셀에서 $p=2.0$ 및 $\omega=122.0$ 조건에서 최고의 수렴 인자 $\kappa=0.6006753$ 를 확보하였다.
• 최적의 루빈 매개변수 $p$ 는 교차점에서 약 $O(1/h)$ 비율로 척도화되며, 이는 이전의 이론적 추정과 일치한다. 이는 일반적인 경계점에서의 $O(1/\sqrt{h})$ 와는 다름을 나타낸다.
• 루빈 조건의 질량 럼핑은 교차점에서도 탄성 도함수의 이산화가 필요 없이 대수적으로 고차원의 벤츠셀 전달 조건을 동일하게 생성한다.