[논문 리뷰] A new embedded variable stepsize, variable order family of low computational complexity

Victor DeCaria, Ahmet Guzel외 2명|arXiv (Cornell University)|2018. 10. 15.

Numerical methods for differential equations참고 문헌 31인용 수 5

한 줄 요약

이 논문은 ODE용 임베디드 변수단계크기·변수차수(VSVO) 방법인 MOOSE234를 소개한다. 이는 BDF3와 유사한 계산 비용으로 2, 3, 4차 정확도를 달성한다. BDF3 해에 선형 시간 필터를 적용하여 고차수 근사와 오차 추정기를 낮은 추가 비용으로 생성함으로써, 레거시 CFD 코드에 비침습적이고 병렬화 가능한 효율적인 구현을 가능하게 한다.

ABSTRACT

Variable Stepsize Variable Order (VSVO) methods are the methods of choice to efficiently solve a wide range of ODEs with minimal work and assured accuracy. However, VSVO methods have limited impact in timestepping methods in complex applications due to their computational complexity and the difficulty to implement them in legacy code. We introduce a family of implicit, embedded, VSVO methods that require only one BDF solve at each time step followed by adding linear combinations of the solution at previous time levels. In particular, we construct implicit and linearly implicit VSVO methods of orders two, three and four with the same computational complexity as variable stepsize BDF3. The choice of changing the order of the method is simple and does not require additional solves of linear or nonlinear systems.

연구 동기 및 목표

• 유체역학과 같은 복잡한 응용 분야에서 VSVO 방법의 보편적 채택이 제한된 이유인 높은 계산 복잡도와 구현 난이도를 해결하기 위해.
• 기존 VSVO 방법에서 타임스텝당 다수의 비선형 해법이 요구되는 장벽을 최소한의 추가 계산 비용으로 극복하기 위해.
• 고차수 정확도를 추가 해법 없이 분리함으로써, 기존 CFD 코드에 비침습적이고 병렬화 우호적인 통합을 가능하게 하기 위해.
• 단일 BDF3 해를 기반으로 후처리 필터를 적용하여 오차 추정기와 고차수 해를 생성하는 임베디드 방법의 가족을 개발하기 위해.
• 2차 변형에서 A-안정성과 G-안정성을 확보하면서도 변수단계크기 및 변수차수 선택에 대해 낮은 계산 복잡도를 유지하기 위해.

제안 방법

• BDF3 해에 선형 시간 필터를 적용하여 4차 정확도를 달성하는 필터링된 BDF3 방법(FBDF4)을 구성한다.
• BDF3 해에 두 번째 필터인 BDF3-Stab를 적용하여 2차 정확도를 가지며 A-안정성과 G-안정성을 확보한 근사값을 생성한다.
• 필터링된 해의 차이(예: $y^2_{n+1} - y^3_{n+1}$)를 사용하여 적응형 타임스텝을 위한 임베디드 오차 추정기로 활용한다.
• 일단계법(OLM)으로 방법을 제작하여 변수단계크기 전략과의 호환성 및 유리한 안정성 특성을 확보한다.
• 과거 해값의 선형 조합을 사용하여 필터를 설계함으로써 낮은 계산 비용과 SIMD 병렬 처리를 유지한다.
• 섹션 4에서 단일 BDF3 해를 기반으로 임베디드 오차 추정기를 사용해 동적으로 단계크기와 차수를 선택하는 변수단계크기 알고리즘을 구현한다.

실험 결과

연구 질문

• RQ1 단일 BDF3 해에서 낮은 추가 비용으로 고차수, 임베디드, 안정성 확보가 가능한 시간 필터링 접근법을 사용할 수 있는가?
• RQ2 백워드 오일러에서 출발할 때 선형 시간 필터링을 통해 달성 가능한 최대 차수는 얼마이며, 이는 설계 접근법에 어떤 제약을 끼치는가?
• RQ3 BDF3와 유사한 계산 복잡도를 유지하면서도 정확도와 안정성을 확보하는 변수단계크기·변수차수(VSVO) 알고리즘을 구성할 수 있는가?
• RQ4 비선형 시스템의 추가 해법 없이도 오차 추정기를 해법 과정에 임베디드할 수 있는가?
• RQ5 유사한 코드 변경 최소화와 병렬 아키텍처 지원을 통해 레거시 CFD 코드에 효율적으로 통합할 수 있는가?

주요 결과

• 해당 방법은 타임스텝당 단일 BDF3 해를 사용하여 2, 3, 4차 정확도를 달성하며, 고차수 해와 오차 추정기는 선형 시간 필터를 통해 얻는다.
• 2차 변형(BDF3-Stab)은 A-안정성과 G-안정성이 증명되어 있으며, 굳은 문제에 대한 강건성을 보장한다.
• 전체 임베디드 방법의 계산 비용은 단일 BDF3 해에 의해 지배되며, 필터링 단계는 무시할 만큼 낮은 오버헤드를 갖는다.
• 오차 추정기 $Est_2 = y^2_{n+1} - y^3_{n+1}$ 및 $Est_3 = y^4_{n+1} - y^3_{n+1}$ 는 임베디드되어 있으며 추가 기능 평가가 필요하지 않다.
• MOOSE234는 반복적 비선형 해법 없이 바나 더 폴 문제에서 우수한 성능을 보이며, 압축성 없는 나비에-스토크스 방정식의 선형화된 암시적 형태를 성공적으로 해결한다.
• 이 방법은 비침습적이고 SIMD 우아하며 병렬화 가능한 방식으로 레거시 CFD 코드에 통합 가능하며, VSVO 방법의 보급을 방해하는 주요 장벽을 극복한다.