[논문 리뷰] A second-order time-stepping scheme for simulating ensembles of parameterized flow problems

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

Numerical methods for differential equations참고 문헌 7인용 수 3

한 줄 요약

이 논문은 다양한 물리적 파라미터, 초기 조건, 경계 조건, 외력 함수를 가진 다수의 비압축성 스토크스 유동 문제를 효율적으로 시뮬레이션하기 위해 이차 정확도를 가지며, 앙상블 기반의 시간 스텝 방법을 제안한다. 반응형, 앙상블 평균 형식을 사용함으로써, 각 시간 스텝에서 다수의 우변을 가진 단일 선형 시스템을 풀게 되어 계산 비용을 크게 감소시키면서도 이차 정확도를 유지하고 안정성 및 오차 추정을 입증한다.

ABSTRACT

We consider settings for which one needs to perform multiple flow simulations based on the Navier-Stokes equations, each having different values for the physical parameters and/or different initial condition data, boundary conditions data, and/or forcing functions. For such settings, we propose a second-order time accurate ensemble-based method that to simulate the whole set of solutions, requires, at each time step, the solution of only a single linear system with multiple right-hand-side vectors. Rigorous analyses are given proving the conditional stability and error estimates for the proposed algorithm. Numerical experiments are provided that illustrate the analyses.

연구 동기 및 목표

• 다양한 물리적 파라미터, 초기 조건, 경계 조건, 외력 함수를 가진 다수의 유동 시뮬레이션에서 발생하는 높은 계산 비용을 해결하기 위해.
• 다양한 유동 앙상블에서 효율성을 유지하면서도 시간 정확도와 이차 정확도를 갖춘 수치적 방법을 개발하기 위해.
• 제안된 앙상블 방법의 안정성과 엄밀한 오차 추정을 확보하기 위해.
• 데이터 융합, 저차원 모델링, 민감도 분석과 같은 응용 분야를 위한 대규모 효율적 시뮬레이션을 가능하게 하기 위해.

제안 방법

• 속도와 압력 갱신에 대해 삼단계 역차분 공식을 사용한 이차 정확도를 가지는 반응형 시간 이산화를 적용한다.
• 앙상블 평균 속도 $ \overline{u}^n $ 과 공통의 점성도 $ \overline{\nu} $ 를 도입하여, 모든 앙상블 구성원에서 동일한 계수 행렬을 사용할 수 있도록 한다.
• 비선형 대류 항을 앙상블 평균 및 변동 항의 조합을 통해 분리함으로써 일致성과 안정성을 확보한다.
• 이산 시스템은 $ J $ 개의 우변을 가진 하나의 선형 시스템으로 해석되며, 여기서 $ J $ 는 앙상블 구성원의 수이다. 이는 효율적인 LU 분해 또는 블록 반복 해법을 가능하게 한다.
• 각 시간 스텝에서 압축성 조건 $ \nabla \cdot u_j^{n+1} = 0 $ 를 유지한다.
• 에너지 추정과 이산 그론월라 부등식을 사용한 안정성 및 오차 분석을 통해 $ L^2 $ 및 $ H^1 $ 노름에서 최적의 수렴 속도를 도출한다.

실험 결과

연구 질문

• RQ1 다양한 파라미터와 입력 조건을 가진 다수의 스토크스 유동 시뮬레이션에 대해 높은 효율성을 유지하면서도 이차 정확도를 갖는 시간 정확도 앙상블 방법을 구성할 수 있는가?
• RQ2 다수의 우변을 가진 단일 선형 시스템을 사용하여 정확도나 안정성을 희생시키지 않고 다수의 유동 앙상블을 시뮬레이션할 수 있는가?
• RQ3 이차 수렴을 보장하기 위해 앙상블 평균과 개별 구성원 보정 사이의 최적 균형은 무엇인가?
• RQ4 다양한 점성도 및 외력 함수 조건 하에서 제안된 방법의 이론적 오차 한계와 안정성 조건은 무엇인가?
• RQ5 각 앙상블 구성원에 대해 별도의 시스템을 풀는 것과 비교해 볼 때, 이 방법의 계산 비용은 어떻게 스케일링되는가?

주요 결과

• 적절한 정규성 조건 하에서, 속도 및 압력 근사에 대해 시간 기반 이차 정확도를 달성한다.
• 앙상블 점성도와 해의 정규성에 따라 제한되는 시간 스텝 제약 조건 하에서 조건부 안정성을 확보한다.
• 오차 추정은 최적의 수렴 속도를 보이며, 속도의 경우 $ H^1 $-노름에서 $ \mathcal{O}(h^{k+1}) $, 압력의 경우 $ L^2 $-노름에서 $ \mathcal{O}(h^{s+1}) $ 이다. 여기서 $ k $ 와 $ s $ 는 유한요소 다항식 차수와 관련된다.
• 각 시간 스텝에서 $ J $ 개의 우변을 가진 단일 선형 시스템을 풀기 때문에, $ J $ 개의 별도 시스템을 풀어야 하는 것보다 계산 비용이 감소한다.
• 오차 한계에는 $ |\nu_j - \overline{\nu}| $ 에 의존하는 항이 포함되어 있어 점성도 변화에 민감함을 나타내지만, 이는 $ \overline{\nu}^{-1} $ 으로 유한하게 제한되어 있어 작은 상대적 차이에 대해서는 안정성을 보장한다.
• 수치 실험은 이론적 오차 추정을 확인하고, 특히 $ J $ 가 클 경우 실용적 효율성 향상을 입증한다.