[논문 리뷰] The vortex method for 2D ideal flows in exterior domains

Diogo Arsénio, Emmanuel Dormy외 1명|arXiv (Cornell University)|2017. 07. 05.

Navier-Stokes equation solutions참고 문헌 24인용 수 4

한 줄 요약

이 논문은 부드러운 경계를 가진 외부 영역에서 2차원 이상성 유체 흐름을 시뮬레이션하기 위한 복소법의 엄밀한 정당성을 제시한다. 여기서 코어시티는 경계에 따라 점 소용리로 근사되어 접선 조건과 순환량 보존을 확보한다. 복소법 근사의 수렴성을 입증하고, 수치적 조건을 개선하며 안정성을 높이는 유체 전하 방법을 도입한다.

ABSTRACT

The vortex method is a common numerical and theoretical approach used to implement the motion of an ideal flow, in which the vorticity is approximated by a sum of point vortices, so that the Euler equations read as a system of ordinary differential equations. Such a method is well justified in the full plane, thanks to the explicit representation formulas of Biot and Savart. In an exterior domain, we also replace the impermeable boundary by a collection of point vortices generating the circulation around the obstacle. The density of these point vortices is chosen in order that the flow remains tangent at midpoints between adjacent vortices and that the total vorticity around the obstacle is conserved. In this work, we provide a rigorous justification of this method for any smooth exterior domain, one of the main mathematical difficulties being that the Biot--Savart kernel defines a singular integral operator when restricted to a curve (here, the boundary of the domain). We also introduce an alternative method---the fluid charge method---which, as we argue, is better conditioned and therefore leads to significant numerical improvements.

연구 동기 및 목표

• 부드러운 외부 영역에서 2차원 이상성 흐름에 대한 복소법의 엄밀한 수학적 정당성을 제공한다. 이 경우 경계는 이산적인 점 소용리로 표현된다.
• 경계 곡선에 제한된 비에-사바르 커널에 의해 발생하는 특이 적분 연산자 문제를 다룬다.
• 복소법 근사가 옐러 방정식의 진정한 해로 향해 동적 및 정적 수렴성을 확립한다.
• 수치적 조건과 안정성을 향상시키는 대안 방법인 유체 전하 방법을 도입하고 분석한다.
• 특히 리만 합 이산화의 맥락에서 복소법과 유체 전하 근사의 수렴 속도를 증명한다.

제안 방법

• 경계 곡선에 점 소용리를 놓아 외부 영역의 코어시티를 유한한 합으로 근사함으로써 접선 흐름과 총 순환량 보존을 확보한다.
• 비에-사바르 법칙을 사용하여 코어시티 분포로부터 속도장을 재구성한다. 이 경우 경계는 밀도가 유체의 비율 조건을 만족하도록 선택된 소용리판으로 간주된다.
• 경계 적분 표현에 기반한 속도의 특이 적분 연산자를 분석하기 위해 플레멜지 및 푸앵카레-버트랑 공식을 적용한다.
• 균일하거나 준균일 메쉬를 사용하여 경계 적분 방정식을 이산화하고, 이산 연산자의 유계성과 스펙트럼 수렴성을 증명한다.
• 복소법의 대안으로 유체 전하 방법을 도입한다. 이 경우 경계는 점 소용리가 아니라 유체 전하(이중층 포텐셜)로 표현되며, 더 잘 조절된 선형 시스템을 이끈다.
• 리만 합의 오차를 정량화하기 위해 스펙트럼 분석과 리만 합의 수렴 추정치를 사용한다. 특히 경계에서의 매끄럽고 주기적인 함수에 대해 고려한다.

실험 결과

연구 질문

• RQ1 비에-사바르 커널이 경계에서 특이성을 띠는 바에다가도, 부드러운 외부 영역에서 2차원 이상성 흐름에 대한 복소법이 엄밀히 정당화될 수 있는가?
• RQ2 이러한 영역에서 복소법 근사가 옐러 방정식의 진정한 해로 수렴하는 속도는 어떠한가?
• RQ3 유체 전하 방법은 전통적인 복소법에 비해 수치적 조건과 안정성 측면에서 어떻게 비교되는가?
• RQ4 이산화된 경계 적분 연산자의 스펙트럼 특성은 무엇이며, 이는 수치적 방법의 수렴성에 어떻게 영향을 미치는가?
• RQ5 경계에서 매끄럽고 주기적인 함수에 대해 사용되는 리만 합의 수렴성이 최적의 속도로 정량화될 수 있는가?

주요 결과

• 복소법은 외부 영역에서 2차원 옐러 방정식의 진정한 해로 균일하게 수렴하며, 코어시티의 정(regularity)과 메쉬 품질에 따라 오차 한계 $ \mathcal{O}(N^{-\beta}) $를 가진다. 여기서 $ \beta > 0 $이다.
• 유체 전하 방법은 복소법에 비해 더 잘 조절된 선형 시스템을 제공하여 시뮬레이션에서 수치적 안정성과 성능을 크게 향상시킨다.
• 이산화된 경계 적분의 리만 합 수렴은 적분함수의 매끄러움 조건 하에 $ \mathcal{O}(N^{-\beta}) $의 속도로 정량화되며, 여기서 $ \beta \geqslant 2 $이다.
• 복소법에서 유도된 이산 시스템은 Poincaré–Bertrand 항등식의 이산 형태를 만족하며, 이는 일致성과 안정성 증명에 필수적이다.
• 적절한 메쉬 조건 하에서 복소법과 관련된 이산 연산자의 스펙트럼은 0으로부터 멀리 떨어져 있으므로 선형 시스템의 잘 정의됨이 보장된다.
• 유체 전하 근사의 동적 수렴성이 $ L^∞ $ 노름에서 컴acts 시간 간격에 대해 $ \mathcal{O}(N^{-\beta}) $의 수렴 속도로 확립된다. 여기서 $ \beta > 0 $이다.