[논문 리뷰] Analysis of a low complexity, time-accurate discretization of the Navier-Stokes equations

Victor DeCaria, William Layton외 1명|arXiv (Cornell University)|2018. 10. 15.

Numerical Methods and Algorithms참고 문헌 1인용 수 2

한 줄 요약

이 논문은 표준 후방 오일러 이산화에 두 단계 선형 시간 필터를 적용하여 나이비에-스토크스 방정식에 대한 저복잡도이자 시간 정확도를 갖는 방법을 제안한다. 이 방법은 계산 비용을 증가시키지 않으면서 정확도를 향상시키고 과다 감쇠를 제거하며, 무조건적인 에너지 안정성을 유지한다. 수치 실험에서는 항력과 양력 예측이 향상된다.

ABSTRACT

This report presents a low computational and cognitive complexity, stable, time accurate and adaptive method for the Navier-Stokes equations. The improved method requires a minimally intrusive modification to an existing program based on the fully implicit / backward Euler time discretization, does not add to the computational complexity, and is conceptually simple. The backward Euler approximation is simply post-processed with a two-step, linear time filter. The time filter additionally removes the overdamping of Backward Euler while remaining unconditionally energy stable, proven herein. Even for constant stepsizes, the method does not reduce to a standard / named time stepping method but is related to a known 2-parameter family of A-stable, two step, second order methods. Numerical tests confirm the predicted convergence rates and the improved predictions of flow quantities such as drag and lift.

연구 동기 및 목표

• 나이비에-스토크스 방정식에 대해 시간 정확도, 안정성 및 계산 효율성이 뛰어난 시간 이산화 방법을 개발하기 위해.
• 기존의 후방 오일러 솔버를 수정할 때 인지적 및 계산적 오버헤드를 최소화하기 위해.
• 후방 오일러에 내재된 과다 감쇠를 제거하면서도 무조건적인 에너지 안정성을 유지하기 위해.
• 항력과 양력과 같은 주요 유동 양의 예측 정확도를 향상시키기 위해.
• 구현 시 기존 코드베이스에 최소한의 영향을 주며 단순성을 유지하기 위해.

제안 방법

• 두 단계 선형 시간 필터를 후방 오일러 해의 후처리 단계로 적용한다.
• 필터는 수치 감쇠를 줄이면서도 안정성을 유지하도록 시간 이산화를 수정한다.
• 완전히 암시적인 후방 오일러 방법의 변형으로 유도되며, 시스템 행렬 조립에 대한 변경이 필요 없다.
• 필터는 논문에서 증명된 바와 같이 무조건적인 에너지 안정성을 갖도록 설계된다.
• 이 방법은 두 개의 매개변수를 갖는 알려진 A-안정성, 두 번째 차수, 두 단계 방법의 가문과 관련이 있다.
• 최소한의 코드 변경으로 구현 가능하여 기존의 레거시 솔버에 통합하기에 적합하다.

실험 결과

연구 질문

• RQ1 간단한 후처리 필터가 계산 비용 증가 없이 후방 오일러의 정확도를 향상시킬 수 있는가?
• RQ2 필터링된 방법이 수치 과다 감쇠를 줄이면서도 무조건적인 에너지 안정성을 유지하는가?
• RQ3 이 방법이 흐름 시뮬레이션에서 항력과 양력 예측을 얼마나 향상시키는가?
• RQ4 수렴성과 안정성 측면에서 표준 두 번째 차수 시간 적분기와 비교해 볼 때 어떻게 성능을 내는가?
• RQ5 기존의 후방 오일러 코드베이스에 최소한의 변경으로 이 방법을 구현할 수 있는가?

주요 결과

• 필터링된 방법은 비록 표준 이름이 부여된 방법으로 내재되지 않더라도, 일정한 시간 간격을 사용함에도 불구하고 시간 기준으로 두 번째 차수 정확도를 달성한다.
• 수치 실험은 시간 이산화의 예측된 수렴률을 확인한다.
• 표준 후방 오일러에 비해 계산된 항력과 양력 계수의 정확도가 크게 향상된다.
• 시간 필터는 후방 오일러에 존재하는 과도한 수치 감쇠를 성공적으로 제거하면서도 무조건적인 에너지 안정성을 유지한다.
• 논문에서 엄밀하게 증명된 바와 같이 방법은 여전히 무조건적인 에너지 안정성을 유지한다.
• 계산 복잡도의 증가 없이도 최소한의 코드 기여로 구현 가능하다.