[논문 리뷰] Uniform convergence on a Bakhvalov-type mesh using the preconditioning approach: Technical report

Thái Anh Nhan, Relja Vulanović|arXiv (Cornell University)|2015. 04. 16.

Differential Equations and Numerical Methods참고 문헌 9인용 수 3

한 줄 요약

이 논문은 Bakhvalov 유형의 메esh에서 경계층을 가진 특수하게 섭동된 대류-확산 경계값 문제에 대해 균일 수렴을 달성하기 위한 전조건화 방법을 제시한다. 이산 시스템에 전조건자를 적용함으로써, 섭동 매개변수 ε에 관계없이 조건수(condition number)가 유계로 유지되며, 이는 ε에 대해 균일하게 O(N⁻¹)의 최적 순서의 점별 수렴을 가능하게 한다. 이는 경계층이 존재하는 상황에서도 성립한다.

ABSTRACT

The linear singularly perturbed convection-diffusion problem in one dimension is considered and its discretization on a Bakhvalov-type mesh is analyzed. The preconditioning technique is used to obtain the pointwise convergence uniform in the perturbation parameter.

연구 동기 및 목표

• 경계층을 가진 특수하게 섭동된 대류-확산 문제에 대해 균일 수렴하는 수치적 방법을 개발한다.
• 섭동 매개변수 ε이 작을 경우 정확도와 안정성을 유지하는 데 도전하는 문제를 다룬다.
• Bakhvalov 유형의 메esh에서 이산 시스템에 전조건화를 적용하여 ε에 대해 균일 수렴을 보장한다.
• 전조건화된 시스템의 조건수가 ε에 독립적으로 유계로 유지됨을 증명한다.
• 이산 해가 연속 해에 대해 ε에 대해 균일하게 최적 순서의 점별 수렴을 보임을 확립한다.

제안 방법

• 일차원 특수하게 섭동된 대류-확산 문제를 Bakhvalov 유형의 메esh에서 유한차분 업윈드 스킴으로 이산화한다.
• 이산 시스템의 조건수를 향상시키기 위해 스케일링 전조건자 행렬을 정의한다.
• 전조건자를 사용하여 원래 시스템을 조건수가 유계인 새로운 시스템으로 변환한다.
• 일致성-안정성 원리를 적용하여 일치 오차가 O(N⁻¹)임을 증명하고, 전조건화된 시스템의 유계로 유지되는 역행렬을 확보한다.
• 해를 매끄럽고 경계층 영역으로 분해하여 서로 다른 메쉬 영역에서의 절단 오차를 추정한다.
• 이산 최대원리(discrete maximum principle)를 통해 유계 조건수와 일치 오차 추정치를 조합하여 균일 수렴을 확립한다.

실험 결과

연구 질문

• RQ1 전조건화는 Bakhvalov 유형의 메esh에서 대류-확산 문제에 대해 균일 수렴을 달성하는 데 사용될 수 있는가?
• RQ2 전조건화를 사용할 경우 이산 시스템의 조건수가 섭동 매개변수 ε에 독립적으로 유지되는가?
• RQ3 Bakhvalov 유형의 메esh에서의 유한차분 스킴의 일치 오차는 ε에 대해 균일하게 O(N⁻¹)인가?
• RQ4 이 방법을 통해 해를 ε에 대해 균일하게 최적 순서 O(N⁻¹)의 정확도로 근사할 수 있는가?
• RQ5 전조건자는 경계층 존재 시 이산 시스템의 안정성과 조건수에 어떤 영향을 미치는가?

주요 결과

• 전조건화된 시스템의 조건수는 ε에 관계없이 O(N²)로 유계로 유지된다.
• Bakhvalov 유형의 메esh에서의 유한차분 스킴의 일치 오차는 모든 ε에 대해 균일하게 O(N⁻¹)로 유계로 유지된다.
• a > 4/β를 만족할 경우, 이산 해는 정확한 해로 향해 점별 오차 O(N⁻¹)로 수렴하며, 이는 ε에 대해 균일하다.
• 전조건자는 작은 ε에 의해 유도되는 악조건화를 효과적으로 완화하여 강건한 수렴을 보장한다.
• 분석 결과, Bakhvalov 유형의 메쉬와 전조건화의 조합은 날카로운 경계층이 존재하는 상황에서도 최적 수렴 속도를 달성함을 확인한다.
• 메쉬와 전조건자 설계 덕분에 경계층 영역에 특별한 처리 없이도 균일 수렴을 달성한다.