[논문 리뷰] Robust preconditioning of monolithically coupled multiphysics problems

Karl Erik Holter, Miroslav Kuchta외 1명|arXiv (Cornell University)|2020. 01. 15.

Advanced Numerical Methods in Computational Mathematics참고 문헌 53인용 수 6

한 줄 요약

이 논문은 분수 및 가중 치환 라플라시안 연산자를 인터페이스에서 사용하여 단일화된 결합 다물리컬 문제에 대해 강건하고 매개변수에 강건한 조건부 프레임워크를 제안한다. 분수 스오볼레프 공간에서의 연산자 조건부를 활용함으로써, Darcy-Stokes 및 유체-구조 상호작용과 같은 도전적인 문제들에 대해 메esh에 의존하지 않는 수렴성을 달성하였으며, 수치 실험을 통해 다양한 물리적 매개변수와 메쉬 정밀도에서 최적의 수렴 속도와 유한한 조건수를 확인하였다.

ABSTRACT

In many applications, one wants to model physical systems consisting of two different physical processes in two different domains that are coupled across a common interface. A crucial challenge is then that the solutions of the two different domains often depend critically on the interaction at the interface and therefore the problem cannot be easily decoupled into its subproblems. Here, we present a framework for finding robust preconditioners for a fairly general class of such problems by exploiting operators representing fractional and weighted Laplacians at the interface. Furthermore, we show feasibility of the framework for two common multiphysics problems; namely the Darcy-Stokes problem and a fluid--structure interaction problem. Numerical experiments that demonstrate the effectiveness of the approach are included.

연구 동기 및 목표

• 서로 다른 물리적 하위문제들이 인터페이스에서 강하게 결합된 경우에 대해 단일화된 형식의 다물리컬 문제에 대한 강건한 조건부 전략을 개발하기.
• 점성도, 투과도, Beavers-Joseph-Saffman 계수와 같은 물리적 매개변수의 변화에 대해 조건부가 매개변수에 강건하도록 보장하기.
• 가중 및 분수 스오볼레프 공간을 사용하여 인터페이스 결합 문제에 대한 연산자 조건부 이론을 확장하기.
• Darcy-Stokes 유동 및 유체-구조 상호작용과 같은 두 가지 표준 문제에 대해 프레임워크를 수치적으로 검증하기.
• 메쉬 정밀도 증가 조건에서도 조건수는 유한하고 최적의 수렴 속도를 유지하는지 확인하기.

제안 방법

• 인터페이스 연속성 조건을 강제하기 위해 라그랑주 승수를 사용하는 단일화된 약한 형태로 결합된 다물리컬 문제를 수립하기.
• 분수 및 가중 스오볼레프 공간에서의 리에츠 사상에 기반한 조건부를 구성하여 인터페이스 변수의 올바른 척도를 반영하기.
• 분수 라플라시안 연산자를 인터페이스에서 사용하여 전달 조건을 모델링하고, Beavers-Joseph-Saffman 조건을 포함하기.
• 속도, 압력 및 라그랑주 승수에 대한 노름을 유도하여 물리적 척도를 반영하고 결합 시스템의 잘 정의됨을 보장하기.
• MINRES 반복 해법 내에서 조건부를 적용하여 메쉬 수렴성과 매개변수에 강건한 수렴성 확보하기.
• 다양한 메쉬 크기와 물리적 매개변수를 가진 2D Darcy-Stokes 및 유체-구조 상호작용 문제에 대한 수치 실험을 통해 방법을 검증하기.

실험 결과

연구 질문

• RQ1 점성도 및 투과도와 같은 물리적 매개변수의 변화에 대해 강건한 단일화된 조건부 프레임워크를 결합 다물리컬 문제에 적용할 수 있는가?
• RQ2 가중 및 분수 스오볼레프 공간을 어떻게 활용하여 인터페이스 변수의 올바른 척도를 반영하는 조건부를 설계할 수 있는가?
• RQ3 Darcy-Stokes 및 유체-구조 상호작용 문제와 같은 경우에 메쉬 정밀도 증가 조건에서도 제안된 조건부가 조건수를 유한하게 유지하고 최적의 수렴 속도를 확보하는가?
• RQ4 해결 시간과 반복 횟수 측면에서 단일화된 조건부는 도메인 분할 방법 대비 상대적으로 어떤 성능을 보이는가?
• RQ5 적절한 가정 하에 하위문제의 잘 정의됨이 고전적으로 확립되지 않은 문제에도 이 프레임워크를 적용할 수 있는가?

주요 결과

• 제안된 조건부는 메쉬 정밀도에 관계없이 조건수를 유한하게 유지하며, 예를 들어 Stokes-Navier 문제에서 h = 2⁻⁸일 때 약 12.72의 조건수를 기록하고 최적의 수렴 속도를 확보하였다.
• 수치 실험 결과, 도메인 분할 방법 대비 단일화된 조건부가 우수한 성능을 보였으며, Darcy-Stokes 문제의 경우 반복당 약 47초, Stokes-Navier 문제의 경우 약 116초의 해결 시간을 기록했고, 유사한 성능을 보이는 도메인 분할 방법 대비 60~150초의 시간을 소요하였다.
• Beavers-Joseph 계수의 극단적인 값(k = 10⁻⁶에서 1까지) 변화에 대해서도 수렴성이 안정적으로 유지되며, 매개변수에 강건함을 입증하였다.
• MINRES 솔버의 수렴은 메쉬에 독립적이며 물리적 매개변수에 대해 강건하며, 예를 들어 Stokes-Navier 문제에서 h = 2⁻⁸일 때 반복 횟수가 일관되게 47회로 유지되었다.
• 분석 [24]에서 유도된 이론적 노름이 라그랑주 승수를 포함하도록 성공적으로 확장되어 강건한 조건부를 가능하게 하였다.
• Darcy-Stokes 및 유체-구조 상호작용 문제와 같은 두 가지 벤치마크 문제에서 프레임워크가 검증되었으며, 오차 수렴 속도는 이론적 추정치와 일치하였다.