[논문 리뷰] A Reynolds-robust preconditioner for the Reynolds-robust Scott-Vogelius discretization of the stationary incompressible Navier-Stokes equations

Patrick E. Farrell, Lawrence Mitchell외 2명|arXiv (Cornell University)|2020. 04. 20.

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

한 줄 요약

이 논문은 삼등분 메시로 정렬된 메시에서 비압축성 스토크스-바울레 방정식의 스코트-보겔리우스 유한요소 해법에 대해 보강된 라그랑주 승수 전처리기를 소개한다. 메시의 구조를 활용하여 기울기-발산 항의 핵을 분해함으로써, 메시의 구조를 활용하여 기울기-발산 항의 핵을 분해함으로써, 레이놀즈 수에 관계없이 수렴 성능과 오차 추정치가 둘 다 레이놀즈 수에 영향을 받지 않는 성능을 달성하였으며, 2차원 및 3차원 수치 실험을 통해 검증되었다.

ABSTRACT

Augmented Lagrangian preconditioners have successfully yielded Reynolds-robust preconditioners for the stationary incompressible Navier-Stokes equations, but only for specific discretizations. The discretizations for which these preconditioners have been designed possess error estimates which depend on the Reynolds number, with the discretization error deteriorating as the Reynolds number is increased. In this paper we present an augmented Lagrangian preconditioner for the Scott-Vogelius discretization on barycentrically-refined meshes. This achieves both Reynolds-robust performance and Reynolds-robust error estimates. A key consideration is the design of a suitable space decomposition that captures the kernel of the grad-div term added to control the Schur complement; the same barycentric refinement that guarantees inf-sup stability also provides a local decomposition of the kernel of the divergence. The robustness of the scheme is confirmed by numerical experiments in two and three dimensions.

연구 동기 및 목표

• 고레이놀즈 수에서 오차 추정치가 악화되는 문제를 해결하기 위해 스코트-보겔리우스 해법에 대해 레이놀즈 수에 강건한 전처리기가 부족한 문제를 해결한다.
• 레이놀즈 수에 영향을 받지 않는 수렴 속도를 유지하는 전처리기를 개발한다.
• 레이놀즈 수가 증가함에 따라 해석 오차가 안정적이고 잘 제어됨을 보장한다.
• 삼등분 정련을 활용하여 기울기-발산 항의 핵에 대한 국소적 공간 분해를 가능하게 한다.
• 고레이놀즈 수 조건에서도 수치적 강건성과 최적의 오차 추정치를 동시에 확보한다.

제안 방법

• 스웨어드-포인트 시스템에서 발생하는 슈어 컴플리멘트를 안정화하기 위해 보강된 라그랑주 프레임워크를 활용한다.
• 삼등분 정련 기반으로 특수화된 공간 분해를 구성하여 기울기-발산 항의 핵을 포착한다.
• 정련이 인프-서프 안정성을 보장하고, 비압축성 핵의 국소적 분해를 가능하게 한다.
• 전처리기는 반복적 해법 성능과 오차 추정치 양쪽 모두에서 강건성을 유지하도록 설계된다.
• 메시의 기하학적 구조를 활용하여 핵의 효과적인 근사가 가능하다.
• 이론적 주장의 검증을 위해 2차원 및 3차원에서 수치 실험을 수행한다.

실험 결과

연구 질문

• RQ1 스코트-보겔리우스 해법에 대해 보강된 라그랑주 전처리기를 설계할 수 있는가? 이는 레이놀즈 수에 관계없이 수렴 성능을 보장할 수 있는가?
• RQ2 삼등분 정련의 사용이 기울기-발산 항의 핵에 대한 안정적이고 강건한 공간 분해를 가능하게 하는가?
• RQ3 해당 방법이 레이놀즈 수에 강건한 반복적 성능과 오차 추정치를 동시에 확보할 수 있는가?
• RQ4 2차원 및 3차원에서 고레이놀즈 수 흐름에서 전처리기는 어떻게 성능을 발휘하는가?
• RQ5 삼등분 정련의 어떤 구조적 특성이 레이놀즈 수에 강건한 성능을 확보하기 위해 필요한 핵 분해를 가능하게 하는가?

주요 결과

• 제안된 전처리기는 레이놀즈 수에 관계없이 선형 시스템의 반복적 해법에 대해 레이놀즈 수에 강건한 수렴 성능을 달성한다.
• 레이놀즈 수가 증가함에 따라 오차 추정치가 악화되지 않으며, 이는 강건한 오차 추정치 덕분에 유지된다.
• 삼등분 정련은 기울기-발산 항의 핵을 효과적으로 포착하는 국소적 공간 분해를 가능하게 한다.
• 수치 실험을 통해 2차원 및 3차원 모두에서 전처리기가 강건한 성능을 유지함을 확인하였다.
• 보강된 라그랑주 안정화와 메시 기반 핵 분해의 조합이 안정성과 정확성을 모두 보장한다.
• 기존 방법이 고레이놀즈 수에서 비강건하거나 오차 제어가 부족한 문제를 성공적으로 해결하였다.