[논문 리뷰] High resolution methods for time dependent problems with piecewise smooth solutions

Eitan Tadmor|ArXiv.org|2003. 05. 01.

Computational Fluid Dynamics and Aerodynamics참고 문헌 20인용 수 7

한 줄 요약

이 논문은 조각별로 스무스한 해를 가지는 시간에 의존하는 대류지배 문제를 위한 고해상도 수치 방법을 제시한다. 국소적 방법으로는 중심 스킴을, 전역적 방법으로는 적응형 스펙트럼 점증성(스펙트럼 점증성, SV)을 사용한다. 가장자리 탐지, 적응형 재구성, 스펙트럼 점증성을 통해 지수 수렴을 입증하여 충격파와 불연속성을 진동 없이 정확하게 해석한다.

ABSTRACT

A trademark of nonlinear, time-dependent, convection-dominated problems is the spontaneous formation of non-smooth macro-scale features, like shock discontinuities and non-differentiable kinks, which pose a challenge for high-resolution computations. We overview recent developments of modern computational methods for the approximate solution of such problems. In these computations, one seeks piecewise smooth solutions which are realized by finite dimensional projections. Computational methods in this context can be classified into two main categories, of local and global methods. Local methods are expressed in terms of point-values (-- Hamilton-Jacobi equations), cell averages (-- nonlinear conservation laws), or higher localized moments. Global methods are expressed in terms of global basis functions. High resolution central schemes will be discussed as a prototype example for local methods. The family of central schemes offers high-resolution ``black-box-solvers'' to an impressive range of such nonlinear problems. The main ingredients here are detection of spurious extreme values, non-oscillatory reconstruction in the directions of smoothness, numerical dissipation and quadrature rules. Adaptive spectral viscosity will be discussed as an example for high-resolution global methods. The main ingredients here are detection of edges in spectral data, separation of scales, adaptive reconstruction, and spectral viscosity.

연구 동기 및 목표

• 대류 지배 문제에서 충격파와 꺾임과 같은 비스무스한 특징을 정확하게 계산하는 데 도전하는 것.
• 불연속성이 존재하는 상황에서도 정확성과 안정성을 유지하는 고해상도 수치 방법을 개발하는 것.
• 조각별로 스무스한 성질을 이용하여 이론적 수렴 속도와 관측된 계산 성능 사이의 격차를 메우는 것.
• 스펙트럼 정확도를 유지하면서 엔트로피 소산을 강제하는 적응형 스펙트럼 점증성을 안정화 기법으로 도입하는 것.
• 적응형 모이피어와 필터링을 통해 물리 공간과 푸리에 공간에서의 재구성에서 지수 수렴을 입증하는 것.

제안 방법

• 점 값, 셀 평균, 비진동 재구성 기반의 고해상도 중심 스킴을 사용하여 충격파를 진동 없이 포착한다.
• 스펙트럼 방법의 안정화와 기브스 진동 방지를 위해 필터 $\sigma_p(\xi) = e^{\beta \xi^p / (\xi^2 - 1)}$ 를 사용하는 적응형 스펙트럼 점증성(SV)을 적용한다.
• 조각별로 스무스한 영역에서 지수 수렴을 달성하기 위해 $\theta = d(x)$ 와 $p \sim d(x)N$ 인 적응형 모이피어 $\Psi_{\theta,p}$ 를 사용한다.
• 정밀도 향상을 위해 딜리클레 커널 $D_p(y)$ 와 게브레-클래스 모이피어 $\rho_\beta$ 를 활용하여 국소화 및 상쇄 효과를 확보한다.
• 가장자리 탐지와 스펙트럼 필터링을 조합하여 스무스한 방향에서 해를 재구성함으로써 불연속성 근처의 오차를 최소화한다.
• SV 정규화를 적용한 수정된 갈레르킨 방법을 도입한다: $\partial_t v_N + \mathcal{P}_N \nabla_x \cdot f(v_N) = -N \sum_{|k|\leq N} \sigma(|k|/N) \hat{v}_k \phi_k(x)$, 이는 엔트로피 소산을 보장한다.

실험 결과

연구 질문

• RQ1 대류 지배 문제에서 조각별로 스무스한 해를 정확하게 해석하기 위한 고해상도 수치 방법을 어떻게 설계할 수 있는가?
• RQ2 적응형 스펙트럼 점증성은 스펙트럼 정확도를 훼손하지 않으면서 스펙트럼 방법을 어떻게 안정화하는가?
• RQ3 스펙트럼 또는 푸리에 사영으로부터 조각별로 스무스한 함수의 재구성에서 지수 수렴을 달성할 수 있는가?
• RQ4 모이피어 파rameter $\theta$ 와 $p$ 의 적응적 선택이 불연속성 근처의 해상도를 어떻게 향상시키는가?
• RQ5 불연속성까지의 거리 $d(x)$, 모드 수 $N$, 그리고 스펙트럼 방법의 수렴 속도 사이의 관계는 무엇인가?

주요 결과

• 조각별로 스무스한 함수의 적응형 재구성에서 지수 수렴이 달성된다: $|\Psi_{\{d(x),d(x)N\}}*\mathcal{P}_N v(x) - v(x)| \leq C (d(x)N)^r e^{-\gamma \sqrt{d(x)N}}$.
• 고해상도 스킴의 $L^1$-오차는 수렴 속도 $\|v^h - u\|_{L^1} \leq \alpha_t h |\log h|$ 로서 일반적인 BV 해의 이론적 $O(h^{1/2})$ 에 비해 훨씬 빠르게 수렴한다.
• 적응형 스펙트럼 점증성은 엔트로피 소산을 강제함으로써 스펙트럼 근사의 강한 수렴을 보장하여 임의의 진동을 방지한다.
• 소규모이고 국소적인 점증성만을 고주파 성분에 추가함으로써 형식적 스펙트럼 정확도를 유지하며, $s \in (0, \infty)$ 가 소산 척도를 제어한다.
• 수치 실험은 물리 공간과 푸리에 공간 양쪽에서 지수 수렴 속도를 확인하였으며, 그 결과는 그림 5–8에 나타나 있다.
• 이론적 분석은 비선형 보존법칙에서 $p > 1$ 인 경우 $L^p$-안정성이 실패함을 보여주며, 이는 BV 및 조각별로 스무스한 정규성 프레임워크의 필요성을 정당화한다.