[논문 리뷰] Generalised Sobolev Stable Flux Reconstruction

Will Trojak|arXiv (Cornell University)|2018. 04. 04.

Computational Fluid Dynamics and Aerodynamics참고 문헌 32인용 수 5

한 줄 요약

이 논문은 기존의 OSFR 및 ESFR 방법을 확장하여 수정된 소볼레프 노름을 통합한 고차수 플럭스 재구성(FLUX RECONSTRUCTION, FR) 기반의 일반화된 소볼레프 안정성 흐름 재구성(GSFR)을 소개한다. 이 방법은 향상된 시간 안정성과 감소된 소산을 가능하게 하며, 분석적 및 수치적 검증을 통해 선형 및 비선형 편미분방정식(PDE), 특히 유동역학 방정식인 오일러 방정식과 소용루르 운반 문제에서 더 높은 CFL 한계와 향상된 성능을 보여준다.

ABSTRACT

A new set of symmetric correction functions is presented for high-order flux reconstruction, that expands upon, while incorporating, all previous correction function sets and opens the possibility for improved performance. By considering FR applied to the linear advection equation, and through modification to the Sobolev norm, criteria are presented for a wider set of correction functions. Legendre polynomials are then used to fulfil these criterion and realise functions for third to fifth order FR. The sufficient conditions for the existence of the modified norms are also explored, before Fourier and Von Neumann analysis are applied to analytically find temporal stability limits for various Runge-Kutta temporal integration schemes. For all cases, correction functions are found that extend the temporal stability of FR. Two application-inspired investigations are performed that aim to explore the effect of aliasing and non-linear equations. In both cases unique correction functions could be found that give good performance, compared to previous FR schemes, while also improving upon the temporal stability limit.

연구 동기 및 목표

• 플럭스 재구성(FR)에서 안정적인 보정 함수를 위한 일반화된 프레임워크를 개발하여 기존의 OSFR 및 ESFR 기법을 초월한다.
• 에너지 안정성 분석에 사용되는 소볼레프 노름을 수정하여 FR 기법의 시간 안정성을 향상시킨다.
• 선형 및 비선형 PDE의 고차수 시뮬레이션에서 수치적 소산과 앨리아싱 오차를 감소시킨다.
• 선형 운반 문제와 오일러 방정식을 포함한 벤치마크 문제에 대해 새로운 보정 함수를 검증한다.
• 기존의 FR 기법에 비해 GSFR가 더 높은 CFL 한계를 달성하고 소용루르 시뮬레이션에서 운동에너지의 보존을 향상시킬 수 있음을 보여준다.

제안 방법

• 손상된 소볼레프 공간에서 에너지 안정성을 확보하기 위해 수정된 소볼레프 노름을 사용한 일반화된 보정 함수 프레임워크를 유도한다.
• 3차에서 5차 FR 기법의 안정성 기준을 만족시키기 위해 레전드르 다항식을 적용한다.
• 푸리에 및 본 네움 분석을 통해 다양한 룬게쿠타 기법의 시간 안정성 한계를 분석적으로 도출한다.
• 텐서 곱 확장 기법을 사용하여 1차원 FR 수식을 3차원 헥사에드럴 요소를 기반으로 일반화한다.
• 비선형 테스트 케이스에서 수치적 플럭스 계산을 위해 루사노프 플럭스와 특성 속도를 활용한다.
• 이른바 비점성 소용루르(ICV) 테스트 케이스와 선형 비균질 PDE를 사용하여 앨리아싱과 소산을 평가한다.

실험 결과

연구 질문

• RQ1 기존의 OSFR 및 ESFR 기법을 포함하고 개선하는 일반화된 보정 함수 세트를 도출할 수 있는가?
• RQ2 에너지 안정성을 확보하고 새로운 보정 함수를 가능하게 하기 위해 소볼레프 노름에 어떤 수정이 필요한가?
• RQ3 새로운 GSFR 프레임워크는 룬게쿠타 기법에 대해 더 높은 시간 안정성 한계(CFL 수치)를 제공하는가?
• RQ4 OSFR 및 ESFR에 비해 GSFR는 선형 및 비선형 PDE에서 소산과 앨리아싱 오차를 어떻게 감소시키는가?
• RQ5 오일러 방정식과 소용루르 운반 문제와 같은 복잡한 비선형 문제에서 GSFR는 안정성과 정확도를 유지할 수 있는가?

주요 결과

• p = 3일 때, GSFR는 RK44 시간 적분에서 OSFR보다 더 높은 CFL 한계를 확보하며, 분석적 안정성 분석을 통해 향상된 시간 안정성이 확인된다.
• 도표적 교차 분석(그림 11)을 통해 GSFR 보정 함수는 OSFR 또는 ESFR 집합에 포함되지 않으며, 고유하다.
• ICV 테스트 케이스에서 GSFR는 OSFR에 비해 운동에너지 감쇠가 감소한 것으로 나타났지만, 저주파수 영역에서 더 낮은 소산을 보이는 DG가 가장 우수한 성능을 보였다.
• GSFR 기법은 다항식 차수에 관계없이 안정성을 유지하고 수치적 순서 정확도(OOA) 복구를 달성한다.
• 분산 및 소산도 그래프(그림 6)를 통해 GSFR 보정 함수는 고주파수 영역에서 더 낮은 소산을 보이며, 고차수 시뮬레이션에서 향상된 성능을 발휘한다.
• 선형 비균질 PDE 테스트에서 GSFR는 높은 CFL 한계와 낮은 오차를 가진 안정 영역을 확보하여 앨리아싱 효과에 대한 강건성을 보였다.