[논문 리뷰] Dispersion-dissipation condition for finite difference schemes

Xiao Hu, V. K. Tritschler외 2명|arXiv (Cornell University)|2012. 04. 23.

Computational Fluid Dynamics and Aerodynamics참고 문헌 10인용 수 26

한 줄 요약

이 논문은 유한차분 스킴에서 잡음이 섞인 고파르수수파를 감쇠시키기 위해 필요한 최소 수치적 소산을 결정하는 분산-소산 조건을 제안한다. 수정된 파수 분석을 통해 군속도와 소산률을 연결함으로써, 선형 및 비선형 스킴에서 분산과 소산을 최적화한다. 이 조건은 WENO-CU6-M2 스킴에 적용되어 난류 흐름의 정밀한 미해결 해석 시뮬레이션, 특히 자기유사 에너지 감쇠와 난류로의 전이를 향상시켰다.

ABSTRACT

A general dispersion-dissipation condition for finite difference schemes is derived by analyzing the numerical dispersion and dissipation of explicit finite-difference schemes. The proper dissipation required to damp spurious high-wavenumber waves in the solution is determined from a physically motivated relation between group velocity and dissipation rate. The application to a previously developed low-dissipation weighted essentially non-oscillatory scheme (WENO-CU6-M2) [X. Y. Hu and N. A. Adams; Scale separation for implicit large eddy simulation, J. Comput. Phys. 230 (2011) 7240-7249] demonstrates that this condition can serve as general guideline for optimizing the dispersion and dissipation of linear and non-linear finite-difference schemes. Moreover, the improved WENO-CU6-M2 which satisfies the dispersion-dissipation condition can be used for under-resolved simulations. We demonstrate this capability by considering transition to turbulence and self-similar energy decay of the three-dimensional Taylor-Green vortex. Simulations of the inviscid and the viscous Taylor-Green vortex at Reynolds numbers ranging from $Re=400$ to $Re=3000$ show a significant improvement over the classical dynamic Smagorinsky model and demonstrate competitiveness with state-of-the-art implicit LES models, while preserving shock-capturing properties.

연구 동기 및 목표

• 유한차분 스킴에서 분산과 소산 간의 정량적 연관성을 확립하여 비물리적 고파르수수파를 제어한다.
• 해결된 척도의 정확도를 떨어뜨리지 않으면서 비물리적 진동을 억제하기 위해 최적의 소산 수준을 결정하는 오랜 도전 과제를 해결한다.
• 특히, 해상도가 낮은 시뮬레이션을 위한 선형 및 비선형 스킴에서 분산과 소산을 최적화하기 위한 일반적 지침을 제공한다.
• WENO-CU6-M2 스킴을 분산-소산 조건을 만족하도록 개선하여, 난류 흐름의 정밀한 대규모 유동 해석을 가능하게 한다.

제안 방법

• 선형 선형성의 이동 방정식을 사용하여 선형 유한차분 스킴의 수정된 파수로부터 분산-소산 조건을 유도한다.
• 위상 오차(|dξ̃R/dξ − 1|)와 소산률(−ξ̃I)의 비율을 정의하여, 경계에 가까운 해상도가 낮은 파동에 대한 충분한 감쇠를 보장한다.
• WENO-CU6-M2 스킴에 이 조건을 적용하기 위해 재구성 절차를 수정하여, 고해상도 특성을 유지하면서 소산을 조절한다.
• 스펙트럼 분석을 통해, 특히 니크비스트 절단 주파수 근처에서 파장 수준의 분산 및 소산 특성을 평가한다.
• Sod의 충격파 튜브, 이중 폭발파 상호작용, Taylor-Green 소용루와 같은 벤치마크 문제를 통해 수정된 스킴을 검증한다.
• 정밀한 직접 시뮬레이션(DNS) 데이터 및 기존의 LES 모델(Dynamic Smagorinsky, ALDM)과의 비교를 통해 에너지 감쇠 및 소산률 예측 정확도를 평가한다.

실험 결과

연구 질문

• RQ1 유한차분 스킴에서 비물리적 고파르수수파를 효과적으로 감쇠시키기 위해 필요한 최소 수치적 소산은 얼마인가?
• RQ2 분산과 소산을 어떻게 균형 있게 조절하여 해결된 척도의 정확도를 유지하면서 비물리적 진동을 억제할 수 있는가?
• RQ3 분산-소산 조건은 WENO-CU6-M2와 같은 비선형 충격 포착 스킴에 적용될 수 있는가?
• RQ4 개선된 WENO-CU6-M2 스킴은 전통적인 LES 모델과 비교해, 난류 흐름의 해상도가 낮은 시뮬레이션에서 어떻게 성능을 발휘하는가?
• RQ5 분산-소산 조건을 만족함으로써, Taylor-Green 소용루에서 자기유사 에너지 감쇠 예측 능력은 어느 정도 향상되는가?

주요 결과

• 분산-소산 조건을 만족하는 수정된 WENO-CU6-M2 스킴은 3차원 점성 없는 Taylor-Green 소용루의 자기유사 에너지 감쇠를 정확히 포착하여 이론적 t−1.2 감쇠율과 일치한다.
• Re=400에서 Re=3000까지의 점성 있는 Taylor-Green 소용루 시뮬레이션에서, 개선된 스킴은 기존의 동적 Smagorinsky 모델보다 소산률 예측에서 뚜렷한 향상을 보였다.
• Re=400 및 Re=800에서, 스킴의 소산률은 최신 기술의 암시적 LES 모델인 적응형 국소 탈구조 모델(ALDM)과 유사한 수준이었다.
• Kolmogorov 관성 영역에서 k−5/3 스케일링이 절단 주파수까지 잘 유지되어, 해결된 척도에서 최소한의 수치적 소산이 발생했음을 시사한다.
• Sod의 충격파 튜브 및 이중 폭발파 문제를 통해, 스킴은 비물리적 파동 전파를 줄이면서도 정확한 충격 포착 성능을 유지함을 확인했다.
• 분산-소산 조건은 선형 및 비선형 스킴에서 수치적 소산을 체계적으로 최적화할 수 있도록 하여, 난류 흐름의 고정밀 시뮬레이션을 위한 일반적인 프레임워크를 제공한다.