[논문 리뷰] An efficient algorithm for simulating ensembles of parameterized flow problems

Max Gunzburger, Nan Jiang외 1명|arXiv (Cornell University)|2017. 05. 25.

Fluid Dynamics and Turbulent Flows참고 문헌 16인용 수 3

한 줄 요약

이 논문은 다양한 점성 계수, 초기 조건, 외력 조건을 가진 다수의 매개변수화된 비압축성 스토크스 유동을 시뮬레이션하기 위한 새로운 반암시적 앙상블 알고리즘을 제안한다. 공통 계수 행렬을 사용하여 다중 우변을 가진 단일 선형 시스템을 구성함으로써, 시간에 대한 일阶 정확도와 최적의 공간 수렴성을 유지하면서도 상당한 계산 비용 절감을 달성한다. 이는 시간 간격과 매개변수 편차 제약 조건 하에서 엄밀한 안정성 및 오차 추정을 포함한다.

ABSTRACT

Many applications of computational fluid dynamics require multiple simulations of a flow under different input conditions. In this paper, a numerical algorithm is developed to efficiently determine a set of such simulations in which the individually independent members of the set are subject to different viscosity coefficients, initial conditions, and/or body forces. The proposed scheme applied to the flow ensemble leads to need to solve a single linear system with multiple right-hand sides, and thus is computationally more efficient than solving for all the simulations separately. We show that the scheme is nonlinearly and long-term stable under certain conditions on the time-step size and a parameter deviation ratio. Rigorous numerical error estimate shows the scheme is of first-order accuracy in time and optimally accurate in space. Several numerical experiments are presented to illustrate the theoretical results.

연구 동기 및 목표

• 다른 점성 계수, 초기 조건, 외력 조건을 가진 매개변수화된 비압축성 유동 문제의 앙상블 시뮬레이션을 위한 효율적인 수치 알고리즘을 개발한다.
• 각 앙상블 구성원에 대해 별도의 시스템을 풀지 않고, 다중 우변을 가진 단일 선형 시스템을 풀어 계산 비용을 절감한다.
• 시간 간격과 매개변수 편차 비율 제약 조건 하에서 비선형적이고 장기적인 안정성을 확립한다.
• 포괄적인 오차 추정을 통해 시간에 대한 일계 정확도와 공간에 대한 최적 정확도를 증명한다.
• 기준 흐름 문제에 대한 수치 실험을 통해 이론적 결과를 검증한다.

제안 방법

• 비선형 대류 항은 명시적으로, 점성 항과 압력 항은 암시적으로 처리하는 반암시적 시간 이산화를 스토크스 방정식에 적용한다.
• 다른 점성 계수, 초기 조건, 또는 외력 함수에 관계없이 모든 앙상블 구성원에서 동일한 계수 행렬을 유지하는 새로운 공간 이산화 방법을 도입한다.
• 각 시간 단계에서 유도된 선형 시스템은 단일 계수 행렬과 J개의 우변을 가지며, 각 앙상블 구성원에 대해 하나씩 구성되며, 이는 블록 솔버를 통한 효율적 해법을 가능하게 한다.
• 이산 에너지 방법을 활용한 안정성 분석을 수행하며, 오차 전파를 통제하기 위해 수정된 이산 그론월라 부등식을 도입한다.
• 다양한 노름, 특히 $L^2$, $H^1$, 및 $L^\rho$ 유형의 노름에서 오차 추정을 유도하며, 이는 점성 계수 편차, 시간 간격 크기, 메쉬 크기에 따라 달라진다.
• 변동하는 점성 계수를 처리하기 위해 기준 점성 계수 $\overline{\nu}$ 를 도입하고, 안정성 및 오차 추정에서 $|\nu_j - \overline{\nu}|$ 의 편차를 분석한다.

실험 결과

연구 질문

• RQ1 다른 점성 계수, 초기 조건, 외력 조건을 가진 스토크스 유동의 앙상블 시뮬레이션을 위해 다중 우변을 가진 단일 선형 시스템을 구성할 수 있는가?
• RQ2 시간 간격 크기와 점성 계수 편차 비율에 어떤 조건이 충족되어야 비선형적이고 장기적인 안정성이 확보되는가?
• RQ3 제안된 앙상블 알고리즘의 시간 및 공간 수렴 속도는 얼마이며, 이는 최적인가?
• RQ4 오차 추정은 점성 계수 편차, 시간 간격, 메쉬 크기와 어떻게 스케일링되는가?
• RQ5 각 시뮬레이션을 별도로 풀는 것과 비교해 제안된 방법이 상당한 계산 비용 절감을 이룰 수 있는가?

주요 결과

• 제안된 알고리즘은 시간에 대해 일계 정확도를 확보하고, 유속 및 압력 근사에 대해 공간에서 최적의 $O(h^k)$ 수렴성을 보인다. 여기서 $k$ 는 유한요소 공간의 다항식 차수이다.
• 시간 간격 크기와 점성 계수 편차 비율이 $\overline{\nu}$ 및 $|\nu_j - \overline{\nu}|$ 를 포함한 특정 bound를 만족할 경우, 이 방법은 조건부로 비선형적으로 안정되고 장기적으로 안정하다.
• 오차 추정은 $O(\Delta t^2 + h^{2k})$ 의 스케일을 가지며, 이는 $\overline{\nu}^{-1}$, $\|\nabla u_{j,t}\|_{L^2}$, 및 해의 고차 미분항에 의존하는 상수를 포함한다.
• 계산 비용은 각 시간 단계에서 $J$ 개의 별도 선형 시스템을 풀던 것을 단 한 번의 시스템을 $J$ 개의 우변으로 풀도록 줄여 효율성을 크게 향상시킨다.
• 수치 실험은 이론적 오차 추정을 확인하며, 다양한 매개변수 변화에 걸쳐 방법의 강건성과 효율성을 입증한다.
• 점성 계수가 앙상블 구성원 간에 크게 다를 경우에도, 편차 비율 조건을 충족하면 안정성과 정확도를 유지한다.