[논문 리뷰] A general approach to enhance slope limiters on non-uniform grids

Xianyi Zeng|arXiv (Cornell University)|2013. 01. 06.

Computational Fluid Dynamics and Aerodynamics참고 문헌 19인용 수 3

한 줄 요약

이 논문은 비균일 격자에서 고해상도 스킴, 특히 MUSCL 방법에 대해 기울기 제한 함수를 향상시키기 위한 일반적인 프레임워크를 제시한다. 이는 제2차 정확도, TVD 안정성, 대칭성 유지에 대한 충분조건을 유도함으로써 이루어지며, 고전적 재구성 이론과 하르텐의 안정성 이론을 비균일 1차원 격자로 확장한다. 1차원 및 2차원 수치 예제를 통해 검증되었으며, 기존의 제한 함수를 사용할 경우 비정규 격자에서 상실되는 제2차 수렴 속도를 복원한다.

ABSTRACT

Abstract. A general approach to study and enhance the slope limiter functions on non-uniform grids is presented. Slope limiters are preferred in high-resolutions schemes in general and MUSCL in particular to solve hyperbolic conservation laws. However, most 1D limiters are developed assuming uniform meshes in space, which are shown to be inadequate on non-uniform grids. Especially, second-order convergence is shown to be lost when the conventional limiters are applied on irregular grids in the case of smooth solutions. A methodology based on the classical reconstruct-evolve-project approach and Harten’s stability theory is presented to study the slope limiters on 1D non-uniform computational grids. Sufficient conditions for the limiters to lead to formal second-order spatial accuracy, total-variational-diminishing stability and symmetry-preserving property are derived. The analysis and results extend naturally to cell-centered finite volume methods in multiple dimensions. Several most widely used conventional limiters are enhanced to satisfy these conditions, and their performances are illustrated by various 1D and 2D numerical examples.

연구 동기 및 목표

• 고해상도 스킴에서 기존의 기울기 제한 함수를 비균일 격자에 적용할 경우 발생하는 제2차 공간 정확도의 손실 문제를 해결하기 위해.
• 재구성-진행-투영 프레임워크와 하르텐의 안정성 이론을 기반으로 한 체계적인 방법론을 개발하여 비균일 1차원 격자에서 기울기 제한 함수를 분석하기 위해.
• 비균일 이산화에서 형식적 제2차 정확도, 총변동감소(Total-variation-diminishing, TVD) 안정성, 대칭성 유지 성질을 보장하는 충분조건을 도출하기 위해.
• 다차원에서의 세포 중심 유한체적 방법으로 분석을 확장하기 위해.
• 유 widely 사용되는 기존의 제한 함수를 도출된 조건을 만족하도록 개선하고, 그 성능을 수치적으로 검증하기 위해.

제안 방법

• 비균일 격자에서 기울기 제한 함수를 분석하기 위한 기초로 재구성-진행-투영 접근법을 채택하기 위해.
• 비균일 메쉬의 맥락에서 TVD 안정성과 대칭성 유지에 대한 충분조건을 도출하기 위해 하르텐의 안정성 이론을 적용하기 위해.
• 비균일 격자에서 형식적 제2차 공간 정확도를 보장하기 위해 제한 함수에 대한 수학적 제약 조건을 설정하기 위해.
• 기존의 제한 함수(예: minmod, van Leer, Superbee)를 수정하여 도출된 조건을 만족하도록 제한 함수를 변형하기 위해.
• 비정규 격자에서 정확도와 안정성을 테스트하기 위해 1차원 및 2차원 유한체적 스킴에 개선된 제한 함수를 구현하기 위해.
• 정밀한 해를 사용하여 수렴 속도를 정량화하고, 제2차 정확도가 복원되었는지 확인하기 위해.

실험 결과

연구 질문

• RQ1 비균일 격자에서 1차원 초구형 보존 법칙에 대해 기울기 제한 함수가 제2차 공간 정확도를 유지하기 위해 만족해야 할 조건은 무엇인가?
• RQ2 하르텐의 안정성 이론은 비균일 계산 격자에서 기울기 제한 함수를 분석하기 위해 어떻게 적응시킬 수 있는가?
• RQ3 비균일 메쉬에서 TVD 안정성과 대칭성 유지 보장을 보장하기 위해 기존의 제한 함수에 어떤 수정이 필요한가?
• RQ4 개선된 제한 함수가 비정규 격자에서 부드러운 해에 대해 제2차 수렴을 얼마나 잘 복원하는가?
• RQ5 제안된 프레임워크는 다차원 세포 중심 유한체적 방법으로 확장될 수 있는가?

주요 결과

• 기존의 기울기 제한 함수는 비균일 격자에 적용되더라도 부드러운 해에 대해서조차 제2차 공간 정확도를 상실한다.
• 유도된 충분조건은 비균일 1차원 격자에서 기울기 제한 함수에 대해 형식적 제2차 정확도, TVD 안정성, 대칭성 유지 보장을 보장한다.
• 일반적으로 사용되는 제한 함수들(예: minmod, van Leer, Superbee)의 개선된 버전은 비균일 격자에서 제2차 수렴을 성공적으로 복원한다.
• 이 방법론은 다차원 세포 중심 유한체적 방법으로 자연스럽게 일반화된다.
• 1차원 및 2차원 수치 예제를 통해 개선된 제한 함수가 비정규 격자에서 고해상도 성질을 유지하면서 최적의 정확도를 달성함을 확인하였다.
• 이 접근법은 실용적인 유체역학 계산 응용 분야에서 비균일 격자에 적합한 강력하고 정확한 제한 함수를 설계하기 위한 체계적인 프레임워크를 제공한다.